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Abstract 

The  focus  of  this  research  was  to  acquire  definitive  experimental  data  on  pre¬ 
dominant  point  defects  in  three  important  wide-band-gap  semiconductors.  Hyperhne 
interactions  in  electron  paramagnetic  resonance  spectra  were  used  to  characterize 
the  neutral  nitrogen  acceptor  in  zinc  oxide,  to  identify  a  silicon  interstitial  impurity 
in  titanium  dioxide,  and  to  determine  the  electronic  structure  of  the  singly  ionized 
sulfur  vacancy  in  stannous  hexathiohypodiphosphate  (SPS).  Research  on  the  basic 
properties  of  these  technologically  important  materials  plays  a  crucial  role  in  the 
development  of  advanced  optical  and  electronic  systems. 

Zinc  oxide  is  an  electro-optic  material  with  the  potential  to  produce  high- 
performance  electronics  and  also  ultraviolet  detectors  and  emitters.  The  angular 
dependencies  of  axial  and  basal  6'Zn  hyperhne  lines  were  used  to  extract  6'Zn  hy¬ 
perhne  interaction  parameters  for  the  neutral  nitrogen  acceptor  in  zinc  oxide.  These 
data  establish  the  exact  nature  of  the  ground-state  wave  function  associated  with  this 
important  acceptor.  The  experimental  information  obtained  in  this  study  about  the 
ground-state  wave  function  will  allow  computational  theorists  to  model  the  nitrogen 
acceptor  and  determine  whether  it  can  successfully  lead  to  p-type  conduction  in  zinc 
oxide.  The  search  for  p-type  zinc  oxide  is  one  of  the  most  important  basic  questions 
in  present-day  research  on  wide-band-gap  semiconductors. 

Titanium  dioxide  is  an  energy-conversion  material  being  developed  for  hydrogen 
generation  and  use  as  an  electrode  in  lithium-ion  batteries.  The  angular  dependencies 
of  an  1  =  0  line  and  two  sets  of  29 Si  hyperhne  lines  were  used  to  calculate  g  param¬ 
eters  and  29Si  hyperhne  interaction  parameters  for  a  silicon-related  S  =  |  defect  in 
rutile  titanium  dioxide.  These  data  were  used  to  identify  the  defect  as  a  deep  donor 
composed  of  two  silicons,  with  one  silicon  at  a  titanium  site  and  the  other  silicon  at  a 
neighboring  tetrahedral  site.  The  model  established  in  this  study  for  silicon  in  rutile 


IV 


titanium  dioxide  is  unique  and  confirms  that  silicon  is  an  ubiquitous  impurity  in  this 
material.  Its  presence  may  have  a  critical  effect  in  photocatalytic  applications. 

SPS  is  a  high-speed  photorefractive  material  that  is  being  developed  for  use 
as  a  protective  element  in  visible  and  infrared  sensor  systems  to  guard  against  laser 
induced  damage  and/or  jamming.  The  angular  dependencies  of  two  sets  of  11 7Sn  and 
119Sn  hyperfine  lines,  and  two  sets  of  31P  hyperfine  lines  were  used  to  determine  g 
parameters  and  hyperfine  interaction  parameters  for  an  intrinsic  S  =  |  defect  in  SPS. 
These  data  were  used  to  identify  the  defect  as  a  singly  ionized  sulfur  vacancy.  This 
identification  establishes  the  singly  ionized  sulfur  vacancy  as  the  dominant  intrinsic 
defect  (i.e.,  donor)  in  SPS  and  has  direct  implications  for  enhancing  the  photorefrac¬ 
tive  properties  of  this  emerging  material. 
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HYPERFINE  INTERACTIONS  IN  THE  ELECTRON 
PARAMAGNETIC  RESONANCE  SPECTRA  OF  POINT 
DEFECTS  IN  WIDE-BAND-GAP  SEMICONDUCTORS 


I.  Introduction 

The  origins  of  magnetic  resonance  date  back  to  the  Nobel-prize-winning  Stern-Gerlach 
experiment  [1].  In  this  experiment,  neutral  silver  atoms  from  a  furnace  are  directed 
through  a  static  magnetic  field  and  captured  on  a  screen  to  measure  their  deflection 
angle.  The  axis  of  the  magnetic  field  is  normal  to  the  initial  path  of  the  silver 
atoms.  If  the  silver  atom  has  a  magnetic  moment,  then  it  will  deflect  according  to  the 
orientation  of  the  magnetic  moment  when  it  enters  the  magnetic  field.  If  the  magnetic 
moment  of  the  silver  atom  is  due  to  a  classical  angular  momentum  and  thus  zero, 
then  the  angular  distribution  of  the  silver  atoms  will  have  only  one  centroid  elongated 
along  the  axis  of  the  magnetic  field.  The  silver  atom  was  observed  to  have  a  magnetic 
moment,  and  the  angular  distribution  was  observed  to  have  two  centroids  displaced 
along  the  axis  of  the  magnetic  field.  This  result  was  explained  by  Uhlenbeck  and 
Goudsmit  [2]  as  a  quantum  interaction  between  the  unpaired  electron  in  the  5s  orbital 
of  the  silver  atom  and  the  magnetic  field.  This  quantum  interaction  only  allows  for 
the  magnetic  moment  of  the  electron  to  align  parallel  or  antiparallel  to  the  magnetic 
held.  Since  an  electron  in  an  s  orbital  has  no  orbital  angular  momentum,  the  magnetic 
moment  is  entirely^  due  to  an  angular  momentum  that  is  intrinsic  to  the  electron. 
Because  of  this,  the  Stern-Gerlach  experiment  established  the  concept  of  intrinsic 

1  Silver  has  two  isotopes  with  I  =  \  that  make  up  100%  of  all  naturally  occurring  silver  atoms 
mm-  Therefore,  the  magnetic  moment  of  a  neutral  silver  atom  is  not  entirely  due  to  the  5s 
electron.  As  seen  in  Section  |2.1[  the  magnetic  moment  of  the  silver  nucleus  is  much  weaker  than 
that  of  the  electron.  The  Stern-Gerlach  experiment  lacked  the  resolution  necessary  to  detect  the 
hyperfine  interaction  between  the  5s  electron  and  the  silver  nucleus.  Breit  and  Rabi  0  described 
this  hyperfine  interaction  in  the  context  of  the  Stern-Gerlach  experiment,  and  Rabi  et  al.  [6]  were 
the  first  to  directly  measure  nuclear  magnetic  moments  with  this  beam  technique  [7]  . 
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quantum  angular  momentum  or  spin  angular  momentum.  A  discussion  of  classical  and 
quantum  angular  momentum  is  included  in  Appendix  [Aj  This  establishment  of  spin 
angular  momentum  gave  rise  to  a  variety  of  magnetic  resonance  techniques,  including 
electron  paramagnetic  resonance,  nuclear  magnetic  resonance,  and  electron-nuclear 
double  resonance.  Both  pulsed  and  continuous-wave  magnetic  resonance  techniques 
have  been  developed.  The  experimental  technique  used  in  the  present  dissertation  is 
cw  electron  paramagnetic  resonance  spectroscopy. 


Electron  paramagnetic  resonance  was  first  observed  by  Zavoisky  in  1944  [7]. 
After  this  discovery,  significant  contributions  by  Abragam,  Bleaney,  Pryce,  and  Van 
Vleck  in  the  early  1950s  led  to  well-developed  theories  of  electron  paramagnetic  reso¬ 
nance  and  well-developed  techniques  of  electron  paramagnetic  resonance  spectroscopy 
[7] .  Electron  paramagnetic  resonance  is  observed  when  defects  or  impurities  with  one 
or  more  electrons  that  do  not  share  a  wave  function  with  another  electron  of  opposite 
spin  are  present.  These  electrons  are  called  unpaired  electrons  or  unpaired  spins. 
As  discussed  in  Chapter  [TT|  these  unpaired  electrons  interact  with  an  applied  mag¬ 
netic  held  creating  magnetic-field-dependent  energy  levels.  Electron  paramagnetic 
resonance  spectroscopy  is  a  collection  of  experimental  techniques  used  for  detecting 
transitions  between  these  energy  levels  in  insulating  and  semiconducting  crystals, 


amorphous  materials,  and  free  radicals.  As  discussed  in  Chapter  III,  the  technique 
used  in  this  dissertation  research  involves  the  observation  of  magnetic-field-dependent 
absorptions  of  microwaves  in  a  resonant  cavity  containing  the  sample.  The  energy  of 
the  microwaves  (i.e. ,  their  frequency)  is  held  constant  while  an  applied  magnetic  held 
splits  the  energy  levels.  An  absorption  occurs  whenever  the  energy  of  the  microwaves 
equals  the  energy  difference  between  two  energy  levels  of  an  allowed  transition.  An 
electron  paramagnetic  resonance  spectrum  in  this  case  is  a  recording  of  the  rehected 
power  from  the  cavity  as  the  magnetic  held  is  swept  across  a  specihed  range.  With 
this  technique,  the  observation  of  hyperhne  interactions  is  particularly  useful  in  the 
identification  of  paramagnetic  defects. 
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As  discussed  in  Section  \2.4{  the  hyperfine  interaction  is  dependent  upon  the 
ground-state  wave  function  of  an  unpaired  electron  or  unpaired  electrons  in  the  pres¬ 
ence  of  nuclei  that  possess  a  magnetic  moment.  Because  of  this,  hyperfine  interactions 
provide  information  about  the  local  physical  environment  of  the  defect,  and  contain 
information  about  the  ground-state  electronic  structure  of  the  defect  in  its  param¬ 
agnetic  state.  In  this  dissertation,  observations  of  hyperfine  interactions  are  used  to 
further  characterize  or  identify  point  defects  in  three  wide-band-gap  semiconductors 
of  particular  military  and  commercial  importance.  These  wide-band-gap  semicon¬ 
ductors  are  zinc  oxide  (ZnO),  titanium  dioxide  (TiC^),  and  stannous  hexathiohy- 
podiphosphate  (S112P2S6).  All  of  these  materials  have  properties  that  are  dependent 
on  the  presence  of  point  defects.  Therefore,  it  is  essential  to  identify  and  characterize 
these  point  defects  in  order  to  understand  the  processes  that  give  rise  to  the  desired 
material  properties. 


1.1  Point  Defects 

A  point  defect  is  an  interruption  in  the  periodicity  of  a  crystal  lattice  that  is 
localized  within  a  few  lattice  sites.  A  crystal  lattice  is  composed  of  a  Bravais  lattice 
which  is  a  periodic  set  of  sites  defined  by  a  set  of  three  vectors.  At  each  Bravais 
lattice  site,  there  is  a  basis  which  is  an  arrangement  of  atoms  about  the  site.  Because 
of  this,  a  crystal  lattice  is  composed  of  sublattices  which  are  Bravais  lattices  formed 
by  each  atom  in  the  basis.  The  simplest  point  defects  are  characterized  as  intrinsic 
or  extrinsic  and  substitutional  or  interstitial.  An  intrinsic  point  defect  involves  atoms 
that  are  native  to  the  crystal,  while  extrinsic  point  defects  involve  atoms  that  are 
not  native  to  the  crystal.  A  substitutional  point  defect  is  located  at  a  crystal  lattice 
site,  while  an  interstitial  point  defect  is  not  located  at  a  crystal  lattice  site.  Intrinsic 
substitutional  point  defects  include  vacancies  and  antisites.  A  vacancy  is  the  absence 
of  an  atom  at  a  crystal  lattice  site,  and  an  antisite  is  a  native  atom  that  is  located 
at  a  sublattice  site  to  which  it  does  not  belong.  Intrinsic  interstitial  point  defects  are 
called  self-interstitials,  and  are  native  atoms  that  are  not  located  at  a  crystal  lattice 


3 


site.  Extrinsic  substitutional  point  defects  are  called  substitutional  impurities,  and 
are  atoms  that  are  not  native  to  the  crystal  and  are  located  at  a  crystal  lattice  site. 
Extrinsic  interstitial  point  defects  are  called  interstitial  impurities,  and  are  atoms 
that  are  not  native  to  the  crystal  and  not  located  at  a  crystal  lattice  site.  These 
simple  point  defects  can  couple  to  form  complex  point  defects.  The  point  at  which  a 
complex  point  defect  becomes  an  extended  defect  can  be  ambiguous.  At  some  point 
as  the  complex  gets  larger,  it  will  take  on  its  own  physical  properties  that  are  separate 
from  that  of  the  crystal,  i.e.,  an  inclusion  is  formed.  As  the  defect  concentration  gets 
even  larger,  all  the  point  defects  will  couple  changing  the  physical  properties  of  the 
entire  crystal,  with  the  result  being  an  alloy. 

A  characteristic  of  any  semiconductor  crystal  is  the  existence  of  a  band  gap 
between  two  broad  electronic  energy  bands.  For  a  flawless  semiconductor  crystal  at 
a  temperature  of  0  K,  one  band  called  the  conduction  band  is  unoccupied  by  elec¬ 
trons  and  the  other  called  the  valence  band  is  completely  occupied  with  electrons.  In 
this  flawless  frozen  state,  the  semiconductor  crystal  is  not  electronically  conductive. 
At  higher  temperatures,  electrons  can  become  thermally  excited  into  the  conduction 
band,  thus  creating  two  current-carrying  particles  and  becoming  electronically  con¬ 
ductive.  One  current-carrying  particle  is  the  electron  in  the  conduction  band  and  the 
other  is  a  pseudo-particle  called  a  hole  in  the  valence  band.  This  happens  at  room 
temperature  for  narrow-band-gap  semiconductors  such  as  silicon  and  germanium. 
For  wide-band-gap  semiconductors,  this  effect  is  minuscule  at  room  temperature. 
However,  point  defects  in  wide-band-gap  semiconductor  crystals  can  create  either 
conduction-band  electrons  or  valence-band  holes  when  ionized  by  thermal  excitation, 
electromagnetic  irradiation,  or  particle  irradiation. 

A  point  defect  that  can  supply  a  conduction-band  electron  when  ionized  is  called 
a  donor.  An  example  of  a  donor  is  a  Group  III  atom  such  as  gallium  substituting 
for  zinc  in  zinc  oxide.  Gallium  which  has  valence  electrons  in  a  4s24p  configuration 
has  one  additional  valence  electron  than  zinc  which  has  valence  electrons  in  a  4s2 
configuration.  When  gallium  is  substituted  into  the  zinc  sublattice  as  a  neutral  donor, 


4 


this  additional  unpaired  electron  occupies  a  state  with  a  highly  diffuse  wave  function 
that  is  centered  on  the  gallium.  This  is  often  referred  to  as  an  effective-mass  or 
hydrogenic  state.  Evidence  for  this  behavior  is  found  in  the  electron  paramagnetic 
resonance  spectrum  of  the  neutral  shallow  donor  in  zinc  oxide.  The  spectrum  consists 
of  a  single  nearly  isotropic  electron  paramagnetic  resonance  line  with  principal  values 
of  the  g  matrix  around  1.96  m  ia  mg.  Since  both  naturally  occurring  isotopes  of 
gallium  have  /  =  §  mm,  the  unresolved  and  weak  hyperhne  interaction  suggests 
that  the  electron  is  not  localized  at  the  gallium,  but  is  spread  out  (i.e. ,  shared)  over 
many  shells  of  ions  around  the  gallium  impurity. 

A  point  defect  that  can  capture  a  valence-band  electron  creating  a  hole  is  called 
an  acceptor.  An  example  of  an  acceptor  is  a  Group  V  atom  such  as  nitrogen  sub¬ 
stituting  for  oxygen  in  zinc  oxide.  Nitrogen  which  has  valence  electrons  in  a  2s22p 3 
configuration  has  one  less  valence  electron  than  oxygen  which  has  valence  electrons 
in  a  2s22p 4  configuration.  When  nitrogen  is  substituted  into  the  oxygen  sublattice 
as  a  neutral  acceptor,  this  lack  of  an  electron  leaves  an  unpaired  electron  in  a  state 
that  is  highly  localized  about  the  nitrogen.  Evidence  for  this  is  given  in  Chapter  |IV| 
In  general,  the  wave  functions  of  shallow  acceptors  are  more  localized  than  the  wave 
functions  of  shallow  donors. 

1.2  Zinc  Oxide 

Zinc  oxide  is  a  versatile  material  that  has  been  used  for  centuries  in  paints, 
pharmaceuticals,  and  ceramics  and  since  the  middle  of  the  previous  century,  it  has 
been  investigated  for  use  as  a  semiconductor.  Its  refractive  properties  and  fine  particle 
size  when  produced  by  the  “American”  or  “French”  methods  give  paints  and  phar¬ 
maceuticals  such  as  sunscreen  and  cosmetics  high  “covering  power”  and  give  ceramic 
glazes  “brilliance”  m  By  the  1930s,  some  physical  and  optical  properties  of  zinc 
oxide  were  being  explored  [Qj.  By  the  1960s  it  was  being  investigated  for  its  semi¬ 
conductor  properties  ra-  Properties  of  interest  include  a  wide  direct  band  gap  and 
the  possibility  of  band-gap  engineering  by  alloying  with  cadmium  and  magnesium. 
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The  band  gap  of  zinc  oxide  is  about  3.3  eV  [2]  which  classifies  it  as  a  wide- 
band-gap  semiconductor.  This  property  is  useful  in  high-voltage  electronics,  since 
wide-band-gap  semiconductors  are  less  susceptible  to  breakdown  where  high  voltages 
can  cause  charge  carriers  to  create  electron-hole  pairs  through  impact  ionization  or 
to  tunnel  through  p-n  junctions.  A  wide  band  gap  also  reduces  the  requirement 
that  circuits  be  insulated  from  exposure  to  light  because  visible  photons  lack  enough 
energy  to  create  electron-hole  pairs.  Since  silicon-based  devices  have  a  band  gap 
of  about  1.1  eV  [15],  visible  and  infrared  photons  can  produce  excess  electron-hole 
pairs  in  the  device  which  can  significantly  affect  device  performance.  Due  to  its 
wide  band  gap,  zinc  oxide  is  transmissive  to  visible  and  infrared  photons  making  it 
useful  as  a  transparent  conducting  window  for  use  in  solar  cells  and  displays,  and 
absorptive  to  ultraviolet  photons  making  it  useful  as  an  ultraviolet  photodiode  for 
use  in  ultraviolet  detectors.  Another  aspect  of  the  band  gap  in  zinc  oxide  is  that  it  is 
direct  where  the  conduction-band  minimum  and  the  valence-band  maximum  appear 
at  the  same  location  in  momentum  space.  This  makes  luminescent  transitions  more 
probable  and  with  proper  band-gap  engineering,  zinc  oxide  is  suitable  for  ultraviolet 
light  emitting  devices  such  as  light  emitting  diodes  and  laser  diodes. 

Since  zinc  oxide  is  a  II- VI  semiconductor,  band-gap  engineering  is  possible  by 
alloying  with  Group  IIA  and  IIB  metals,  and  Group  VI  nonmetals  and  metalloids. 
The  most  promising  elements  for  band-gap  engineering  are  magnesium  and  cadmium 
[HI  EH  US].  Magnesium  oxide  has  a  band  gap  of  about  8.7  eV  na  ,  so  the  incorporation 
of  magnesium  increases  the  band  gap  of  zinc  oxide  [18].  Cadmium  oxide  has  a  band 
gap  of  about  2.3  eV  [19],  so  the  incorporation  of  cadmium  decreases  the  band  gap 
of  zinc  oxide  [T5].  This  ability  to  engineer  the  band  gap  is  essential  to  the  successful 
production  of  zinc  oxide  based  heterostructures  which  are  the  building  blocks  of  light 
emitting  diodes,  laser  diodes,  transistors,  and  many  other  electronic  and  photonic 
devices. 
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1.3  Titanium  Dioxide 


Like  zinc  oxide,  titanium  dioxide  has  found  many  uses  in  low-tech  and  high- 
tech  applications.  Titanium  dioxide  has  been  used  as  a  white  pigment  in  paints  since 
the  early  part  of  the  20th  century  because  of  its  high  refractive  index  [20] .  It  is  also 
nontoxic  making  it  a  replacement  for  lead-based  “white  lead”  paints  [20].  It  has  also 
found  a  use  as  a  diamond  simulant  in  jewelry  dJ  Recently,  research  on  titanium 
dioxide  has  focused  on  its  photocatalytic  [ZB  122]  and  ferromagnetic  [231, 123]  properties, 
and  use  as  a  memristor  [25]  and  battery  terminal  [26]. 

Of  particular  use  to  the  military  are  the  photocatalytic  properties  of  titanium 
dioxide  and  use  as  a  battery  terminal  in  lithium  ion  batteries.  Joint  terminal  air 
controllers,  combat  controllers,  and  combat  weathermen  rely  on  compact  and  light¬ 
weight  power  systems  to  operate  essential  equipment  such  as  radios,  designators, 
sensors,  and  beacons.  These  operators  routinely  find  themselves  in  remote  and  diffi¬ 
cult  terrain  where  they  must  carry  their  equipment  and  the  power  sources  to  operate 
them.  Therefore,  any  improvement  in  the  capacity  or  efficiency  of  these  power  sources 
will  increase  the  operators  ability  to  complete  their  mission.  As  a  photocatalyst,  ti¬ 
tanium  dioxide  has  the  ability  to  hydrolyze  water  in  the  presence  of  ultraviolet  light 
EH-  Hydrogen  gas  produced  by  this  method  can  provide  power  through  a  fuel  cell, 
thus  increasing  capacity  as  long  as  there  is  water  and  sunlight  to  fuel  the  process. 
Since  titanium  dioxide  has  better  lithium-transport  properties  than  current  carbon- 
based  materials,  efficiency  gains  can  be  made  by  replacing  the  carbon  electrodes  in 
lithium  ion  batteries  with  titanium  dioxide  electrodes  [26]. 

1.4  Stannous  Hexathiohypodiphosphate 

Recent  studies  of  stannous  hexathiohypodiphosphate  (SPS)  have  focused  on  its 
ferroelectric  S3  and  photorefractive  [251  [2D,,  0Q]  properties.  Although  these  proper¬ 
ties  have  been  well  established  in  SPS,  the  point  defects  participating  in  these  effects 
are  not  well  established.  Therefore,  experimental  techniques  such  as  electron  param¬ 
agnetic  resonance  spectroscopy  are  needed  to  identify  the  point  defects  and  correlate 
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their  presence  and  properties  to  ferroelectric  and  photorefractive  effects  in  the  mate¬ 
rial.  This  information  is  subsequently  used  to  monitor  these  defects  for  production 
of  material  with  the  desired  properties. 

Of  particular  use  to  the  military  are  the  high-speed  photorefractive  properties 
of  SPS.  High-sensitivity  visible  and  infrared  sensor  systems  have  become  an  essential 
component  in  most  weapons  platforms.  These  sensors  are  used  for  missions  such  as 
surveillance,  reconnaissance,  targeting,  tracking,  and  guidance.  Because  these  sensor 
systems  are  composed  of  highly  photosensitive  components,  they  are  prone  to  laser- 
induced  damage  and  jamming.  The  proliferation  of  highly  portable  high-power  visible 
and  infrared  laser  systems  has  led  to  the  need  to  protect  these  sensor  systems.  One 
method  of  protection  is  to  direct  the  damaging/jamming  laser  radiation  away  from 
the  photosensitive  components.  This  is  accomplished  with  photorefractive  materials, 
because  the  refractive  index  of  these  materials  is  dependent  upon  the  intensity  of 
light  illuminating  them.  An  intensity  pattern  created  by  the  interference  between 
the  damaging/jamming  laser  and  a  local  laser  “writes”  a  diffraction  grating  into  the 
photorefractive  material.  This  diffraction  grating  then  deflects  the  damaging/jam¬ 
ming  laser  radiation  away  from  the  photosensitive  components  of  the  sensor  system, 
thus  protecting  the  photosensitive  components  from  damage,  and/or  eliminating  the 
jamming  signal. 


II.  Principles  of  Electron  Paramagnetic  Resonance 

Spectroscopy 

In  general,  spectroscopy  is  the  detection  of  discrete  transitions  between  energy  levels 
present  within  a  material.  These  transitions,  or  resonances,  are  detected  by  expos¬ 
ing  the  material  to  a  source  of  photons  (e.g.,  ultraviolet  light,  visible  light,  infrared 
light,  or  microwaves)  and  measuring  either  a  reflection  or  absorption  of  these  photons 
in  the  material.  Electron  paramagnetic  resonance  spectroscopy  involves  transitions 
between  energy  levels  created  by  the  interactions  of  one  or  more  unpaired  electrons 
with  a  magnetic  held  and  any  other  magnetic  moments  present  in  the  material.  These 
interactions  are  described  by  a  spin  Hamiltonian  and  each  interaction  adds  its  own 
component  to  the  spin  Hamiltonian.  With  a  knowledge  of  the  specific  interactions 
present  in  a  particular  spin  system,  a  total  spin  Hamiltonian  is  constructed  as  a  sum 
of  each  interaction’s  contribution  to  the  spin-dependent  energy. 

2. 1  Magnetic  Moments 

For  a  classical  electrodynamic  system  [31] ,  the  magnetic  held  created  by  a  gen¬ 
eral  current  distribution  is 


B(r)  =  x  [  J^\d3r'  (1) 

4vr  y  |r  —  r'| 

where  /io  is  the  permeability  of  free  space,  J  is  the  current  density,  r  is  the  position 
vector  of  the  magnetic  held,  r'  is  the  position  vector  of  the  current  density,  and  d3 r' 
indicates  integration  over  the  volume  containing  the  current  density.  For  a  system  of 
localized  current  density  placed  about  the  origin,  the  denominator  of  Equation  (|TJ)  is 
expanded  in  a  Maclauren  series  with  respect  to  the  primed  coordinates  to  approximate 
the  magnetic  held  when  |r|  >>  |r7 1 , 

1  1  r  •  r' 

I - 7?  =  pi  +  -ryy  +  •  •  •  (2) 
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Substituting  this  expansion  into  Equation  ((l])  results  in 

B(r)  =  £vx  J(r,)d3r,+ J(r')(r  ■  r')dV  +  . . .  V  (3) 

The  first  term  in  this  equation  is  called  the  monopole  term  and  the  second  is  called  the 
dipole  term.  Since  the  current  density  is  assumed  to  be  localized  with  no  divergence, 
the  monopole  term  is  zero  and  the  dipole  term  can  be  written  as 

j^s  /  J(r')(r-r')dV=  (4) 

where 


m  =  -  j  r'  x  J(r/)d3r/ 


(5) 


is  called  the  magnetic  moment.  With  these  transformations,  the  magnetic  field  outside 
of  a  localized  current  density  is  approximated  as 


B(r)  = 


47t  |r|d 
fM) 


47T  |  r  | 


2  rn  —  3- 


r  x  m  x  r 


(6) 


471 1  r 


p0  /  Q  fr  •  mjr 

3 - r-rr; - III 


An  example  of  a  simple  magnetic  moment  is  that  of  a  charged  particle  in  a 
circular  orbit.  The  current  generated  by  the  orbiting  particle  is 


I 


cjv 

27r|r| 


(7) 


and  the  orbital  angular  momentum  is 


L  =  r  x  mv 


(8) 
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where  q  is  the  charge  of  the  particle,  v  is  the  velocity  of  the  particle,  r  is  the  orbital 
position  of  the  particle,  and  m  is  the  mass  of  the  particle.  Since 


J(r')d3r'  — >  I(r')dr' 


(9) 


for  hlamental  currents  in  a  closed  circuit,  substituting  Equations  ([T])  and  ([8])  into 
Equation  (J5j)  results  in 

1  q 

m  =  -  (r  x  I)  27T | r  |  =  L.  (10) 

This  is  the  classical  relation  between  magnetic  moment  and  orbital  angular  momen¬ 
tum. 

When  this  concept  is  applied  to  quantum  mechanical  systems,  the  magnetic 
moment  is  viewed  as  a  magnetic  moment  operator, 


m  = 


— L 

2  m 


(11) 


The  magnetic  moment  due  to  spin  angular  momentum  can  differ  from  that  of  the 
classical  result.  This  is  due  to  other  quantum  mechanical  and  relativistic  effects  not 
discussed  here.  These  effects  are  accounted  for  with  a  g  factor, 


m  =  q - S. 

y  2  m 


(12) 


For  a  free  electron,  the  magnetic  moment  operator  is 


m  =  —ge 


W S 

2  me 


-ge^B  S 


(13) 


where  ge  is  the  free  electron  g  factor  and  is  known  to  be  2.0023  [32],  e  is  the  elementary 
charge,  me  is  the  electron  mass,  and  /ib  is  the  Bohr  magneton  and  is  known  to  be 
9.2740  x  10~24  J/T  [32].  For  a  nucleus,  the  magnetic  moment  operator  is 


Zmn 


(14) 
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where  gn  is  the  g  factor  of  a  particular  nucleus  and  is  on  the  order  of  unity,  mp  is  the 
proton  mass,  and  //jv  is  the  nuclear  magneton  and  is  known  to  be  5.0508  x  10“ 2 '  J/T 
[32 j .  Since  the  order  of  the  nuclear  g  factor  is  unity  and  the  nuclear  magneton  is 
three  orders  of  magnitude  less  than  the  Bohr  magneton,  the  electron  has  a  magnetic 
moment  that  is  significantly  greater  than  all  known  nuclei.  The  nuclear  g  factor 
typically  carries  the  sign  of  the  magnetic  moment  because  neutral  particles,  like  the 
neutron,  have  a  magnetic  moment  even  though  they  carry  no  charge  and  some  nuclei 
have  a  magnetic  moment  that  is  in  opposition  to  its  spin  even  though  they  carry  only 
positive  charge. 

2. 2  Spin  Hamiltonian 

A  quantum  mechanical  system  can  be  in  any  number  of  quantum  mechanical 
states  represented  by  the  kets  \ijjn)  with  the  requirement  that  the  kets  obey  the  time 
independent  Schrodinger  equation. 

H  | i>n)  =  En  | i>n)  (15) 

where  H  is  the  Hamiltonian  operator  and  En  is  the  energy  of  the  state.  The  Hamilto¬ 
nian  includes  terms  that  account  for  the  kinetic  energy  of  the  system  and  the  potential 
energy  of  the  system. 

H  =  T  +  U.  (16) 

In  this  equation,  T  accounts  for  kinetic  energy  and  U  accounts  for  potential  energy. 
The  potential  energy  can  be  further  separated  into  parts  that  account  for  potentials 
that  are  purely  a  function  of  spatial  coordinates  and  potentials  that  are  dependent 
on  the  spin  angular  momenta  of  the  particles  that  make  up  the  system. 

H  =  f  +  Ur  +  Us  =  Hr  +  Ua  (17) 
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where  Ur  is  the  spatial  potential  and  only  contains  purely  spatial  contributions  to 
the  potential  energy,  Us  is  the  spin  potential  and  only  contains  potentials  that  are 
dependent  on  spin,  and  Hr  is  the  spatial  Hamiltonian  and  only  contains  purely  spatial 
contributions  to  the  total  energy. 

The  spin  potential  can  be  treated  as  a  perturbation  to  the  ground  state  of  the 
spatial  Hamiltonian.  The  first-order  perturbation  to  the  ground-state  energy  is 

S0  =  4r)  +  (4r)|^|4r)}  (18) 

where  represents  the  ground  state  of  the  spatial  Hamiltonian  and  is  the 

ground-state  energy  of  the  spatial  Hamiltonian  such  that 

itr  |dr))  =  -Ef  Id'1)  •  (19) 

Since  the  spin  of  a  particle  is  an  intrinsic  characteristic,  it  is  reasonable  to  assume 
that  the  spin  operators  commute  with  the  spatial  Hamiltonian,  therefore  sharing  a 
common  set  of  eigenstates.  Because  of  this,  the  ground  state  can  be  represented  by 
the  kets  rris% ,  ^ ,  and  the  spatial  part  of  the  ground  state  can  be  treated  separately 

from  the  spin  part  of  the  ground  state. 

To  illustrate  this  concept,  suppose  that  a  spin-dependent  potential  has  the  form 

Us  =  f(r)Sz.  (20) 

The  first-order  perturbation  of  this  operator  to  the  ground-state  energy  is 

( ms ,ipo]  Sgf{  r)  ms,ipo)S}  =  ms  J  (r)/(r)^r)  (r)d3r.  (21) 


13 


Since  the  spatial  part  of  the  ground  state  can  be  considered  separately  from  the  spin 
part,  the  spin  Hamiltonian  for  this  operator  is  defined  as 

Hs  =  Sz  J  '0o')* (r ) / (r ) V'cf) (r ) ^3r  (22) 

which  is  the  original  operator  with  the  spatial  parts  averaged.  Therefore,  the  spin 
Hamiltonian  is  defined  as  the  spin-dependent  potentials  of  the  Hamiltonian  with  the 
spatial  parts  averaged  according  to  the  ground  state  of  the  spatial  Hamiltonian. 

The  spin  Hamiltonian  represents  perturbations  to  the  ground  state  of  the  sys¬ 
tem.  In  most  cases,  differences  in  energy  (resonances)  are  measured  in  an  experiment. 
Because  of  this,  the  ground-state  energy  can  be  ignored  and  the  spin  Hamiltonian 
can  be  used  to  explain  transitions  between  energy  levels  in  a  spin  system.  These  phe¬ 
nomena  include  but  are  not  limited  to  the  Zeeman  effect,  fine  interactions,  hyperhne 
interactions,  and  nuclear  electric  quadrupole  interactions.  All  of  these  phenomena  are 
appropriate  for  describing  energy  transitions  in  a  spin  system  with  a  single  unpaired 
electron  (i.e.,  an  S  =  \  system),  and  are  described  in  the  following  sections. 

2. 3  Zeeman  Effect 

The  Zeeman  effect  is  caused  by  a  change  in  the  potential  energy  of  a  magnetic 
moment  due  to  the  presence  of  a  magnetic  held.  It  is  a  quantum  mechanical  effect,  and 
gives  rise  to  electron  paramagnetic  resonance  phenomena  when  applied  to  electrons 
and  to  nuclear  magnetic  resonance  phenomena  when  applied  to  nuclei.  For  a  classical 
electrodynamic  system  [31j ,  the  potential  energy  of  a  magnetic  moment  in  a  magnetic 
held  is 

U  =  — m  •  B.  (23) 

The  magnetic  moment  in  this  equation  is  replaced  with  the  magnetic  moment  operator 
of  an  electron  or  nucleus  to  arrive  at  the  spin  Hamiltonian. 
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The  spin  Hamiltonian  of  a  free  electron  is 


Hs  =  •  B.  (24) 

Since  the  electron  is  a  half-spin  particle,  the  eigenstates  of  this  spin  Hamiltonian  are 
represented  by  the  kets  | mg)  where  mg  G  {+|,  —  Selection  of  the  projection  axis 
of  the  spin  operator  is  arbitrary,  therefore  it  is  common  to  place  the  projection  axis 
along  the  magnetic  field  to  simplify  the  equation.  Doing  this  results  in 

Hs  =  geiiBBSz  (25) 


where  B  is  the  magnitude  of  the  magnetic  field.  The  action  of  the  spin  Hamiltonian 
on  the  eigenstates  results  in 


H„ 


H„ 


de^B  B  S ' z 
ge^BBSz 


9^bB 

f  2 

_9e^BB 

2 


and 


(26) 


where  E+i  and  E_i  are  the  two  energy  levels  of  the  system.  Electron  paramagnetic 
resonance  is  observed  by  detecting  the  absorption  of  photons  due  to  a  transition  from 
the  lower  energy  |  —  |)  state  to  the  higher  energy  |  +  ^)  state.  The  energy  of  this 
transition  is 

E  =  E+ 1  —  E_  i  =  geHBB  =  hi y  (27) 

where  vr  is  the  frequency  of  the  photon  necessary  for  absorption.  The  energy  levels 
as  a  function  of  magnetic  field  are  presented  in  Figure  [T|  for  a  free  electron. 

The  spin  Hamiltonian  of  a  nucleus  is 


Hs  fiti.P'ivi  ■  B. 


(28) 
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Figure  1.  Energy  levels  due  to  the  Zeeman  effect  for  a  free  electron  as  a  func¬ 
tion  of  magnetic  field.  The  vertical  double  arrow  line  corresponds  to  an  electron 
paramagnetic  resonance  of  10  GHz. 

Selection  of  the  projection  axis  along  the  magnetic  field  results  in 


Hs  = 


dnl^N  B I z 


(29) 


The  eigenstates  of  this  spin  Hamiltonian  are  represented  by  kets  | mj)  where  rrij  e 
{— /,  —I  +  1, . . . ,  /  —  1, 1}  and  /  G  { 0,  1,  §, . . .  }.  Since  a  nucleus  can  have  a  spin 

quantum  number  greater  than  one,  there  are  2/  +  1  energy  levels  with  values  of 


h'//(  i  Snl^N  B  TTlj . 


(30) 


Nuclear  magnetic  resonance  is  observed  by  detecting  the  absorption  of  photons  due 
to  a  transition  from  a  lower  energy  state  to  a  higher  energy  state.  The  selection  rule 
for  this  process  is  that  the  absolute  value  of  the  difference  between  the  spin  projection 
quantum  number  of  the  final  state  and  the  spin  projection  quantum  number  of  the 
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initial  state  must  equal  one.  Because  of  the  selection  rule,  the  energy  of  this  transition 

is 

E  —  \EmiJr\  Emi  |  \dn\ UnE  —  hvr.  (31) 

In  general,  the  electron  is  not  free  and  can  be  bound  within  an  atom,  molecule, 
or  point  defect.  In  this  case,  the  electron  can  have  orbital  angular  momentum  in 
addition  to  spin  angular  momentum.  The  spin  Hamiltonian  becomes 

Hs  =  iiB(geS  +  L)  •  B  +  AS  •  L  =  gefiB S  •  B  +  ^/i^B  +  Asj  •  L  (32) 

where  A  is  a  spin-orbit  coupling  coefficient.  It  is  necessary  to  write  this  spin  Hamil¬ 
tonian  in  terms  of  the  spin  operator  alone.  This  is  accomplished  by  applying  a 
second-order  perturbation  [7j.  Doing  this  results  in 

Hs  =  ge/iB S  •  B  +  (/iB)2  B  A  B  +  2A/iBS  ■  A  •  B  +  A2S  ■  A  ■  S.  (33) 


The  elements  of  the  A  matrix  are 


A«  =  ~EEE 


( r ) 


^n]  )  (  $ 


'>) 


4] 


i(r) 


n=l 


E, 


(r) 


(34) 


where  i,j  G  {x,y,z},  is  the  ground  state  of  the  spatial  Hamiltonian,  is  the 

(r) 

ground-state  energy  of  the  spatial  Hamiltonian,  'ijjn  '  is  an  excited  state  of  the  spatial 
Hamiltonian,  and  En  ^  is  an  excited-state  energy  of  the  spatial  Hamiltonian.  The 


second  term  in  Equation  (33)  adds  the  same  amount  of  energy  to  all  energy  levels. 


Because  of  this,  it  can  be  ignored  when  observing  energy  transitions.  The  last  term  in 


Equation  ( 33 )  adds  the  same  amount  of  energy  to  all  energy  levels  whenever  the  spin 


quantum  number  is  equal  to  Because  of  this,  it  can  be  ignored  when  observing 
energy  transitions  in  a  half-spin  system.  Taking  this  into  consideration,  the  spin 
Hamiltonian  becomes 

Hs  =  yB S  ■  g  ■  B  +  S  ■  D  ■  S  (35) 
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where 

g  =  ge  1  +  2AA,  and  D  =  A2A.  (36) 

Matrix  g  is  called  a  g  matrix,  matrix  D  is  called  a  fine  interaction  matrix,  and  matrix 
1  is  the  identity  matrix. 


2.4  Hyperfine  Interaction 

The  hyperfine  interaction  is  an  interaction  between  magnetic  moments  and  the 
magnetic  fields  that  they  generate.  It  is  a  quantum  mechanical  effect,  and  is  observed 
in  electron  paramagnetic  resonance  spectroscopy  when  the  magnetic  moment  of  an 
electron  interacts  with  that  of  a  nucleus.  When  the  magnetic  field  is  generated  by 
another  magnetic  moment  placed  at  some  origin,  Equation  (]6j)  is  substituted  into 


Equation  (23)  to  arrive  at 


U  = 


ho 


47T I  r  | 


in,  ■  rn0  —  3 


1  •  mo)(m,  ■  r) 


(37) 


where  r  is  the  position  vector  of  magnetic  moment  mr,  and  m0  is  the  magnetic 
moment  at  the  origin.  Since  the  hyperfine  interation  is  between  the  magnetic  moments 
of  an  electron  and  a  nucleus,  the  magnetic  moment  at  the  origin  is  replaced  with  the 
magnetic  moment  operator  of  a  nucleus,  and  the  magnetic  moment  at  r  is  replaced 
with  the  magnetic  moment  operator  of  an  electron  to  arrive  at 


U  = 


hO  hehshnhW 
47T 


(r  •  I)(S  •  r)  SI 


(38) 


It  is  convenient  to  place  the  nucleus  at  the  origin  and  the  electron  at  r,  since  in 
most  cases  the  nucleus  is  relatively  stationary  while  the  electron  can  be  at  many 
other  positions  with  a  probability  defined  by  its  ground-state  spatial  wave  function. 


Expanding  Equation  ( 38 )  results  in 


jj  _  l^09e^Bgn^N  ST''  t  ^ j 


47T 


*  j 


h  \  sj 


1  J 


(39) 


18 


where  i,j  G  {x,y,z},  and  Sij  is  a  Kronecker  delta.  In  order  to  arrive  at  the  spin 
Hamiltonian,  the  parenthetical  expression  in  Equation  (39)  must  be  replaced  by 


trS-rrr)  /  4r)*W  ( 3r‘rj 


r|5  Wor)(rM3r 


(40) 


where  i/jq'1  is  the  ground-state  spatial  wave  function  of  the  electron.  After  accom¬ 
plishing  this,  the  spin  Hamiltonian  becomes 


if,  =  S  ■  A  •  I 


(41) 


where  A  is  called  the  hyperhne  interaction  matrix  and  its  elements  are 


A,,  =  W^“■l‘l,  /  4r>*(r)  [  3 


nr>  ~  rk )  df’M dh. 


(42) 


Equation  (40)  becomes  problematic  whenever  r  is  a  null  vector.  This  occurs  whenever 


the  electron  occupies  an  s  orbital  about  a  nucleus.  Because  of  this,  the  hyperhne 
interaction  requires  the  addition  of  a  contact  interaction  to  account  for  an  electron 
in  an  s  orbital.  The  potential  energy  of  this  interaction  is 


U  =  --/io|V,or)(0)|2me-m,1 


(43) 


where  |V’o^(0)|2  is  the  probability  that  the  electron  is  at  the  nucleus,  me  is  the  mag¬ 
netic  moment  of  the  electron,  and  m„  is  the  magnetic  moment  of  the  nucleus  [33]  - 
Replacing  the  magnetic  moments  with  their  appropriate  quantum  expressions  results 
in  the  spin  Hamiltonian 


Hs  =  ^po5,eAiB5,nhAr|'0of)(O)|2S  i  =  AS  1 


(44) 


The  coefficient  A  is  called  the  isotropic  hyperhne  interaction  parameter  or  the  Fermi 
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Figure  2.  Energy  levels  due  to  magnetic  interactions  between  an  electron  and  a 
proton  in  a  hydrogen  atom  as  a  function  of  magnetic  field.  The  vertical  double  arrow 
lines  correspond  to  electron  paramagnetic  resonances  of  10  GHz.  The  two  high-energy 
levels  are  more  closely  spaced  than  the  two  low-energy  levels  due  in  part  to  the  nuclear 
Zeeman  effect.  The  two  high-energy  levels  will  eventually  cross  at  about  16.655  T 
where  the  nuclear  Zeeman  effect  becomes  equivalent  to  the  hyperhne  interaction. 


contact  term.  Addition  of  Equation  (44)  to  the  spin  Hamiltonian  results  in 


Hs  =  AS  •  I  +  S  •  A  •  i  =  S  ■  [A1  +  A]  •  i  — ►  S  ■  A  •  I  (45) 


where  A1  is  the  isotropic  hyperhne  interaction  matrix. 

An  example  of  a  spin  system  with  a  hyperhne  interaction  is  the  hydrogen  atom. 
The  spin  Hamiltoniam  for  this  spin  system  is 


Hs  =  ge^EiS  B  +  AS  i  -  gpfiNl  ■  B 


(46) 


where  gp  is  the  free  proton  g  factor  and  is  known  to  be  5.5857  [32]  and  A  is  the 
isotropic  hyperhne  interaction  parameter  of  hydrogen  and  is  known  to  be  1420.4  MHz 
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[34].  The  eigenstates  of  the  spin  operators  are  represented  by  kets  | ms,  mi)  where 
7ns,iTij  G  {  — |,  +  4}.  The  spin  Hamiltonian  is  represented  in  this  spin  basis  as  a  4  x  4 
matrix.  This  spin  Hamiltonian  matrix  is  diagonalized  to  extract  the  four  energy  levels 
and  four  eigenstates.  A  plot  of  the  energy  levels  as  a  function  of  magnetic  field  is 
presented  in  Figure  [2] 


2.5  Nuclear  Electric  Quadrupole  Interaction 

Since  all  nuclei  except  for  1 H  are  composed  of  multiple  protons  and  neutrons, 
these  nuclei  will  have  a  definitive  charge  density.  If  this  charge  density  is  not  spher¬ 
ically  symmetric,  then  the  electric  potential  is  dependent  on  the  orientation  of  the 
nucleus.  This  orientation  is  defined  by  the  spin  angular  momentum  of  the  nucleus. 
This  gives  rise  to  a  spin  dependent  potential  that  must  be  included  in  the  spin  Hamil¬ 
tonian.  The  electric  potential  in  terms  of  the  spherical  harmonics  is 


1  Yim( r) 

7°h,h  p+wf 


yL(r')\v'fp(r')d 


'V 


(47) 


where  Y(m  is  a  spherical  harmonic  and  p  is  the  charge  density  [31] .  Due  to  symmetry 
rules  in  the  standard  model  of  particle  physics,  only  potentials  with  even  values  of 
l  that  satisfy  l  <  2 1  are  nonvanishing  [35J.  The  l  =  0  term  is  called  the  monopole 
term  and  gives  rise  to  a  spherically  symmetric  potential  that  is  independent  of  the 
spin  angular  momentum.  Because  of  this,  the  monopole  term  is  not  included  in 
the  spin  Hamiltonian.  Therefore,  the  most  significant  terms  in  the  electric  potential 
that  contribute  to  the  spin  Hamiltonian  are  the  l  =  2  terms.  These  terms  are  called 
quadrupole  terms.  The  quadrupole  potential  can  be  put  into  terms  of  the  spin  angular 
momentum  operators  by  evoking  the  Wigner-Eckart  theorem  [33],  This  results  in 


Hs  =  i  Q  i 


(48) 
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where  in  the  principal  axis  basis  of  the  matrix 


Q  — y 


e2qQ 

2/(2/  -  1) 


-1 

0 

0 

1 

0 

0 

0 

-1 

0 

+  F 

0 

-1 

0 

0 

0 

2 

0 

0 

0 

(49) 


In  this  equation,  F  is  an  anisotropy  factor,  eq  is  the  electric-held  gradient  at  the 
nucleus,  and 

e< Q  =  2\[y  J  ^20 ( r 7 )  I r 7 1 2 P ( r 7 ) r 7  (50) 

where  Q  is  called  the  nuclear  electric  quadrupolc  moment  and  is  usually  given  in 
barns  (1CT28  m2). 
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III.  Instrumentation  of  Electron  Paramagnetic  Resonance 

Spectrometers 

An  electron  paramagnetic  resonance  spectrometer  is  designed  to  detect  energy  tran¬ 
sitions  in  a  spin  system.  It  consists  of  a  source  of  microwaves,  a  microwave  detector, 
and  a  large  electromagnet.  Typical  water-cooled  electromagnets  can  create  magnetic 
fields  of  up  to  1.5  T.  For  a  free  electron,  this  electromagnet  can  produce  splittings 
in  the  energy  levels  corresponding  to  frequencies  of  up  to  42  GHz.  Electromagnetic 
radiation  with  frequencies  of  1  —  100  GHz  covers  the  L,  S,  X,  Q,  and  W  microwave 
bands  of  the  electromagnetic  spectrum  [36].  Transitions  within  a  spin  system  caused 
by  microwaves  can  be  measured  by  absorption  where  a  material  is  placed  within  a 
waveguide  that  connects  the  microwave  source  and  detector,  or  by  reflection  where  a 
material  is  placed  within  a  resonant  cavity  that  reflects  microwave  radiation  when¬ 
ever  the  resonant  frequency  or  the  Q  of  the  cavity  is  changed  by  absorption  in  the 
material. 

The  instrument  used  in  the  present  dissertation  research  is  a  Bruker  EMX  spec¬ 
trometer  with  an  Oxford  Instruments  helium  gas  flow  cryostat  for  temperature  control 
down  to  4.2  K.  This  electron  paramagnetic  resonance  spectrometer  operates  in  the 
microwave  X  band  (8  —  12  GHz),  and  detects  transitions  via  reflections  from  a  reso¬ 
nant  microwave  cavity.  These  reflections  are  detected  and  recorded  while  sweeping  the 
magnitude  of  an  applied  magnetic  field  across  a  specified  range  to  generate  electron 
paramagnetic  resonance  spectra. 

3.1  Microwave  Bridge  and  Cavity 

Microwaves  are  produced,  detected,  and  controlled  by  a  collection  of  compo¬ 
nents  referred  to  as  the  microwave  bridge.  A  schematic  of  a  microwave  bridge  is 
shown  in  Figure  [3j  The  source  is  a  solid-state  Gunn  diode  which  excites  resonant 
microwave  frequencies  within  a  tunable  source  cavity.  The  microwaves  are  coupled 
into  a  wave-guide  which  guides  the  microwaves  to  a  splitter.  At  this  point,  part  of 


23 


the  microwave  radiation  is  guided  along  a  reference  arm,  and  part  is  guided  to  the 
microwave  cavity.  The  reference  arm  is  composed  of  an  attenuator  and  a  phase  shifter 
and  serves  as  a  means  to  bias  the  detector  diode.  Biasing  the  detector  diode  is  neces¬ 
sary  to  ensure  that  it  operates  in  a  linear  mode  (i.e. ,  a  relative  change  in  microwave 
power  leads  to  the  same  relative  change  in  signal).  The  attenuator  in  the  reference 
arm  is  adjusted  to  create  the  desired  amount  of  bias  in  the  detector  diode  and  the 
phase  shifter  is  adjusted  to  match  the  phase  of  the  reference  arm  microwaves  to  the 
phase  of  the  microwaves  reflected  from  the  cavity. 


Reference  Arm  Reference  Arm 
Attenuator  Phase  Shifter 


Figure  3.  Schematic  of  a  microwave  bridge. 


Microwaves  that  are  directed  to  the  cavity  first  go  through  an  attenuator  which 
varies  the  power  of  the  microwaves  according  to  the  needs  of  the  experiment,  then  to  a 
circulator  which  guides  the  microwaves  to  the  cavity.  To  achieve  maximum  sensitivity, 
the  microwaves  are  critically  coupled  into  the  cavity.  This  means  that  almost  all  of 
the  microwave  energy  enters  the  cavity  independent  of  the  incident  amount  of  energy. 
This  critical  coupling  is  achieved  by  inserting  or  withdrawing  a  teflon  screw  with  a 
small  piece  of  metal  on  the  tip  in  front  of  the  cavity’s  iris  (Figure  [d]).  Once  inside  the 
cavity,  the  microwaves  form  a  standing  wave  pattern.  Figure  [5]  shows  the  magnetic 
and  electric  fields  of  the  standing  wave  pattern  in  a  circular  cavity.  Since  magnetic 
moments  in  the  sample  interact  with  the  magnetic  held  of  the  standing  wave,  the 
cavity  is  designed  to  place  the  magnetic  held  maximum  of  the  microwaves  at  the 
location  of  the  sample.  In  Figure  |5j  the  sample  is  inserted  into  the  cavity  through  an 
access  hole  in  the  top.  The  cavity  is  positioned  within  the  gap  of  the  electromagnet 
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such  that  the  magnetic  field  of  the  microwaves  is  normal  to  the  applied  magnetic  field 
generated  by  the  electromagnet. 


Microwave  Cavity 

Figure  4.  Schematic  of  microwave  energy  coupling  into  a  resonant  cavity. 


Sample  Access 


Side  Veiw  Axial  View 

Magnetic  Fields  Electric  Fields 


Figure  5.  Schematic  showing  magnetic  fields  in  a  side  view  and  electric  fields  in  an 
axial  view  within  a  cylindrical  resonant  cavity. 


Microwaves  are  reflected  from  the  cavity  whenever  the  resonant  frequency  or 
the  Q  of  the  cavity  is  changed  by  absorption  of  microwaves  in  the  sample  under 
study.  Reflected  microwaves  go  back  to  the  circulator  which  guides  these  microwaves 
to  the  detector  diode.  The  detector  is  a  Schottky  diode  that  generates  a  signal  when 
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microwaves  are  incident.  This  detector  signal  is  measured  by  a  lock-in  amplifier  which 
requires  the  signal  to  be  modulated.  Since  the  intensity  of  the  reflected  microwaves 
is  a  function  of  magnetic  field,  modulation  of  the  detector  signal  is  accomplished  by 
adding  a  small  modulation  to  the  applied  magnetic  field.  This  modulates  the  intensity 
of  the  microwaves  incident  upon  the  detector  diode  creating  a  modulated  signal.  The 
small  modulation  in  the  applied  magnetic  field  is  generated  by  a  set  of  coils  that 
are  embedded  in  the  sides  of  the  microwave  cavity.  The  modulation  frequency  is 
typically  100  kHz.  The  lock-in  amplifier  mixes  and  averages  the  detector  signal  with 
a  single-frequency  reference  signal  to  produce  another  signal  that  is  proportional 
to  the  modulation  amplitude  of  the  detector  signal.  This  signal  is  at  a  maximum 
whenever  the  diode  signal  and  reference  signal  have  the  same  phase.  This  method 
filters  out  noise  that  has  a  frequency  or  phase  that  differs  from  that  of  the  reference 
signal.  This  “filtered”  signal  and  the  averaging  increase  the  signal-to-noise  ratio 
of  the  spectrometer.  Since  the  detector  signal  is  modulated  by  applying  a  small 
modulation  to  the  applied  magnetic  field  and  the  microwave  power  reflected  by  the 
microwave  cavity  is  a  function  of  applied  magnetic  field,  the  signal  produced  by  the 
lock-in  amplifier  is  proportional  to  the  first  derivative  (slope)  of  the  microwave  power 
reflected  from  the  microwave  cavity  (Figure  [6]). 

The  electron  paramagnetic  resonance  spectrometer  at  the  Air  Force  Institute 
of  Technology  is  capable  of  detecting  as  few  as  5  x  1010  unpaired  spins  at  10  K  [37], 
Knowing  this  sensitivity,  the  concentration  of  paramagnetic  defects  present  in  the 
sample  can  be  estimated  with 

(5  x  1010)(AH)2(SA)(AL)(T/10)(1/K)  (51) 

where  A B  is  the  line-width  of  the  spectrum  in  Gauss,  SN  is  the  signal  to  noise  ratio, 
NL  is  the  number  of  lines  associated  with  the  resonance,  T  is  the  temperature  in 
Kelvin,  and  V  is  the  volume  of  the  sample  in  cm3  03- 
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Figure  6.  Figure  illustrating  the  effects  of  modulating  the  magnetic  held  of  a 
magnetic-held-dependent  signal.  The  modulation  amplitude  of  a  magnetic-field- 
dependent  signal  is  proportional  to  the  slope  of  the  signal.  When  this  modulation  is 
measured  by  a  lock-in  amplifier,  the  output  is  proportional  to  the  first  derivative  of 
the  signal. 

3.2  Magnetic  Field  Measurements  and  Calibration 

The  applied  magnetic  held  from  the  electromagnet  is  measured  with  a  Hall 
probe  when  collecting  a  spectrum,  and  with  a  teslameter  and  calibration  sample  when 
performing  precise  measurements  of  the  magnetic  held  at  an  electron  paramagnetic 
resonance  line.  The  Hall  probe  has  a  quick  response  time,  but  is  prone  to  direct- 
current  drift,  since  it  relies  on  an  applied  voltage,  and  is  dependent  on  the  magnetic 
held  direction.  Therefore,  it  is  more  suitable  for  controlling  the  applied  magnetic  held 
when  sweeping  the  applied  magnetic  held  during  the  measurement  of  a  spectrum. 
The  teslameter  is  a  more  precise  instrument  since  it  relies  on  the  principles  of  nuclear 
magnetic  resonance,  but  has  a  slow  response  time.  It’s  measurement  of  the  magnetic 
held  is  independent  of  the  magnetic  held  direction.  Therefore,  it  is  more  suitable 
for  precise  measurement  of  the  magnetic  held  when  measuring  the  magnetic  held 
necessary  to  observe  an  electron  paramagnetic  resonance  line. 
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The  Hall  probe  makes  use  of  the  Hall  effect  to  measure  the  applied  magnetic 
held.  The  Hall  effect  is  the  creation  of  an  electric  potential  in  an  electrically  conduct¬ 
ing  material  due  to  the  presence  of  applied  electric  and  magnetic  fields.  This  electric 
potential  is  proportional  to  the  product  of  the  applied  electric  held,  applied  magnetic 
held,  and  the  sine  of  the  angle  between  them.  The  proportionality  constant  is  a 
function  of  several  variables  including  the  conductivity  and  shape  of  the  electrically 
conducting  material.  If  the  applied  electric  held  is  held  constant,  the  magnitude  of 
the  applied  magnetic  held  can  be  measured.  The  Hall  probe  is  mounted  on  one  side 
of  the  electromagnet’s  gap.  Because  of  this,  the  difference  between  the  magnetic  held 
measured  by  the  Hall  probe  and  the  magnetic  held  at  the  sample  can  be  on  the  order 
of  1  mT. 


After  completing  a  sweep,  the  Hall  probe  is  used  to  set  the  electromagnet  at  a 
magnetic  held  where  a  resonance  was  detected.  The  teslameter  is  then  used  to  make 
a  precise  measurement  of  the  magnetic  held.  The  teslameter  measures  the  transition 


energy  of  the  nuclear  magnetic  resonance  of  protons.  As  seen  in  Section  2.3,  this 
energy  is  proportional  to  the  magnitude  of  the  magnetic  held.  The  teslameter  probe 
is  in  the  form  of  a  long  rod.  Within  the  tip  of  the  rod  is  a  material  that  contains 
protons  and  a  radio  frequency  coil  for  excitation  and  detection  of  the  nuclear  magnetic 
resonance.  The  tip  of  this  rod  can  be  placed  next  to  the  cavity  to  get  a  closer 
measurement  of  the  magnetic  held  at  the  sample,  but  can  still  differ  on  the  order  of 
0.01  mT.  Because  of  this,  a  calibration  sample  with  a  known  electron  paramagnetic 
resonance  line  is  used  to  correct  the  magnetic  held  measurement. 


Assuming  that  the  teslameter  will  measure  a  magnetic  held  of  zero  when  the 
magnetic  held  at  the  sample  is  zero  and  that  the  correlation  is  linear,  the  calibrated 
magnetic  held  is 

hl>cal 

Qca  l  HB Bcal 

where  ucai  is  the  microwave  frequency  and  Bcai  is  the  magnetic  held  at  which  the 
calibration  sample’s  known  electron  paramagnetic  resonance  line  is  detected,  gcai  is 


the  g  factor  of  the  calibration  sample’s  known  electron  paramagnetic  resonance  line, 
and  B  is  the  uncorrected  magnetic  field.  Both  of  the  magnetic  fields  in  this  equation 
are  measured  with  the  teslameter,  and  with  the  calibration  sample  and  experimental 
sample  in  the  same  physical  setup.  The  calibration  sample  used  in  this  research  is  a 
single  crystal  of  magnesium  oxide  that  has  been  doped  with  chromium.  Chromium 
Cr3+  substitutes  for  magnesium  Mg2+  and  creates  a  paramagnetic  defect  that  has  an 
isotropic  g  matrix  with  a  principal  value  of  1.9800  [38] . 
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IV.  Nitrogen  Substitutional  Impurity  in  Zinc  Oxide  (ZnO) 

Nitrogen  substitutes  for  oxygen  in  zinc  oxide  and  forms  an  optically  active  acceptor. 
Figure  [7] shows  a  nitrogen  ion  substituting  for  an  oxygen  ion  in  the  zinc  oxide  structure 
(see  Appendix  0i.  The  nitrogen  has  four  nearest-neighbor  zinc  ions  (one  axial  and 
three  basal)  and  twelve  neighboring  oxygen  ions  (divided  into  three  sets  of  3,  6  and 
3  with  the  middle  set  of  6  sightly  more  distant  from  the  nitrogen  than  the  upper 
and  lower  sets  of  3).  In  n-type  material,  nitrogen  acceptors  are  normally  in  the 
singly  ionized  charge  state  with  no  unpaired  electrons  and  no  electron  paramagnetic 
resonance  signal.  However,  when  these  samples  are  illuminated  at  low  temperature 
with  near-band-edge  laser  light,  a  portion  of  the  singly  ionized  nitrogens  are  converted 
to  the  neutral  charge  state,  and  thus  become  paramagnetic.  Using  bulk  n-type  zinc 
oxide  crystals,  Carlos  et  al.  [35J  and  Garces  et  al.  [3D|  were  the  first  to  identify 
the  photoinduced  electron  paramagnetic  resonance  spectrum  from  neutral  nitrogen 
acceptors.  As  expected,  the  spectrum  is  uniquely  distinguished  by  a  well-resolved 
three-line  14N  hyperhne  pattern.  This  neutral  nitrogen  acceptor  spectrum  is  also  seen 
in  powder  and  polycrystalline  zinc  oxide  mm-  Recent  investigations  (42  >  331  03] 
show  that  prominent  forbidden  lines  appearing  in  the  electron  paramagnetic  resonance 
spectrum  are  due  to  a  14N  nuclear  electric  quadrupole  interaction. 

The  question  as  to  whether  an  isolated  nitrogen  is  a  deep  or  shallow  acceptor 
in  zinc  oxide  has  received  considerable  attention  in  recent  years.  Advanced  first- 
principles  calculations  [45, 1307,  08J  have  suggested  that  nitrogen  is  a  deep  acceptor 
with  an  ionization  energy  approaching  1.3  eV.  On  the  experimental  side,  the  wave¬ 
length  dependence  of  the  photo-conversion  of  singly  ionized  nitrogen  acceptors  to  the 
neutral  charge  state  has  an  onset  near  2.0  eV  and  provides  evidence  that  nitrogen  is 
a  deep  acceptor  in  zinc  oxide  m-  Additional  support  for  the  deep  acceptor  model  of 
nitrogen  comes  from  the  observation  of  a  broad  photoluminescence  band  peaking  near 
730  nm  in  nitrogen-doped  bulk  crystals  [3U].  In  contrast,  a  donor-acceptor  pair  peak 
is  observed  in  photoluminescence  studies  of  nitrogen-doped  zinc  oxide  epitaxial  layers 
and  is  assigned  to  a  shallow-donor  to  shallow-acceptor  transition  with  an  acceptor 
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Figure  7.  Schematic  of  the  crystal  structure  of  zinc  oxide  as  seen  along  the  [1010] 
direction.  A  nitrogen  substitutes  for  an  oxygen  and  has  four  neighboring  zinc  cations 
and  twelve  neighboring  oxygen  anions. 

binding  energy  of  around  160  meV  and  a  donor  binding  energy  near  60  meV  [50].  It 
is  not  known  if  the  participating  shallow  acceptor  is  an  isolated  nitrogen  or  a  larger 
complex  involving  one  or  more  nitrogen  ions.  The  important  point,  however,  is  that 
these  donor-acceptor  pair  results  [50]  show  shallow  acceptors  are  formed  in  zinc  oxide 
by  nitrogen  doping.  The  conflicting  nature  of  the  reports  thus  far  suggests  a  need  for 
additional  fundamental  studies  of  nitrogen-doped  zinc  oxide. 

In  this  dissertation,  electron  paramagnetic  resonance  spectroscopy  is  used  to 
investigate  the  67Zn  hyperhne  interactions  associated  with  the  neutral  nitrogen  ac¬ 
ceptor  in  zinc  oxide.  These  results  have  been  recently  published  in  the  Journal  of 
Applied  Physics  [51].  Sets  of  less  intense  6'Zn  (4.1%  natural  abundance  [3])  hyper¬ 
hne  lines  are  observed  around  each  of  the  primary  14N  (99.6%  natural  abundance 
[3])  hyperhne  lines  in  the  photoinduced  spectra  taken  at  5  K  (using  442  or  633  nm 
laser  light).  They  are  assigned  to  axial  and  basal  67Zn  nearest  neighbors.  Principal 
values  and  principal-axis  directions  of  the  6'Zn  hyperhne  matrices  are  obtained  from 
the  angular  dependence  of  the  electron  paramagnetic  resonance  spectra.  These  spin- 
Hamiltonian  parameters  provide  information  about  the  distribution  of  the  unpaired 
spin  at  the  zinc  neighbors.  Comparing  the  experimentally  determined  6'Zn  and  14N 
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Figure  8.  Electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  molecular 
nitrogen  acceptor  N2  in  zinc  oxide  with  the  magnetic  held  along  the  [0001]  direction. 
The  data  were  recorded  from  a  nitrogen-doped  sample  at  a  microwave  frequency  of 
9.3957  GHz  and  a  temperature  of  26  K. 

hyperhne  parameters  with  values  predicted  by  density  functional  theory  calculations 
is  expected  to  provide  a  check  on  the  validity  of  the  deep  acceptor  model  of  nitrogen 
substituting  for  oxygen  in  zinc  oxide  with  no  nearby  perturbations  [42,  521  G2]  • 

The  bulk  zinc  oxide  crystal  used  in  the  present  investigation  was  cut  from  a 
larger  boule  grown  at  Eagle-Picher  (Miami,  OK)  by  the  seeded  chemical  vapor  trans¬ 
port  method  with  nitrogen  added  to  the  gas  stream.  Material  from  this  boule  was 
used  in  a  previous  study  [51]  where  the  molecular  nitrogen  acceptor  was  studied.  The 
electron  paramagnetic  resonance  spectrum  of  the  molecular  nitrogen  acceptor  is  de¬ 
tected  in  the  sample  used  in  this  investigation  and  is  shown  in  Figure  [8]  The  as-grown 
boule  appeared  yellow  and  Hall  measurements  at  room  temperature  verified  that  it 
was  n-type.  In  addition  to  the  usual  shallow  donor  signal,  electron  paramagnetic 
resonance  spectroscopy  showed  that  Fe3+  (3d5)  ions  were  present  before  illumination 
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Figure  9.  Electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  iron 
donor  Fe3+  (3d5)  in  zinc  oxide  with  the  magnetic  field  along  the  [0001]  direction. 
The  data  were  recorded  from  a  nitrogen-doped  sample  at  a  microwave  frequency  of 
9.3957  GHz  and  a  temperature  of  26  K.  The  spectrum  of  the  iron  donor  is  detected 
when  the  phase  of  the  lock-in  amplifier  is  90°  from  the  phase  necessary  to  detect  the 
spectra  of  the  nitrogen  acceptors  and  the  shallow  donor. 

at  low  temperature.  The  electron  paramagnetic  resonance  spectrum  of  the  iron  donor 
is  detected  in  the  sample  used  in  this  investigation  and  is  shown  in  Figure  [9}  These 
singly  ionized  iron  donors  act  as  electron  traps  when  electrons  are  pumped  away 
from  the  nitrogen  acceptors  by  near-band-edge  laser  illumination.  Dimensions  of  the 
sample  are  approximately  3x3x1  mm3. 

Helium-cadmium  (442  nm)  and  helium-neon  (633  nm)  lasers  were  separately 
used  to  successfully  convert  the  nitrogen  acceptors  from  the  singly  ionized  charge  state 
to  the  paramagnetic  neutral  charge  state  at  low  temperature.  The  laser  light  moves 
electrons  from  the  valence  band  or  singly  ionized  acceptors  to  shallow  donors  and 
deeper  iron  donors,  where  they  remain  stably  trapped  as  long  as  the  temperature  is 
kept  below  about  100  K.  From  the  intensities  of  the  electron  paramagnetic  resonance 
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signals,  the  largest  concentration  of  neutral  nitrogen  acceptors  photoinduced  in  the 
sample  is  estimated  to  be  approximately  8.0  x  1015  cm-3. 


4-1  Results 


Figure  [lO] shows  the  well-resolved  and  unique  sets  of  three  lines  that  characterize 
the  electron  paramagnetic  resonance  spectra  from  isolated  neutral  nitrogen  acceptors 
in  zinc  oxide  [39l  1401  I4T1  1421  1431  ITT] .  The  magnetic  held  is  along  the  [0001]  direction  in 
Figure  [To][  a)  and  is  in  the  basal  plane  along  the  [1010]  direction  in  Figure  [l0][b).  Both 
spectra  were  taken  at  5  K  after  exposure  at  this  temperature  to  442  nm  laser  light 
(633  nm  laser  light  produced  the  same  spectra).  At  temperatures  above  approximately 
9  K,  these  electron  paramagnetic  resonance  lines  broaden  beyond  recognition  because 
of  short  spin-lattice  relaxation  times.  The  set  of  three  lines  is  the  result  of  a  hyperhne 
interaction  between  an  unpaired  electron  and  a  14N  nucleus  (99.6%  natural  abundance 
[3]  with  1  =  1  [4]). 


Well-resolved  hyperhne  lines  from  67Zn  nearest  neighbors  are  easily  observed 
in  the  electron  paramagnetic  resonance  spectrum  of  the  neutral  nitrogen  acceptor. 
These  6'Zn  nuclei  have  a  natural  abundance  of  4.1%  [3]  and  have  I  =  |  [3].  Because 
of  this  low  natural  abundance,  the  6'Zn  hyperhne  lines  are  much  smaller  in  intensity 


than  the  primary  14N  hyperhne  lines  in  Figure  10  (the  67Zn  hyperhne  lines  are  ap¬ 


proximately  140  times  smaller  than  the  14N  hyperhne  lines).  Figure  11  shows  these 
weaker  hyperhne  lines  in  an  electron  paramagnetic  resonance  spectrum  taken  at  5  K 
after  illuminating  the  sample  at  this  temperature  with  442  nm  laser  light.  The  mag¬ 
netic  held  is  along  the  [0001]  direction  and  the  microwave  frequency  is  9.3973  GHz. 
Compared  to  Figure  [To][ a)  the  vertical  scale  is  expanded  to  show  the  symmetrical 
sets  of  hyperhne  lines.  A  nitrogen  substituting  for  an  oxygen  has  one  zinc  neighbor 
along  the  c  axis  and  three  zinc  neighbors  forming  a  basal  plane  (see  Appendix  |B|. 
With  the  magnetic  held  along  the  [0001]  direction,  the  axial  zinc  neighbor  gives  rise 
to  a  nondegenerate  electron  paramagnetic  resonance  spectrum,  while  the  basal  zinc 
neighbors  give  rise  to  a  triply  degenerate  electron  paramagnetic  resonance  spectrum. 


34 


1 _ I _ t _ I _ t _ I _ I _ I _ I _ I _ I _ I _ t _ I _ I _ I _ I _ , _ I _ I _ t _ I _ t _ I _ I _ I _ 1 

333  335  337  339  341  343 

Magnetic  Field  (mT) 

Figure  10.  Electron  paramagnetic  resonance  spectra  of  the  neutral  nitrogen  accep¬ 
tor  N°  in  zinc  oxide  with  the  magnetic  held  (a)  along  the  [0001]  direction  and  (b)  in 
the  basal  plane  along  the  [1010]  direction.  The  spectra  were  recorded  at  a  frequency 
of  9.3958  GHz  and  a  temperature  of  5  K.  The  spectra  are  photoinduced  with  442  nm 
laser  light  from  a  cw  HeCd  laser. 


This  makes  the  basal  b7Zn  hyperhne  lines  in  the  electron  paramagnetic  resonance 


spectrum  three  times  more  intense  than  the  axial  67Zn  hyperhne  lines.  In  Figure  11 


the  separation  of  adjacent  axial  lines  is  1.3  mT  and  the  separation  of  adjacent  basal 


lines  is  0.7  mT.  Stick  diagrams  above  the  spectrum  in  Figure  [TT]  identify  the  axial 
(red)  and  basal  (blue)  6'Zn  hyperhne  lines  that  surround  each  of  the  three  off-scale 
14N  hyperhne  lines. 


The  electron  paramagnetic  resonance  lines  in  Figure  [l0](a)  shift  to  higher  held, 
but  do  not  split,  as  the  magnetic  held  is  rotated  from  the  [0001]  direction  to  the 
basal  plane.  Forbidden  transitions  arising  from  the  14N  nuclear  electric  quadrupole 
interaction  begin  to  appear  between  the  allowed  lines  when  the  magnetic  held  is 
approximately  60°  away  from  the  [0001]  direction  |44j.  These  forbidden  transitions 
are  most  intense  when  the  held  is  about  86°  from  the  [0001]  direction  (allowed  and 
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Figure  11.  6'Zn  hyperfine  lines  in  the  electron  paramagnetic  resonance  spectrum 
of  the  neutral  nitrogen  acceptor  with  the  magnetic  field  along  the  [0001]  direction. 
The  data  were  recorded  at  a  microwave  frequency  of  9.3973  GHz  and  a  temperature 
of  5  K.  The  spectrum  is  photoinduced  with  442  nm  laser  light  from  a  cw  HeCd 
laser.  Stick  diagrams  show  the  axial  (red)  and  basal  (blue)  6'Zn  hyperfine  lines  that 
surround  each  of  the  three  large  (off-scale)  14N  hyperfine  lines. 


forbidden  lines  have  similar  intensities  near  this  direction).  Two  forbidden  lines  are 
still  present,  but  less  noticeable,  when  the  field  is  in  the  basal  plane.  In  Figure [To|(b) , 
these  forbidden  transitions  are  observed  on  either  side  of  the  middle  allowed  14N 
hyperfine  line. 


In  Figure  11  there  are  three  well- resolved  axial  6'Zn  hyperfine  lines  in  the 
field  region  below  the  14N  hyperfine  line  at  333.6  mT.  These  three  hyperfine  lines 
move  closer  to  the  14N  hyperfine  line  as  the  magnetic  field  is  rotated  from  the  [0001] 
direction  toward  the  [1010]  direction  and  no  splitting  into  components  is  detected. 
The  same  behavior  occurs  when  the  field  is  rotated  from  the  [0001]  direction  toward 
the  [2110]  direction.  These  low-field  axial  67Zn  hyperfine  lines  are  accurately  identified 
and  measured  only  during  the  first  40°  of  rotation  from  the  [0001]  direction.  Beyond 
this  angle,  basal  67Zn  hyperfine  lines  begin  to  interfere.  The  axial  67Zn  hyperfine 
lines  extend  out  farther  from  the  low-field  14N  hyperfine  line  than  the  basal  67 Zn 
hyperfine  lines  when  the  magnetic  field  is  along  the  [0001]  direction  (as  seen  in  Figure 
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0-  However,  as  the  magnetic  field  is  rotated  from  the  [0001]  direction  toward  the 
basal  plane,  this  ordering  reverses  for  angles  beyond  40°.  For  angles  greater  than  40°, 
the  basal  6'Zn  hyperfine  lines  extend  farther  out  from  the  14N  hyperfine  line  than 
the  axial  6'Zn  hyperfine  lines.  In  the  basal  plane,  the  axial  6'Zn  hyperfine  lines  are 
completely  obscured  by  the  inner  basal  6'Zn  hyperfine  lines.  This  interference  from 
basal  6'Zn  hyperfine  lines  limits  the  range  where  the  axial  6'Zn  hyperfine  lines  can 
be  measured  (and  thus  accounts  for  the  lack  of  data  points  beyond  40°  in  Figure 


12).  Fortunately,  the  experimental  results  obtained  from  the  first  40°  of  rotation  are 
sufficient  to  determine  the  axial  67Zn  hyperfine  parameters,  as  described  in  Section 


4.2[  The  experimental  positions  of  these  axial  67Zn  hyperfine  lines  are  plotted  as 


discrete  points  in  Figure  [12} 


Figure  12.  Angular  dependence  of  the  axial  67Zn  hyperfine  lines.  Solid  and  dashed 
lines  are  calculated  at  9.395  GHz  using  the  spin  Hamiltonian  parameters  in  Table  [2] 
Black  dots  are  experimental  data  recorded  at  5  K  and  scaled  to  9.395  GHz.  Solid 
lines  represent  6'Zn  hyperfine  interactions  and  dashed  lines  represent  14N  hyperfine 
interactions. 


In  Figure  11  there  are  three  well- resolved  basal 
field  region  below  the  14N  hyperfine  line  at  333.6  mT. 


67Zn  hyperfine  lines  in  the 
The  lowest-field  basal  6'Zn 
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hyperfine  line  maintains  the  same  spacing  from  the  14N  hyperfine  line  for  the  first 
30°  of  rotation  of  the  magnetic  field  away  from  the  [0001]  direction  toward  the  [1010] 
direction.  Only  a  slight  broadening,  but  no  splitting  into  components  is  observed 
during  this  rotation  of  the  magnetic  field.  Continued  rotation  causes  this  lowest- 
fielcl  basal  6'Zn  hyperfine  line  to  split  into  two  components.  When  the  magnetic  field 
reaches  the  basal  plane,  the  two  components  that  emerged  from  the  original  hyperfine 
line  are  well  separated  and  easily  seen.  Unlike  the  case  of  the  axial  neighbor,  electron 
paramagnetic  resonance  spectra  taken  in  the  basal  plane  provide  detailed  information 


about  the  basal  67Zn  hyperfine  interactions.  Figure  13  shows  two  of  these  spectra 
obtained  at  5  K  and  a  microwave  frequency  of  9.3951  GHz  after  illuminating  the 
crystal  with  442  nm  laser  light.  Only  the  low-field  sides  of  the  spectra  are  shown 


in  Figure  13  because  the  shallow  donor  spectrum  interferes  with  the  high-field  sides 
of  the  spectra.  The  magnetic  field  is  along  the  [1010]  direction  in  Figure  [l3)(a)  and 
along  the  [2110]  direction  in  Figure  [l3|(b).  Interference  from  the  axial  67Zn  hyperfine 
lines  is  not  a  problem  in  the  basal  plane  because  the  basal  67Zn  hyperfine  lines  extend 
farther  from  the  14N  hyperfine  lines  than  the  axial  67Zn  hyperfine  lines.  Also,  in 
the  basal  plane,  the  basal  6'Zn  hyperfine  spacings  are  larger  than  the  14N  hyperfine 
spacings.  This  makes  it  easy  to  identify  individual  basal  67Zn  hyperfine  lines  on  the 


low-field  side  of  the  spectra  in  Figure  [13j  because  they  appear  as  sets  of  three  lines 
that  mimic  the  spacing  of  the  three  14N  hyperfine  lines.  The  lowest-field  three-line 


sets  are  identified  by  stick  diagrams  in  Figure  13 


When  the  magnetic  field  is  along  either  the  [1010]  or  [2110]  directions,  two  of 
the  three  zinc  ions  in  the  basal  plane  have  equivalent  hyperfine  interactions.  Thus, 
there  are  two  sets  of  basal  6'Zn  hyperfine  lines  for  these  directions  of  magnetic  field, 


with  one  set  twice  as  intense  as  the  other  set.  In  Figure  13,  basal  67Zn  hyperfine  lines 
from  a  doubly-degenerate  set  are  labeled  “double”  and  lines  from  a  nondegenerate 
set  are  labeled  “single.”  A  single  set  of  three  lines  (representing  the  interaction  with 
one  zinc  nucleus)  is  lowest  in  field  in  Figure  [l3|(a)  with  the  lowest  line  in  this  set  of 
three  located  at  339.7  mT.  The  lowest  line  of  the  double  set  is  located  at  340.1  mT 
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Figure  13.  Low-field  6'Zn  hyperfine  lines  in  the  electron  paramagnetic  resonance 
spectrum  of  the  neutral  nitrogen  acceptor  with  the  magnetic  field  in  the  basal  plane. 
The  data  were  recorded  at  a  microwave  frequency  of  9.3949  GHz  and  a  temperature 
of  5  K.  The  spectra  are  photoinduced  with  442  nm  laser  light  from  a  cw  HeCd  laser, 
(a)  The  magnetic  field  is  along  the  [1010]  direction,  (b)  The  magnetic  held  is  along 
the  [2110]  direction. 


in  Figure  [l3[a).  The  ordering  is  reversed  in  Figure  13  b)  for  the  other  orientation 
of  magnetic  held  in  the  basal  plane.  The  double  set  of  three  lines  (representing  the 
interaction  with  two  equivalent  zinc  nuclei)  is  now  lowest  in  magnetic  held  with  the 
lowest  line  in  this  set  of  three  located  at  339.8  mT.  The  lowest  line  of  the  single  set 
is  at  340.3  mT  in  Figure  [l3|(b) . 

During  a  30°  rotation  of  the  magnetic  held  from  the  [lOlO]  direction  to  the 
[2110]  direction  in  the  basal  plane,  the  single  set  in  Figure  |l3j]a)  moves  to  higher  held 


until  it  merges  with  another  set  to  form  the  double  set  in  Figure  13  b).  From  the 
perspective  of  a  single  basal  6'Zn  nucleus,  this  represents  rotation  angles  of  0°  —  30° 
from  the  [1010]  direction.  The  double  set  in  Figure  [T3](a)  splits  into  two  sets  with  one 
set  moving  to  lower  held  and  the  other  moving  to  higher  held.  The  set  that  travels 
to  lower  held  merges  with  the  single  set  in  Figure  [l3ta)  to  form  the  double  set  in 
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Figure  14.  Angular  dependence  of  the  basal  67Zn  hyperfine  lines.  Solid  and  dashed 
lines  are  calculated  at  9.395  GHz  using  the  spin  Hamiltonian  parameters  in  Table  [2] 
Black  dots  are  experimental  data  recorded  at  5  K  and  scaled  to  9.395  GHz.  Solid 
lines  represent  6'Zn  hyperfine  interactions  and  dashed  lines  represent  14N  hyperfine 
interactions. 


Figure  |l3|b).  From  the  perspective  of  a  single  basal  67Zn  nucleus,  this  represents 
rotation  angles  of  60°  —  30°  from  the  [1010]  direction.  The  set  that  travels  to  higher 
field  becomes  the  single  set  in  Figure  [l3](b).  From  the  perspective  of  a  single  basal 
6'Zn  nucleus,  this  represents  rotation  angles  of  60°  —  90°  from  the  [lOlO]  direction. 
Therefore,  data  for  90°  of  magnetic  field  rotation  for  a  single  basal  67Zn  hyperfine 
line  are  obtained  with  30°  of  actual  rotation  in  the  basal  plane.  The  experimental 
positions  of  these  low-held  basal  6'Zn  hyperfine  lines  are  plotted  as  discrete  points  in 


Figure  14 


4-2  Spin- Hamiltonian  Analysis 

Spin-Hamiltonian  parameters  describing  the  axial  and  basal  6'Zn  hyperfine  in¬ 
teraction  matrices  are  obtained  from  the  angular  dependence  of  the  67Zn  hyperfine 
lines.  The  width  of  each  electron  paramagnetic  resonance  line  is  about  0.1  mT.  This 
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gives  a  minimum  estimated  error  for  the  g  parameters  of  ±0.0006,  and  for  the  hy- 
perhne  interaction  parameters  of  ±0.6  MHz.  For  many  directions  of  magnetic  held, 
the  majority  of  the  hyperhne  lines  are  strongly  overlapping  and  it  is  difficult  to  make 
specific  assignments  and  extract  useful  information.  Fortunately,  as  described  in  Sec¬ 


tion  4.1,  there  are  angles  where  subsets  of  6'Zn  hyperhne  lines  on  the  low  held  side 
of  the  electron  paramagnetic  resonance  spectra  can  be  assigned  to  specihc  transitions 
and  used  to  determine  precise  values  of  the  parameters.  In  the  following  analysis  of 
the  nitrogen  electron  paramagnetic  resonance  spectra,  the  focus  is  placed  on  these 
low  held  6'Zn  hyperhne  lines  to  determine  the  6'Zn  hyperhne  interaction  parameters. 
Because  of  this  it  is  hrst  necessary  to  determine  the  g  parameters  and  14N  hyperhne 
interaction  parameters  from  the  angular  dependence  of  the  14N  hyperhne  lines. 

Parameters  for  the  g  matrix,  the  14N  hyperhne  interaction  matrix,  and  the  14N 
nuclear  electric  quadrupole  interaction  matrix  have  been  independently  determined 
by  several  groups  [3Hi  3D*  32,  33,  33]-  Their  findings  show  that  the  g  matrix,  14N 
hyperhne  interaction  matrix,  and  14N  nuclear  electric  quadrupole  interaction  matrix 
are  axial  with  the  unique  direction  along  the  [0001]  direction.  A  summary  of  their 
parameters  is  listed  in  Table  [T} 


Table  1.  Spin-Hamiltonian  parameters  for  the  14N  hyperhne  interaction  of  the 
neutral  nitrogen  acceptor  in  zinc  oxide  from  previous  studies.  The  hyperhne  and 
nuclear  electric  quadrupole  interaction  parameters  are  in  terms  of  h  and  are  in  units 
of  MHz.  _ 


Reference 

9± 

011 

Au 

^11 

e2qQ/h 

Carlos  et  al.  pITI] 

1.9633 

1.9953 

9.9 

82.5 

- 

Garces  et  al.  [40] 

1.9632 

1.9948 

9.5 

81.3 

- 

Gallino  et  al.  [42] 

1.9630 

1.9947 

8.6 

81.3 

2.65 

Evans  [43] 

1.9631 

1.9949 

7.95 

81.26 

2.95 

Stehr  et  al.\  T\ 

1.963 

1.995 

8.5 

81.1 

2.95 

The  spin  Hamiltonian 


Hs  =  pbS  •  g  •  B  ±  i  •  Q  •  i  ±  S  •  A  •  I®-  gN/aNi  •  B,  (53) 
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where  gjy  is  the  g  factor  for  14N  and  is  known  to  be  0.40375  [4],  describes  the  angular 
dependence  of  the  14N  hyperhne  lines  and  is  fitted  to  the  14N  hyperhne  data  presented 
in  Figures [12] and [14] and  included  in  Appendix[Ej  The  results  of  a  least-squares  fitting 


routine  are  g  parameters  of  1.9632  and  1.9951,  and  hyperhne  interaction  parameters 
of  8.0  MHz  and  81.1  MHz  for  principal  axes  in  the  basal  plane  and  along  the  [0001] 
direction,  respectively.  The  standard  deviation  of  the  frequency  data  from  the  fit 
is  121  kHz.  This  standard  deviation  translates  into  an  estimated  error  for  the  g 
parameters  of  ±0.00003,  and  for  the  hyperhne  interaction  parameters  of  ±0.2  MHz. 
These  estimates  are  within  the  minimum  estimates.  The  nuclear  electric  quadrupole 
interaction  parameters  of  Evans  [43j  were  used  in  the  fitting  routine.  Gallino  et  al.  S3 
determined  that  the  signs  of  the  hyperhne  interaction  parameters  are  opposite.  The 
hyperhne  interaction  parameter  along  the  [0001]  direction  is  assigned  a  positive  sign 
because  the  nuclear  magnetic  moment  for  14N  is  positive.  From  experiment  [43,  44], 
the  positions  of  forbidden  lines  between  the  middle  and  highest  held  14N  hyperhne 
lines  require  the  signs  of  the  hyperhne  interaction  parameter  in  the  basal  plane  and 
the  nuclear  electric  quadrupole  interaction  parameter  along  the  [0001]  direction  to 
have  opposite  signs.  These  parameters  are  used  in  the  6'Zn  spin  Hamiltonian  to 
extract  the  67Zn  hyperhne  interaction  parameters,  and  are  included  in  Table  [2} 

The  observation  of  a  maximum  hyperhne  splitting  along  the  [0001]  direction  in 
two  planes  of  rotation  for  the  axial  67 Zn  hyperhne  lines  establishes  this  direction  as 
a  principal  axis  direction.  These  hyperhne  lines  are  also  observed  to  be  very  close  to 
axial  with  the  unique  direction  along  the  [0001]  direction.  Using  these  principal-axis 
directions,  the  spin  Hamiltonian 


(54) 


n={N,Zn} 


where  gzn  is  the  g  factor  for  67Zn  and  is  known  to  be  0.35024  [4|,  describes  the 
angular  dependence  of  the  6'Zn  hyperhne  lines  and  is  htted  to  the  axial  67Zn  hyperhne 


data  presented  in  Figure  12  and  included  in  Appendix  [E]  The  spin  Hamiltonian  in 
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Equation  [54]  is  the  same  as  in  Equation  53  except  terms  have  been  added  for  the 
6'Zn  hyperhne  interactions.  The  results  of  a  least-squares  fitting  routine  results  in 
hyperhne  interaction  parameters  of  10.7  MHz  and  36.6  MHz  for  principal  axes  in  the 
basal  plane  and  along  the  [0001]  direction,  respectively.  The  standard  deviation  of 
the  frequency  data  from  the  fit  is  431  kHz.  This  standard  deviation  translates  into 
an  estimated  error  of  ±0.9  MHz.  This  estimate  is  slightly  larger  than  the  minimum 
estimate.  A  plot  of  the  axial  6'Zn  hyperhne  lines  is  shown  in  Figure  12  for  three 
planes  of  rotation.  These  parameters  are  included  in  Table  [2} 


Table  2.  Spin-Hamiltonian  parameters  for  the  neutral  nitrogen  acceptor  in  zinc 
oxide.  The  nuclear  electric  quadrupole  interaction  parameters  are  from  Evans  [43]. 
The  hyperhne  and  nuclear  electric  quadrupole  interaction  parameters  are  in  terms  of 
h  and  are  in  units  of  MHz.  The  estimated  error  in  the  g  parameters  is  ±0.0006  and 
in  the  hyperhne  interaction  parameters  is  ±1  MHz. 


Parameter  Principal  Value  Principal- Axis  Direction 

g  parameters 

g±  1.9632 

g\\  1.9951  [0001] 

Hyperhne  interaction  parameters  for  14N 

A±  -8.0 

,4 1|  81.1  [0001] 

Nuclear  electric  quadrupole  interaction  parameters  for  14N 
Q±  -1.5 

Q[:[  3.0  [0001] 

Hyperhne  interaction  parameters  for  axial  67Zn 
\A±\  10.7 

|A|||  36.6  [0001] 

Hyperhne  interaction  parameters  for  basal  67Zn 
[Ail  14.3  [1210] 

|A2|  18.3  [0001] 

|A3| _ 203 _ [1010] _ 


As  seen  in  Figures  |l3[a)  and  14,  the  lowest-held  (single)  basal  67 Zn  hyperhne 
line  has  a  turning  point  along  the  [1010]  direction.  This  suggests  that  a  principal  axis 
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of  the  basal  6lZn  hyperfine  interaction  matrix  is  along  this  direction.  Rotation  of  the 
magnetic  held  from  the  [1010]  direction  toward  the  [0001]  direction  provides  further 
verification  that  [1010]  is  a  principal-axis  direction.  In  this  second  plane  of  rotation, 
the  lowest-held  basal  6,Zn  hyperhne  line  also  has  a  turning  point  when  the  magnetic 
held  is  along  the  [1010]  direction.  The  crystal  lattice  of  zinc  oxide  has  C$v  symmetry. 
This  symmetry  has  a  six-fold  rotation  axis  along  the  [0001]  direction,  three  {1010} 
mirror  planes  and  three  {2110}  mirror  planes  in  zinc  oxide.  Due  to  this  symmetry,  any 
deviation  of  the  remaining  principal  axes  from  the  [0001]  or  [1210]  (normal  to  [1010]) 
directions  will  cause  additional  splitting  of  the  basal  6'Zn  hyperhne  lines.  Therefore, 
principal  axes  of  the  basal  67Zn  hyperhne  interaction  matrix  must  be  along  the  [0001] 
and  [1210]  directions,  because  additional  splitting  of  the  basal  67Zn  hyperhne  lines  is 
not  observed.  Using  these  principal-axis  directions,  the  spin  Hamiltonian  in  Equation 


(54)  is  htted  to  the  basal  b'Zn  hyperhne  data  presented  in  Figure  14  and  included  in 
Appendix  0  The  results  of  a  least-squares  fitting  routine  are  hyperhne  interaction 
parameters  of  14.3  MHz,  18.3  MHz,  and  20.3  MHz  for  principal  axes  along  the  [1210], 
[0001],  and  [1010]  directions,  respectively.  The  standard  deviation  of  the  frequency 
data  from  the  ht  is  198  kHz.  This  standard  deviation  translates  into  an  estimated 
error  of  ±0.4  MHz.  This  estimate  is  within  the  minimum  estimate.  A  plot  of  the 


basal  6'Zn  hyperhne  lines  is  shown  in  Figure  14  for  three  planes  of  rotation.  These 
parameters  are  included  in  Table  |2} 


4-3  Discussion 

Figure  [15]  is  a  schematic  representation  of  the  neutral  nitrogen  acceptor  in  zinc 
oxide.  The  axial  zinc  neighbor  and  one  basal  zinc  neighbor  are  shown.  Two  of  the 
principal-axis  directions  of  the  basal  67Zn  hyperhne  interaction  matrix  are  in  the 
basal  plane  with  A3  pointing  toward  the  nitrogen.  Having  A3  in  the  basal  plane 
suggests  that  the  nitrogen  ion  has  relaxed  along  the  [0001]  direction  (i.e.  the  nitrogen 
has  moved  away  from  its  axial  zinc  neighbor  to  a  new  equilibrium  position  that  is 
close  to  the  plane  containing  the  three  basal  zinc  neighbors).  This  relaxation  is 
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partially  driven  by  the  electrostatic  repulsion  between  the  “hole”  on  the  nitrogen  and 
the  positive  charge  of  the  axial  zinc  cation.  Density  functional  theory  calculations 
Phi  IE]  have  predicted  this  movement  away  from  the  axial  zinc  neighbor. 

A 

[0001] 
axial  Zn 

[1010] 

— — 


relaxation 

Figure  15.  Model  of  the  neutral  nitrogen  acceptor  in  zinc  oxide  showing  the  direc¬ 
tions  of  the  principal  axes  for  the  6'Zn  hyperhne  interaction  matrices.  These  three 
ions  are  in  the  (1210)  plane. 

Previous  magnetic  resonance  investigations  have  focused  on  the  14N  hyperhne 
interactions  ITf] .  In  the  present  research,  complete  sets  of  6'Zn  hyperhne  param¬ 
eters  have  been  obtained  for  the  one  axial  and  the  three  basal  zinc  neighbors  of  the 
neutral  nitrogen  acceptor.  These  6'Zn  results  establish  that  the  nitrogen  is  at  an 
oxygen  site  and  is  isolated  (i.e. ,  it  is  not  adjacent  to  another  defect).  Together,  the 
14N  and  6'Zn  results  provide  information  about  the  distribution  of  spin  density  in 
the  ground  state  of  this  neutral  defect.  The  67Zn  hyperhne  parameters,  especially, 
address  the  delocalization  of  the  unpaired  spin  and  can  be  used  as  an  experimental 
check  on  advanced  density  functional  theory  calculations. 
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The  14N  hyperfine  interaction  matrix  can  be  separated  into  isotropic  and  anisotropic 


parts, 


A(14n)  A 


1 

0 

0 

2 

5 

0 

0 

0 

1 

0 

+  p 

0 

2 

5 

0 

0 

0 

1 

0 

0 

4 

5 

(55) 


where  Morton  and  Preston  [55]  calculate  A  =  1811  MHz  for  an  electron  in  a  2s  orbital 
and  P  =  138.8  MHz  for  an  electron  in  a  2 p  orbital.  The  14N  hyperfine  interaction 
matrix  in  its  principal  axis  basis  is 


-8.0  0  0 

1  0  0 

-loo 

5 

A<“N>  _> 

0  -8.0  0 

=  21.7 

0  1  0 

+  74.25 

0-^0 

0  0  81.1 

0  0  1 

0  0  I 

5 

.  (56) 


This  hyperfine  interaction  matrix  suggests  that  «  1%  of  the  electron  wavefunction 
is  in  a  nitrogen  2s  orbital  and  ~  53%  of  the  electron  wavefunction  is  in  a  nitrogen 
2 p  orbital.  Thus,  a  significant  portion  of  the  unpaired  spin  (~  54%)  is  localized  on 
the  nitrogen  ion  [39],  02] . 

The  above  analysis  implies  that  approximately  46%  of  the  unpaired  spin  is 
distributed  on  the  zinc  and  oxygen  ions  neighboring  the  nitrogen  ion.  Unfortu¬ 
nately,  hyperfine  lines  from  the  oxygen  ions  are  not  observed  because  the  170  nu¬ 
clei  have  a  natural  abundance  of  only  0.038%  [3].  Thus,  there  is  no  experimental 
information  available  about  the  amount  of  spin  density  on  the  oxygen  ions.  On  the 
other  hand,  the  6'Zn  hyperfine  interaction  matrices  do  provide  information  about 
the  unpaired  spin  density  at  the  zinc  cations.  The  isotropic  parts  of  the  67Zn  ma¬ 
trices  (19.3  MHz  for  the  axial  neighbor  and  17.63  MHz  for  the  basal  neighbors)  indi¬ 
cate  that  the  unpaired  spin  has  only  a  small  contribution  from  4s  orbitals  (less  than 
1%)  on  the  zinc  neighbors  p)5].  Larger  contributions  to  the  unpaired  spin  density 
at  the  zinc  neighbors  from  the  3 d  and  4 p  orbitals  are  suggested  by  the  anisotropic 
parts  of  the  6'Zn  matrices.  Interpreting  these  anisotropic  results  in  terms  of  3d  and 
4 p  orbitals  is,  however,  challenging.  Depending  on  the  choice  made  for  the  signs 
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of  the  6lZn  principal  values,  the  results  of  Morton  and  Preston  [55]  indicate  that  a 
significant  portion  of  the  spin  density  (perhaps  5%  or  more)  may  be  found  on  each 
of  the  neighboring  zinc  ions.  Advanced  density  functional  theory  quantum  modeling 
methods  appear  to  be  the  best  approach  to  extract  electronic  information  from  these 
anisotropic  experimental  67Zn  results. 
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V.  Silicon  Interstitial  Impurity  in  Titanium  Dioxide  (Ti02) 

A  silicon-related  defect  is  identified  in  titanium  dioxide  through  observations  of  4'Ti 
and  49Ti  hyperhne  lines  with  the  magnetic  held  along  the  [001]  direction  and  the 
angular  dependence  of  29Si  hyperhne  lines.  These  observations  suggest  that  the  un¬ 
paired  electron  is  highly  localized  about  two  silicon  nuclei  with  one  silicon  located 
at  or  near  a  titanium  site  and  the  other  in  a  neighboring  tetrahedral  interstitial  site. 
The  unpaired  electron  is  primarily  located  at  the  interstitial  silicon  where  it  inter¬ 
acts  with  two  neighboring  titanium  nuclei.  The  results  of  an  isochronal  pulse  anneal 
experiment  with  an  as-grown  sample  and  observations  of  the  electron  paramagnetic 
resonance  spectrum  in  a  reduced  sample  lead  to  the  assignment  of  an  effective  1+ 
charge  state  for  the  observed  paramagnetic  state. 

The  rutile  titanium  dioxide  crystals  investigated  in  this  dissertation  were  grown 
by  the  Verneuil  process  and  purchased  from  Crystec  (Berlin,  Germany).  A  description 
of  the  rutile  structure  is  provided  in  Appendix  [C}  Samples  were  cut  from  a  10  x  10  x 
2  mm3  plate  into  3x3x2  mm3  pieces.  One  piece  was  left  in  its  as-grown  state. 
Another  piece  was  initially  used  in  a  previous  study  [56]  after  being  lithium-doped  by 
annealing  with  lithium  hydroxide  powder  for  6  hr  at  450  K  in  air.  A  third  piece  was 
initially  used  in  a  previous  study  EH  after  being  reduced  in  nitrogen  for  3  hr  at  650  K. 
The  as-grown  fully  oxidized  titanium  dioxide  crystals  include  several  unintentionally 
incorporated  defects.  These  defects  include  chromium,  copper,  iron,  and  vanadium 
substitutional  impurities,  and  oxygen  vacancies.  When  impurities  such  as  chromium 
and  iron  are  present  in  the  starting  materials,  they  cause  oxygen  vacancies  (for  charge 
compensation)  during  growth  due  to  their  preferred  trivalent  charge  states. 

The  electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  chromium 
acceptor  with  the  magnetic  field  along  the  [001]  direction  is  displayed  in  Figure 
and  shows  the  1  =  0  line  at  136.7  mT  and  the  characteristic  set  of  four  hyperhne 
lines  from  the  I  =  3  (3]  hyperhne  interaction  with  a  9.5%  naturally  abundant  [3]  53Cr 
nucleus.  The  hyperhne  splitting  of  the  outer  most  49Ti  hyperhne  lines  is  0.4  mT. 
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This  5=1  defect  was  characterized  by  Gerritsen  et  al.  m  as  a  substitutional  Cr3+ 
(3d3)  impurity  at  a  titanium  site.  Due  to  a  large  fine  interaction  and  the  limitations 
of  the  spectrometer  used  in  this  dissertation,  only  one  transition  is  observable. 


47Ti,  49Ti  I  1 

53Cr|  I 


Figure  16.  Electron  paramagnetic  resonance  spectrum  of  the  singly  ionized 
chromium  acceptor  Cr3+  (3d3)  in  titanium  dioxide  with  the  magnetic  field  along 
the  [001]  direction.  The  data  were  recorded  from  an  as-grown  sample  at  a  microwave 
frequency  of  9.5485  GHz  and  a  temperature  of  34  K. 

The  electron  paramagnetic  resonance  spectrum  of  the  neutral  copper  acceptor 
with  the  magnetic  field  along  the  [001]  direction  is  displayed  in  Figure  [17],  and  shows 
the  two  characteristic  sets  of  four  hyperfine  lines  from  the  I  —  |  [1]  hyperfine  in¬ 
teraction  with  a  69.2%  naturally  abundant  [3]  63Cu  nucleus  and  a  30.8%  naturally 
abundant  [3J  65Cu  nucleus.  This  S  —  |  defect  was  initially  described  by  Gerritsen 
and  Sabisky  [58j.  The  electron  paramagnetic  resonance  spectrum  shows  significantly 
intense  forbidden  transitions  due  to  a  strong  nuclear  electric  quadrupole  interac¬ 
tion.  The  hyperfine  and  nuclear  electric  quadrupole  interactions  were  characterized 
by  Ensign  et  al.  [59] .  Refined  hyperfine  and  nuclear  electric  quadrupole  interaction 
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parameters,  and  evidence  that  this  copper  related  defect  is  an  acceptor  with  a  Cn2+ 
(3d9)  ion  located  at  a  titanium  site  and  a  doubly  ionized  oxygen  vacancy  immediately 
neighboring  is  provided  by  Brant  et  al.  [BO]. 


forbidden 

63Cu  | 

65Cu  I 


Figure  17.  Electron  paramagnetic  resonance  spectrum  of  the  copper  Cu2+  (3d9) 
substitutional  impurity  next  to  an  oxygen  vacancy  in  titanium  dioxide  with  the  mag¬ 
netic  field  along  the  [001]  direction.  The  data  were  recorded  from  an  as-grown  sample 
at  a  microwave  frequency  of  9.5410  GHz  and  a  temperature  of  20  K. 


The  electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  iron  accep¬ 


tor  with  the  magnetic  field  along  the  [001]  direction  is  displayed  in  Figure  18,  and 
shows  the  1  =  0  line  at  83.7  mT  and  the  characteristic  set  of  two  hyperfine  lines  from 
the  I  —  \  [4]  hyperfine  interaction  with  a  2.1%  naturally  abundant  [3]  57Fe  nucleus. 
The  hyperfine  splitting  of  the  outer  most  49Ti  hyperfine  lines  is  0.2  mT.  This  S  =  | 
defect  was  characterized  by  Carter  and  Okaya  [61]  using  an  electron  paramagnetic 
resonance  spectroscopy  technique  that  is  significantly  different  than  the  technique 
employed  in  this  dissertation.  They  established  the  defect  as  a  substitutional  Fe3+ 


50 


(3d5)  impurity  at  a  titanium  site.  Due  to  a  large  fine  interaction  and  the  limitations 
of  the  spectrometer  used  in  this  dissertation,  only  one  transition  is  observable. 


47Ti 
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Figure  18.  Electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  iron 
acceptor  Fe3+  (3d5)  in  titanium  dioxide  with  the  magnetic  held  along  the  [001]  direc¬ 
tion.  The  data  were  recorded  from  an  as-grown  sample  at  a  microwave  frequency  of 
9.5484  GHz  and  a  temperature  of  34  K. 


The  electron  paramagnetic  resonance  spectrum  of  the  substitutional  vanadium 
isovalent  impurity  with  the  magnetic  held  along  the  [001]  direction  is  displayed  in 
Figure  19,  and  shows  the  characteristic  set  of  eight  hyperhne  lines  from  the  I  —  \  03 


hyperhne  interaction  with  a  100%  [3]  naturally  abundant  51V  nucleus.  This  S  —  \ 
defect  was  characterized  by  Gerritsen  and  Lewis  [6(2]  as  a  substitutional  V4+  (3d1) 
impurity  at  a  titanium  site. 
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Figure  19.  Electron  paramagnetic  resonance  spectrum  of  the  substitutional  vana¬ 
dium  isovalent  defect  V4+  (3d1)  in  titanium  dioxide  with  the  magnetic  field  along  the 
[001]  direction.  The  data  were  recorded  from  an  as-grown  sample  at  a  microwave 
frequency  of  9.5485  GHz  and  a  temperature  of  34  K. 

The  electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  oxygen 
vacancy  with  the  magnetic  field  along  the  [001]  direction  is  displayed  in  Figure 
and  shows  the  1  =  0  line  at  373.7  mT  and  two  characteristic  sets  of  hyperfine  lines. 
One  set  is  a  set  of  eight  hyperfine  lines  from  the  I  —  \  \ 3J  hyperfine  interaction  with 
a  5.4%  |3]  naturally  abundant  49Ti  nucleus.  The  other  set  is  a  set  of  six  hyperfine 
lines  from  the  I  =  |  [|]  hyperfine  interaction  with  a  7.4%  [3]  naturally  abundant  4 ‘  Ti 
nucleus.  The  hyperfine  splitting  of  the  outer  most  49Ti  hyperfine  lines  is  9.3  mT.  This 
S  =  \  defect  and  its  neutral  S  =  1  state  were  initially  described  by  Yang  et  al.  |63|  as 
unpaired  electrons  that  are  each  localized  about  a  single  neighboring  titanium  atom 
forming  Ti'1+  ions.  Of  the  three  titanium  nearest  neighbors  of  the  oxygen  vacancy, 
the  unpaired  electrons  are  localized  on  the  two  that  form  a  line  that  runs  along  the 
[001]  direction.  These  S  =  \  [64]  and  S  =  1  [65]  defects  were  later  characterized  by 
Brant  et  al. 
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Figure  20.  Electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  oxygen 
vacancy  Vq+  in  titanium  dioxide  with  the  magnetic  field  along  the  [001]  direction.  The 
data  were  recorded  from  an  as-grown  sample  at  a  microwave  frequency  of  9.5403  GHz 
and  a  temperature  of  30  K.  The  spectrum  is  photoinduced  with  442  nm  light  from  a 
cw  HeCd  laser. 


5.1  Results 


The  electron  paramagnetic  resonance  spectrum  shown  in  Figure  [21]  is  detected 
in  an  as-grown  Crystec  sample  of  rutile  titanium  dioxide  after  exposure  to  442  nm 
laser  light  with  the  magnetic  field  along  the  [001]  direction.  This  spectrum  is  observed 
in  the  Crystec  sample  after  lithium  diffusion,  the  Crystec  reduced  sample,  and  several 
samples  grown  by  the  Verneuil  process  at  the  University  of  Yamanashi,  Japan.  This 
spectrum  was  also  seen  before  exposure  to  442  nm  laser  light  in  an  as-grown  sample 
as  seen  in  Figure  [19]  but  is  less  intense.  The  central  1  =  0  line  at  352.1  mT  in  Figure 


21  is  created  by  a  single  unpaired  electron.  Surrounding  this  line,  are  six  hyperhne 


lines  overlapping  eight  hyperhne  lines  representing  the  hyperhne  interaction  of  the 
single  unpaired  electron  with  an  /  =  |  and  an  I  —  \  nucleus,  respectively.  There  are 
also  two  hyperhne  lines  representing  the  interaction  of  the  single  unpaired  electron 
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Figure  21.  Electron  paramagnetic  resonance  spectrum  of  the  silicon  defect 
[Sii  —  Sixi]1+  in  titanium  dioxide  with  the  magnetic  held  along  the  [001]  direction.  The 
data  were  recorded  from  an  as-grown  sample  at  a  microwave  frequency  of  9.5490  GHz 
and  a  temperature  of  40  K  after  exposing  the  sample  to  442  nm  light  from  a  cw  HeCd 
laser. 

with  an  I  =  \  nucleus,  and  multiple  hyperhne  lines  within  0.5  rnT  of  the  1  =  0 
line  representing  interactions  of  the  single  unpaired  electron  with  other  more  distant 
nuclei.  Two  of  these  latter  hyperhne  lines  are  identihed  as  part  of  the  six  and  eight¬ 
line  sets,  and  two  more  of  these  hyperhne  lines  are  attributed  to  an  interaction  of  the 
single  unpaired  electron  with  a  29Si  nucleus  [63].  The  six-line  pattern  is  identihed  as 
a  hyperhne  interaction  with  a  4 '  Ti  nucleus,  and  the  eight-line  pattern  is  identihed  as 
a  hyperhne  interaction  with  a  49Ti  nucleus.  These  assignments  are  made  because  4'Ti 
has  I  =  |  jl],  49Ti  has  I  —  \  \ 1]  and  the  intensity  ratio  of  the  high-held  47Ti/49Ti 
hyperhne  line  at  352.9  rnT  to  the  high-held  49Ti  hyperhne  line  at  353.2  mT  is  2.24. 
This  is  close  to  the  predicted  value  of  2.76  calculated  from  natural  abundances  [3j. 
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Comparing  the  intensities  of  the  4 '  Ti  and  49Ti  hyperhnc  lines  to  the  center  1  =  0 
line,  the  estimated  isotopic  abundance  of  the  4,Ti  isotope  is  12.19%  and  of  the  49Ti 
isotope  is  10.23%.  These  values  are  nearly  double  that  of  the  natural  abundances 
which  are  7.4%  for  4'Ti  and  5.4%  for  49Ti  0-  This  strongly  suggests  that  the  un¬ 
paired  electron  is  interacting  equally  with  two  titanium  nuclei.  Using  the  intensity 
of  the  hyperhne  line  at  352.2  mT,  the  estimated  isotopic  abundance  of  the  I  =  \ 
nucleus  previously  attributed  to  29Si  is  8.54%.  Using  the  intensity  of  the  hyperhne 
line  at  353.8  mT,  the  estimated  isotopic  abundance  of  the  unattributed  I  —  \  nucleus 
is  3.28%.  Based  on  known  natural  abundances  [2],  possible  candidates  for  the  previ¬ 
ously  attributed  nucleus  are  '7Se  (7.6%),  125Te  (7.1%),  and  29Si  (4.7%),  and  for  the 
unattributed  nucleus  are  29Si  (4.7%),  and  57Fe  (2.1%).  Since  this  is  an  undoped  sam¬ 
ple,  selenium  and  tellurium  are  not  expected  to  be  present,  but  iron  and  silicon  could 
possibly  have  been  present  in  the  crystal-growth  starting  materials.  This  is  because 
selenium  and  tellurium  are  some  of  the  least  abundant  elements  in  the  Earth’s  crust 
(66th  and  72nd,  respectively  [20]),  and  silicon  and  iron  are  some  of  the  most  abundant 
elements  in  the  Earth’s  crust  (2nd  and  4th,  respectively  [20]).  Iron  is  eliminated  from 
consideration  since  it  is  unlikely  to  form  a  defect  with  a  single  unpaired  electron,  and 
already  forms  a  stable  defect,  as  shown  in  Figure  [l8j  The  subsequent  angular  depen¬ 
dence  of  this  electron  paramagnetic  resonance  spectrum  in  Figure  [21]  establishes  it  as 
the  interaction  of  a  single  unpaired  electron  with  two  silicon  nuclei  and  two  titanium 
nuclei. 


The  angular  dependence  of  the  electron  paramagnetic  resonance  spectrum  is 


shown  in  Figure  [22]  The  data  were  recorded  without  442  nm  illumination  from  a 
rutile  titanium  dioxide  sample  annealed  in  lithium  hydroxide  [55].  The  intensity  of 
the  spectrum  from  this  sample  is  greatest  among  the  three  samples  considered  when 
not  exposed  to  442  nm  light.  Recording  data  with  no  exposure  to  442  nm  light 
eliminates  interference  from  the  oxygen  vacancy  which  has  electron  paramagnetic 
resonance  lines  in  the  same  region  of  magnetic  field  (300  —  400  mT). 
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Figure  22.  Angular  dependence  of  the  1  =  0  line  and  29Si  hyperfine  lines.  Solid 
lines  are  calculated  at  9.569  GHz  using  the  spin- Hamiltonian  parameters  in  Table  [3] 
Black  dots  are  experimental  data  recorded  at  40  K  and  scaled  to  9.569  GHz.  Red 
lines  represent  spectra  with  no  hyperfine  interactions.  Each  set  of  green  and  blue 
lines  represent  spectra  with  a  single  29Si  hyperfine  interaction. 

Upon  rotation  of  the  magnetic  held  from  the  [001]  direction  to  the  [100]  di¬ 
rection  in  the  (100)  plane,  there  is  no  significant  splitting  of  the  central  line  or  of 
the  29Si  hyperfine  lines.  Once  the  magnetic  held  is  rotated  oh  the  [001]  axis,  the 
titanium  hyperfine  lines  become  unidentifiable.  This  may  be  due  to  either  slightly 
inequivalent  hyperhne  interactions  with  each  titanium  nucleus,  the  nuclear  electric 
quadrupole  interaction  becoming  dominant  over  the  hyperhne  interaction  giving  rise 
to  forbidden  lines,  or  a  combination  of  both.  The  hyperhne  splitting  of  the  outermost 
49Ti  hyperhne  lines  is  2.3  mT  with  the  magnetic  held  along  the  [001]  direction.  The 
central  1  =  0  line  moves  to  a  lower  held  reaching  a  minimum  of  351.9  mT  along  the 
[100]  direction,  while  the  hyperhne  splitting  of  the  outer  29Si  hyperhne  lines  increases 
from  3.5  mT  to  8.6  mT,  reaching  a  maximum  along  the  [100]  direction.  The  inner 
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29Si  hyperfine  lines  get  “lost”  in  the  titanium  hyperfine  lines  until  the  magnetic  held 
has  been  rotated  about  60°  away  form  the  [001]  direction  where  they  emerge  from 
the  titanium  hyperfine  lines.  From  this  point,  the  hyperfine  splitting  of  the  inner  29Si 
hyperfine  lines  continues  to  increase,  reaching  a  maximum  of  2.8  rnT  along  the  [100] 
direction  from  an  initial  value  of  0.2  rnT  along  the  [001]  direction. 

Figure  [23]  shows  the  electron  paramagnetic  resonance  spectrum  with  the  mag¬ 
netic  held  along  the  [100]  direction.  Along  this  direction  the  29Si  hyperfine  lines  are 
clearly  observed  as  two  pairs  of  equally  intense  lines.  From  the  intensities  of  the  29Si 
hyperfine  lines  in  this  spectrum,  the  estimated  isotopic  abundance  is  4.7%.  This  is 
essentially  equal  to  the  natural  abundance  of  29Si  [3],  and  thus  conhrming  the  identity 
of  these  I  —  \  nuclei. 
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Figure  23.  Electron  paramagnetic  resonance  spectrum  of  the  silicon  defect 
[Sii  —  Sixi]1+  in  titanium  dioxide  with  the  magnetic  held  along  [100].  The  data  were 
recorded  at  a  microwave  frequency  of  9.5637  GHz  and  a  temperature  of  40  K  from  a 
sample  that  was  annealed  in  lithium  hydroxide  for  6  hr. 

Upon  rotation  of  the  magnetic  held  from  [100]  to  [110]  in  the  (001)  plane,  the 
1  =  0  line  and  the  29Si  hyperfine  lines  split  into  two  spectra.  One  spectrum  moves 
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to  a  higher  field,  reaching  a  maximum  of  357.0  rnT  along  the  [110]  direction.  The 
hyperfine  splitting  of  the  outer  29Si  hyperfine  lines  decreases,  reaching  a  minimum  of 
3.3  rnT  along  the  [110]  direction.  The  hyperfine  splitting  of  the  inner  29Si  hyperfine 
lines  decreases  as  well  until  they  are  lost  in  the  titanium  hyperfine  lines  for  angles 
greater  than  30°  from  the  [100]  direction.  The  other  spectrum  moves  to  a  lower  field, 
reaching  a  minimum  of  347.8  mT  along  the  [110]  direction.  The  hyperfine  splitting 
of  both  sets  of  29Si  hyperfine  lines  increases,  reaching  a  maximum  of  3.8  mT  for  the 
inner  set  and  11.5  mT  for  the  outer  set  along  the  [110]  direction. 

Upon  rotation  of  the  magnetic  field  from  the  [110]  direction  to  the  [001]  direction 
in  the  (110)  plane,  the  high  field  spectrum  moves  to  lower  field  while  the  low  field 
spectrum  moves  to  higher  field  until  the  two  spectra  merge  into  one  when  the  magnetic 
field  is  along  the  [001]  direction.  The  outer  29Si  hyperfine  splitting  of  the  high  field 
spectrum  increases  slightly.  The  inner  29Si  hyperfine  lines  of  the  high  field  spectrum 
are  not  observed  during  this  rotation.  The  outer  29 Si  hyperfine  splitting  of  the  low 
field  spectrum  decreases.  The  inner  29Si  hyperfine  splitting  of  the  low  field  spectrum 
decreases  as  well  until  the  hyperfine  lines  are  lost  in  the  titanium  hyperfine  lines  for 
angles  greater  than  50°  from  the  [110]  direction. 

These  angular  dependence  data  establish  that  a  single  unpaired  electron  is  in¬ 
teracting  with  two  silicon  nuclei  and  are  sufficient  for  determining  the  g  matrix  and 
29Si  hyperfine  interaction  matrices,  but  the  change  in  ionization  state  of  the  defect  due 
to  photo-excitation  is  yet  to  be  determined  by  an  isochronal  pulse  anneal  experiment 
and  an  observation  of  the  spectrum  in  a  reduced  sample  of  titanium  dioxide. 


The  results  of  an  isochronal  pulse  anneal  experiment  are  presented  in  Figure 


24  Data  were  recorded  from  an  as-grown  Crystec  sample.  The  sample  was  exposed 


to  442  nm  light  from  a  cw  HeCd  laser  at  a  temperature  of  20  K  for  5  min  prior 
to  commencing  the  experiment.  The  intensities  of  the  Cr3+,  Cu2+,  Fe3+,  Vq+,  and 
Vq  spectra  were  monitored  during  the  experiment  along  with  the  intensity  of  the 
silicon-related  spectrum.  A  pulse  anneal  consisted  of  raising  the  temperature  in  the 


microwave  cavity  to  a  specified  temperature,  holding  at  that  temperature  for  30  s, 
then  lowering  back  down  to  the  monitoring  temperature  of  20  K  where  the  intensities 
of  the  previously  mentioned  spectra  were  once  again  measured.  The  sample  was  not 
exposed  to  light  between  pulse  anneals. 


Temperature  (K) 

Figure  24.  Results  of  an  isochronal  pulse  anneal  experiment  on  a  sample  of  as- 
grown  titanium  dioxide.  Black  dots  are  experimental  data  recorded  at  a  temperature 
of  20  K  with  the  magnetic  field  along  the  [001]  direction.  Data  connected  by  red  lines 
are  from  the  singly  ionized  oxygen  vacancy,  by  green  lines  are  from  the  neutral  oxygen 
vacancy,  and  by  blue  lines  are  from  the  singly  ionized  silicon  interstitial  complex. 

After  exposure  to  light  but  before  starting  the  anneal,  the  intensity  of  the  Cr3+ 
spectrum  decreased  by  33%,  the  Cu2+  spectrum  decreased  by  74%,  and  the  Fe3+ 
spectrum  decreased  by  30%.  The  oxygen- vacancy-related  spectra  became  observable 
with  the  Vq"  spectrum  64%  less  intense  than  the  Vq  spectrum.  The  intensity  of  the 
silicon-related  spectrum  increased  by  a  factor  of  11.  During  exposure,  electrons  are 
effectively  pumped  from  Cr3+  and  Fe3+  defects  to  Cu2+  and  Vq+  defects.  At  this 
point,  it  is  unclear  if  electrons  are  pumped  to  or  from  the  silicon- related  defect. 
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During  the  isochronal  anneal,  the  Cr3+  spectrum  increased  steadily  until  it 
reached  full  recovery  after  the  40  K  anneal.  The  Cu2+  and  Fe3+  spectra  remained 
relatively  flat  throughout  this  experiment  which  ended  after  a  50  K  anneal,  and  are 
not  expected  to  recover  until  after  an  anneal  above  100  K  when  electrons  are  released 
from  Cu1+  defects  to  be  trapped  once  again  at  Fe3+  defects  [60]. 

As  the  temperature  of  the  anneal  step  increased,  the  intensity  of  the  Vq  spec¬ 
trum  decreased  slightly  until  it  started  decreasing  rapidly  between  the  24  K  and  30  K 
anneals,  and  disappeared  after  the  32  K  anneal.  During  this  drop  in  intensity  of  the 
Vo  spectrum,  the  intensity  of  the  silicon-related  spectrum  decreased  rapidly  until 
it  started  to  level  off  after  the  32  K  anneal,  and  the  intensity  of  the  Vq+  spectrum 
increased  rapidly  until  it  reached  a  peak  after  the  32  K  anneal.  The  decrease  in  the 
intensity  of  the  silicon-related  spectrum  nearly  mirrored  the  increase  in  the  intensity 
of  the  Vq+  spectrum  during  this  stage  of  the  experiment.  Between  the  32  K  and  44  K 
anneals,  the  silicon-related  spectrum  and  the  Vq+  spectrum  increased  their  rate  of 
decay  until  after  the  44  K  anneal  where  they  both  recovered  to  their  pre-exposure 
intensities. 

These  results  suggest  that  for  anneals  between  20  K  and  32  K  electrons  are  being 
released  from  Vq  defects,  which  become  Vq+  defects,  and  returned  to  the  silicon- 
related  defect.  For  anneals  between  32  K  and  44  K,  electrons  are  being  released 
from  Vq+  defects  and  returned  to  the  silicon-related  defect  until  the  defects  reach 
their  pre-exposure  concentrations.  This  establishes  the  silicon-related  defect  as  a 
hole  trap.  Therefore  during  exposure  to  442  nm  light,  electrons  are  being  pumped 
away  from  this  defect.  This  assertion  is  supported  by  an  observation  of  the  spectrum 
without  exposure  to  442  nm  light  in  a  reduced  sample  of  titanium  dioxide,  which 
has  been  annealed  in  nitrogen  for  3  hr  at  650  K.  The  intensity  of  the  spectrum  is 
70%  less  in  the  reduced  sample  compared  to  the  intensity  of  the  spectrum  in  the  as- 
grown  sample.  Since  reduction  decreases  the  intensity  of  the  electron  paramagnetic 
resonance  spectrum,  the  change  in  ionization  state  from  the  non-paramagnetic  state 
to  the  paramagnetic  state  is  positive  |5E,,  OHj- 
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5.2  Spin- Hamiltonian  Analysis 

The  angular  dependence  data  presented  in  Figure  [22]  and  included  in  Appendix 
[F]  are  sufficient  in  determining  the  full  set  of  parameters  necessary  for  describing  the 
g  matrix  and  29Si  hyperfine  interaction  matrices.  These  parameters  are  the  principal- 
axis  directions  and  the  principal  values  associated  with  them.  The  principal- axis 
directions  are  determined  through  observations  of  the  angular  dependence  data.  The 
principal  values  are  determined  by  fitting  the  angular  dependence  data  to  a  set  of 
spin  Hamiltonians.  The  width  of  each  electron  paramagnetic  resonance  line  is  ap¬ 
proximately  0.05  mT.  This  gives  a  minimum  estimated  error  for  the  g  parameters  of 
±0.0003,  and  for  the  hyperfine  interaction  parameters  of  ±0.3  MHz. 

Based  on  the  behavior  of  the  angular  dependence  data,  the  principal  axes  of  the 
g  matrix  and  both  29 Si  hyperfine  interaction  matrices  are  determined.  The  electron 
paramagnetic  resonance  spectrum  has  both  a  high  field  and  low  field  turning  point 
along  the  [110]  direction  in  two  planes  of  rotation.  This  establishes  that  the  g  matrix 
has  two  principal  axes  along  (110)  directions,  and  by  default  the  third  principal  axis 
must  be  along  the  [001]  direction.  The  same  can  be  said  for  the  outer  29Si  hyperfine 
lines,  where  the  hyperfine  splitting  reaches  both  a  maximum  and  a  minimum  along 
the  [110]  direction  in  two  planes  of  rotation.  This  establishes  that  the  outer  29Si 
hyperfine  interaction  matrix  has  two  principal  axes  along  (110)  directions,  and  by 
default  the  third  principal  axis  must  be  along  the  [001]  direction.  The  inner  29Si 
hyperfine  splitting  reaches  a  maximum  along  the  [110]  direction  in  two  planes  of 
rotation.  This  establishes  that  the  inner  29Si  hyperfine  interaction  matrix  has  at 
least  one  principal  axis  along  a  (110)  direction.  The  minimum  of  the  inner  29Si 
hyperfine  splitting  can  not  be  determined  because  the  inner  29 Si  hyperfine  lines  of 
the  high  field  spectrum  about  the  [110]  and  [001]  directions  and  in  the  (110)  plane 
are  not  observable.  A  second  principal  axis  is  determined  for  the  inner  29Si  hyperfine 
interaction  matrix  by  noticing  that  the  spectrum  does  not  split  in  the  (100)  plane  and 
that  the  crystal  has  symmetry.  This  symmetry  has  a  (001)  mirror  plane,  twofold 
rotation  axes  along  the  (110)  directions,  and  a  fourfold  rotation  axis  along  the  [001] 
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direction  in  titanium  dioxide.  Since  the  [110]  direction  has  already  been  established 
as  a  principal  axis,  the  (001)  mirror  plane  will  cause  the  spectrum  to  split  into  two 
components  if  a  principal  axis  is  not  along  the  [001]  direction.  Since  no  splitting  is 
observed  for  the  inner  29Si  hyperfine  lines  when  rotating  the  magnetic  held  in  the 
(100)  plane,  then  a  principal  axis  must  be  along  the  [001]  direction.  By  default, 
this  establishes  the  third  principal  axis  for  the  inner  29Si  hyperfine  interaction  matrix 
along  another  (110)  direction. 

The  principal  values  of  the  g  matrix  are  determined  by  fitting  spin  Hamiltonians 
to  angular  dependence  data  of  the  1  =  0  line  or  29Si  hyperfine  lines,  and  principal 
values  of  the  29Si  hyperfine  interaction  matrices  are  determined  by  fitting  a  spin 
Hamiltonian  to  their  respective  angular  dependence  data.  The  principal  values  of  the 
g  matrix  are  determined  by  fitting  the  spin  Hamiltonian 

Hs  =  ^BS-g-B  (57) 

to  the  angular  dependence  data  of  the  1  =  0  line.  This  line  and  spin  Hamiltonian 
represent  a  configuration  of  the  defect  that  has  only  1  =  0  nuclei  present.  The  results 
of  a  least-squares  fitting  are  g  parameters  of  1.9159,  1.9376,  and  1.9668  for  principal 
axes  along  the  [110],  [001],  and  [110]  directions,  respectively.  The  standard  deviation 
of  the  frequency  data  from  the  fit  is  867  kHz.  This  standard  deviation  translates 
into  an  estimated  error  of  ±0.0002.  This  estimate  is  within  the  minimum  estimate. 
The  principal  values  of  the  g  matrix  and  both  29 Si  hyperfine  interaction  matrices  are 
determined  by  fitting  the  spin  Hamiltonian 

Hs  =  /j,bS  ■  g  •  B  ±  S  •  A  •  i  —  gSifiNI  ■  B,  (58) 

where  gsi  is  the  g  factor  for  29Si  and  is  known  to  be  —1.11052  [1],  to  the  angular  de¬ 
pendence  data  of  the  29Si  hyperfine  lines.  Each  set  of  silicon  hyperfine  lines  and  the 
spin  Hamiltonian  represent  a  configuration  of  the  silicon-related  defect  with  a  single 
29Si  nucleus.  The  results  of  a  least-squares  fitting  to  the  inner  29Si  hyperfine  lines 
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are  g  parameters  of  1.9158,  1.9379,  and  1.9668,  and  hyperfine  interaction  parameters 
of  0.3  MHz,  0.8  MHz,  and  105.3  MHz  for  principal  axes  along  the  [110],  [001],  and 
[110]  directions,  respectively.  The  standard  deviation  of  the  frequency  data  from  the 
fit  is  986  kHz.  This  standard  deviation  translates  into  an  estimated  error  for  the 
g  parameters  of  ±0.0002,  and  for  the  hyperfine  interaction  parameters  of  ±2  MHz. 
The  g  parameter  error  estimate  is  within  the  minimum  estimate,  but  the  hyperhne- 
interaction-parameter  error  estimate  is  larger.  The  results  of  a  least-squares  fitting 
to  the  outer  29Si  hyperfine  lines  are  g  parameters  of  1.9159,  1.9376,  and  1.9668,  and 
hyperfine  interaction  parameters  of  91.4  MHz,  95.4  MHz,  and  316.4  MHz  for  principal 
axes  along  the  [110],  [001],  and  [110]  directions,  respectively.  The  standard  deviation 
of  the  frequency  data  from  the  fit  is  1.56  MHz.  This  standard  deviation  translates 
into  an  estimated  error  for  the  g  parameters  of  ±0.0003,  and  for  the  hyperfine  in¬ 
teraction  parameters  of  ±3  MHz.  The  g  parameter  error  estimate  is  equal  to  the 
minimum  estimate,  but  the  hyperhne-interaction-parameter  error  estimate  is  larger. 
A  consolidation  of  this  data  is  presented  in  Table  [3]  and  is  used  in  a  discussion  of  the 
defect  model. 

Table  3.  Spin-Hamiltonian  parameters  for  the  silicon  interstitial  in  titanium  diox¬ 
ide.  The  hyperfine  interaction  parameters  are  in  terms  of  h  and  are  in  units  of  MHz. 
The  estimated  error  in  the  g  parameters  is  ±0.0003  and  in  the  hyperfine  interaction 
parameters  is  ±3  MHz. 


Parameter 

Principal  Value 

Principal- Axis  Direction 

g  parameters 

9i 

1.9159 

[110] 

92 

1.9377 

[001] 

93 

1.9668 

[110] 

Hyperfine  interaction  parameters  for  29 Si 

lA’I 

0.3 

[110] 

lA’I 

0.8 

[001] 

lA’I 

105.3 

[110] 

l42,l 

91.4 

[110] 

l42,l 

95.4 

[001] 

l42)l 

316.4 

[110] 
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5.3  Discussion 


The  electron  paramagnetic  resonance  data  identify  a  defect  with  a  single  un¬ 
paired  electron  that  interacts  with  at  least  two  silicon  and  two  titanium  nuclei.  The 
29Si  hyperfine  parameters  give  information  about  the  distribution  of  the  electron  wave 
function.  The  intensity  of  the  47Ti  and  49Ti  hyperfine  lines  compared  to  the  intenisty 
of  the  1  =  0  line  along  with  the  known  structure  of  rutile  titanium  dioxide  (see  Ap¬ 
pendix  [C])  provide  information  about  the  location  of  the  silicon  atoms.  The  pulse 
anneal  experiment  suggests  that  this  defect  loses  an  electron  or  captures  a  hole  when 
exposed  to  442  nm  light  from  a  HeCd  laser.  The  spectrum  after  reduction  suggests 
that  the  change  in  ionization  state  from  the  non-paramagnetic  state  to  the  param¬ 
agnetic  state  is  positive.  All  of  this  information  supports  the  proposed  defect  model 
shown  in  Figure  [25]  In  this  model,  one  silicon  is  located  at  a  titanium  site  and  the 
other  silicon  is  located  at  a  neighboring  tetrahedral  interstitial  site.  The  silicon  ions 
are  stabilized  at  their  positions  by  a  double  bond. 


[110] 

— MHO] 


Figure  25.  Model  of  the  silicon  intersititial  impurity  in  titanium  dioxide.  Titanium 
atoms  are  in  turquoise,  oxygen  atoms  are  in  red,  and  silicon  atoms  are  in  gray.  The 
gold  bonds  are  those  normally  present  in  titanium  dioxide,  the  dashed  bonds  are  those 
that  would  exist  between  an  atom  at  a  tetrahedral  site  and  the  neighboring  oxygen 
atoms,  and  the  gray  bond  is  the  double  bond  between  the  two  silicon  atoms. 


Using  the  hyperfine  interaction  parameters  in  Table  [3]  along  with  the  atomic 
predictions  of  Morton  and  Preston  [55],  the  distribution  of  the  electron  wave  function 
among  the  silicon  nuclei  is  estimated.  For  a  29Si  atom,  the  hyperfine  interaction 
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matrix  can  be  broken  into  isotropic  and  anisotropic  parts, 


A'”si>  ->  A 


1 

0 

0 

2 

5 

0 

0 

0 

1 

0 

+  p 

0 

2 

5 

0 

0 

0 

1 

0 

0 

4 

5 

(59) 


where  Morton  and  Preston  [55]  calculate  A  =  —4594  MHz  for  an  electron  in  a  3s 
orbital  and  P  =  —285.5  MHz  for  an  electron  in  a  3 p  orbital.  Within  experimental 
error,  the  inner  silicon  hyperhne  interaction  matrix  in  its  principal-axes  basis  is 


0  0  0 

1 
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o 

-loo 

5 
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0  0  0 

=  -35 

0  1  0 

-87.5 

0  -I  0 
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0  0  1 

0  0  I 

5 

(60) 


where  the  sign  of  the  A^ }  hyperhne  interaction  parameter  is  chosen  to  be  negative, 
since  the  sign  of  the  g  factor  for  29Si  is  negative.  This  hyperhne  interaction  matrix 
suggests  that  ,  ~  1%  of  the  electron  wavefunction  is  in  a  silicon  3s  orbital 


-4549.4 


and 


-87.5 

-285.5 


31%  of  the  electron  wavefunction  is  in  a  silicon  3 p  orbital.  Within 


experimental  error,  the  outer  silicon  hyperhne  interaction  matrix  is 
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A<2>-> 
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=  -167.3 

0  1 

0 

-  185.83 
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(61) 


A  (2) 

where  the  sign  of  the  Hg  J  hyperhne  interaction  parameter  is  chosen  to  be  negative, 
since  the  sign  of  the  g  factor  for  29Si  is  negative.  The  signs  of  the  other  hyperhne 
interaction  parameters  are  chosen  to  be  the  same  because  choosing  signs  such  that 
A^  has  the  opposite  sign  of  A[2)  and  A^  results  in  P  =  —340.83  MHz  which  is 
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significantly  greater  than  the  calculated  P  parameter.  This  hyperfine  interaction 
matrix  suggests  that  4  ~  4%  of  the  electron  wavefunction  is  in  a  silicon  3s 
orbital  and  ~  65%  of  the  electron  wavefunction  is  in  a  silicon  3 p  orbital.  In 

total,  32%  of  the  electron  wavefunction  is  on  one  silicon  nuclei,  and  69%  is  on  the 
other  for  a  total  of  101%.  This  suggests  that  the  unpaired  electron  is  highly  localized 
on  this  pair  of  silicon  nuclei. 


The  intensity  of  the  47Ti/49Ti  hyperfine  lines  compared  to  the  intensity  of  the 


1  =  0  line  in  Figure  [21]  suggests  equivalent  hyperfine  interactions  with  two  titanium 
nuclei.  Based  on  this  and  and  the  tendency  of  silicon  to  form  a  4+  ion,  an  initial 
placement  for  a  silicon  is  at  a  titanium  site  which  has  two  titanium  nearest  neighbors 
along  the  [001]  direction  at  a  distance  of  2.959  A.  As  seen  with  the  chromium  (Figure 
16)  and  iron  (Figure  [l8|  substitutional  impurities,  this  placement  produces  hyperfine 


splittings  of  the  outer  most  i9Ti  hyperfine  lines  that  are  less  than  0.5  mT.  This 
is  more  consistent  with  the  set  of  hyperfine  lines  around  the  1  =  0  line  in  Figure 
The  47Ti/49Ti  hyperfine  pattern  is  more  consistent  with  the  47Ti/49Ti  hyperfine 
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pattern  of  the  oxygen  vacancy  (Figure  [20|  but  has  a  quarter  of  the  hyperfine  splitting. 
Silicon  is  not  likely  to  substitute  for  an  oxygen,  but  a  similar  site  is  the  tetrahedral 
interstitial  site.  This  site  is  surrounded  by  oxygen  anions  at  a  distance  of  1.7  A.  This 
is  a  reasonable  site  for  silicon  because  the  silicon  cations  in  the  various  polymorphs  of 
silicon  dioxide  are  tetrahedrally  coordinated  with  oxygen  anions  at  a  typical  distance 
of  1.6  A  p3Bj.  The  tetrahedral  interstitial  sites  like  the  oxygen  anions  have  three 
nearest-neighbor  titanium  cations.  Two  of  the  titanium  cations  are  located  at  a 
distance  of  1.871  A  and  the  other  is  located  at  a  distance  of  2.104  A.  A  paramagnetic 
silicon  impurity  at  this  site  would  likely  have  a  nondegenerate  4'Ti/49Ti  hyperfine 
pattern  in  addition  to  the  observed  doubly-degenerate  4iTi/49Ti  hyperfine  pattern. 
One  way  of  addressing  this  missing  hyperfine  pattern  is  to  replace  the  titanium  cation 
at  2.104  A  with  the  second  silicon.  With  a  silicon  at  a  tetrahedral  interstitial  site  and 
another  silicon  at  this  titanium  site,  all  of  the  features  of  the  electron  paramagnetic 


resonance  spectrum  in  Figure  21  are  explained.  This  arrangement  also  fits  with  the 
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principal-axis  directions  of  the  g  and  29Si  hyperfine  interaction  matrices,  because  the 
axis  formed  by  the  two  silicons  is  along  the  [110]  direction.  This  model  becomes  even 
more  intriguing  by  noticing  that  the  silicon-to-silicon  distance  is  close  to  the  typical 
bond  length  (2.2  A)  of  the  double  silicon  bond  in  acyclic  disilcnes  m- 

All  of  this  information  taken  together  suggests  that  the  silicon-related  defect  is 
most  likely  a  deep  donor  with  an  effective  charge  of  1+  in  the  paramagnetic  state. 
One  possible  ionic  configuration  of  this  silicon-related  donor  in  its  nonparamagnetic 
neutral  state  is  a  neutral  silicon  atom  at  the  interstitial  site  and  a  Si4+  ion  at  the 
titanium  site.  In  the  paramagnetic  state,  the  captured  hole  is  likely  to  be  localized  on 
the  interstitial  silicon,  essentially  becoming  a  Si1+  ion  which  has  an  unpaired  electron 
in  a  p  orbital.  The  only  problem  with  this  model  is  that  there  is  no  apparent  reason 
why  this  configuration  would  be  stable.  There  is  no  Coulomb  attraction  between 
the  two  silicons,  because  a  Si4+  ion  at  a  titanium  site  has  an  effective  charge  of 
0  according  to  the  crystal  lattice.  This  lack  of  a  stabilizing  force  is  addressed  by 
treating  the  silicons  as  a  disilene  molecule.  This  configuration  of  the  silicon-related 
donor  in  its  nonparamagnetic  neutral  state  has  the  silicons  sharing  a  double  bond. 
This  gives  each  silicon  an  effective  charge  state  of  2+.  In  this  model  the  silicon  atoms 
are  stabilized  in  their  positions  by  the  double  bond  and  a  Coulomb  attraction.  The 
Coulomb  attraction  arises  from  an  effective  charge  of  2+  at  the  interstitial  silicon 
and  an  effective  charge  of  2—  at  the  substitutional  silicon.  For  this  configuration, 
the  paramagnetic  state  likely  arises  from  an  unpaired  electron  in  a  bonding  orbital 
between  the  two  silicons.  Due  to  the  effective  charges  of  the  silicons,  this  unpaired 
electron  is  likely  to  be  repelled  from  the  substitutional  silicon  ion  and  attracted  to  the 
interstitial  silicon  ion.  This  essentially  creates  a  Si1+  ion  at  the  interstitial  site  and 
a  Si4+  ion  at  the  substitutional  site.  These  models  are  somewhat  speculative  which 
is  why  advanced  density  functional  theory  calculations  could  help  confirm  or  exclude 
these  models  as  an  explanation  for  the  observed  features  and  angular  dependence  of 
the  silicon-related  electron  paramagnetic  resonance  spectrum. 
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VI.  Sulfur  Vacancy  in  Stannous  Hexathiohypodiphosphate 

(Sn2P2Sc) 

Previous  electron  paramagnetic  resonance  studies  of  stannous  hexathiohypodiphos- 
phate  (SPS)  have  identified  substitutional  impurities  at  the  Sn2+  site  and  also  an 
intrinsic  defect  at  the  Sn2+  site.  The  substitutional  impurities  are  antimony  Sb2+ 
(5s25p1)  [68]  and  manganese  Mn2+  (3d5)  [[69],  and  the  intrinsic  defect  is  the  small  po- 
laron  Sn3+  CD].  The  intrinsic  small  polaron  was  previously  assigned  to  the  electron 


paramagnetic  resonance  spectrum  in  Figure  26  by  Ruediger  et  al.  mi,  but  the  as¬ 
signment  of  the  intrinsic  small  polaron  to  a  different  electron  paramagnetic  resonance 
spectum  by  Brant  et  al.  caa  has  placed  doubt  on  this  previous  assignment.  In  this 
dissertation,  the  electron  paramagnetic  resonance  spectrum  previously  assigned  to  the 
intrinsic  small  polaron  is  now  shown  to  be  due  to  the  singly  ionized  sulfur  vacancy. 
This  assignment  is  made  through  careful  observations  of  hyperfine  interactions  in  the 
angular  dependence  of  the  electron  paramagnetic  resonance  spectrum. 


6. 1  Results 

As  described  in  Appendix  [Dj  SPS  has  a  monoclinic  crystal  structure  that  only 
sightly  deviates  from  an  orthorhombic  crystal  structure.  Due  to  the  modest  angular 
dependence  of  the  electron  paramagnetic  resonance  spectrum  in  Figure  [26]  this  devi¬ 
ation  from  orthogonality  is  disregarded  and  the  directions  a,  b,  and  c  are  considered 
orthogonal. 

The  electron  paramagnetic  resonance  spectrum  shown  in  Figure  [26]  is  detected 
in  several  samples  of  SPS  after  exposure  to  442  nm  laser  light  at  a  temperature  of  45  K 
with  the  magnetic  field  along  the  b  direction.  After  photo-excitation,  the  spectrum 
is  stable  below  100  K.  The  sample  considered  in  the  current  study  has  dimensions 
of  1.5  x  2  x  3  mm3,  and  was  grown  by  the  Bridgman  process  at  the  Uzhgorod  State 
University,  Ukraine.  The  sample  was  provided  to  the  Air  Force  Institute  of  Technology 
by  Dr.  Dean  Evans  (Air  Force  Research  Laboratory,  Materials  and  Manufacturing 


117Sn,  119Sn 


117Sn,  119Sn 


Figure  26.  Electron  paramagnetic  resonance  spectrum  of  the  singly  ionized  sulfur 
vacancy  V^+  in  SPS  with  the  magnetic  held  along  the  b  direction.  The  data  were 
recorded  at  a  microwave  frequency  of  9.3896  GHz  and  a  temperature  of  45  K.  The 
spectrum  is  photoinduced  with  442  nm  laser  light  from  a  cw  HeCd  laser. 


Directorate).  The  electron  paramagnetic  resonance  spectrum  is  composed  of  a  central 
set  of  lines  at  354.1  mT  surrounded  by  four  sets  of  less  intense  lines.  The  inner  and 
outer  pairs  of  the  four  surrounding  sets  are  each  due  to  a  strong  hyperhne  interaction 
with  an  I  —  \  nucleus. 

A  close-up  of  the  central  set  and  the  high-held  hyperhne  sets  is  shown  in  Figures 


hgures  show  that  the  central  set  is  symmetric  and  composed  of  four  closely  spaced 
lines,  and  that  the  high-held  hyperhne  sets  are  also  symmetric  but  do  not  mirror 
the  structure  of  the  central  set,  as  expected.  With  strong  hyperhne  interactions, 
differences  in  the  nuclear  magnetic  moments  of  a  nucleus  among  its  isotopes  will 
cause  the  associated  hyperhne  lines  to  split.  The  observation  that  the  hyperhne  sets 


27  28j  and  [29]  for  magnetic  helds  along  the  a,  b ,  and  c  directions,  respectively.  These 
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Magnetic  Field  (mT) 

Figure  27.  Hyperfine  lines  in  the  electron  paramagnetic  resonance  spectrum  of  the 
singly  ionized  sulfur  vacancy  in  SPS  with  the  magnetic  held  along  the  a  direction. 
The  data  were  recorded  at  a  microwave  frequency  of  9.3959  GHz  and  a  temperature 
of  45  K.  The  spectrum  is  photoinduced  with  442  nm  laser  light  from  a  cw  HeCcl  laser. 

are  symmetric  indicates  that  the  interacting  nucleus  has  multiple  isotopes  that  have 
significant  (>  1%),  nearly  equal  abundances  and  nuclear  magnetic  moments  that  are 
evenly  distributed.  Nuclei  with  I  —  \  that  meet  these  requirements  include  107 Ag 
and  109 Ag  with  natural  abundances  of  51.8%  and  48.2%  respectively,  mCd  and  113Cd 
with  natural  abundances  of  12.8%  and  12.2%  respectively,  and  117Sn  and  119Sn  with 
natural  abundances  of  7.7%  and  8.6%  respectively  P  4],  Since  these  candidate  nuclei 
have  only  two  isotopes  to  consider,  the  outer  lines  in  a  hyperfine  set  are  associated 
with  the  outer  lines  in  the  central  set.  The  low- held  line  in  a  hyperfine  set  is  associated 
with  one  isotope,  and  the  high-held  line  in  a  hyperfine  set  is  associated  with  the  other 
isotope.  Since  hyperfine  interactions  are  proportional  to  the  nuclear  magnetic  moment 
of  the  nucleus,  the  difference  in  magnetic  held  between  a  line  in  a  hyperfine  set  and 
its  associated  line  in  the  central  set  is  proportional  to  the  nuclear  magnetic  moment 
of  the  associated  isotope.  Therefore,  the  ratio  of  the  nuclear  magnetic  moments  can 
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Figure  28.  Hyperfine  lines  in  the  electron  paramagnetic  resonance  spectrum  of  the 
singly  ionized  sulfur  vacancy  in  SPS  with  the  magnetic  field  along  the  b  direction. 
The  data  were  recorded  at  a  microwave  frequency  of  9.3942  GHz  and  a  temperature 
of  45  K.  The  spectrum  is  photoinduced  with  442  nm  laser  light  from  a  cw  HeCcl  laser. 


be  estimated.  In  addition,  an  estimate  of  the  abundance  of  each  isotope  can  be  made 


from  the  intensity  of  the  lines.  From  data  presented  in  Figures [27|[28|  and  29 ,  the  ratio 


of  the  nuclear  magnetic  moments  is  estimated  as  1.044  and  the  abundance  is  estimated 
as  5.6%.  The  ratio  of  the  known  nuclear  magnetic  moments  is  5jjqf§f  =  1-150  for 
10'Ag  and  109 Ag,  =  1.046  for  mCd  and  113Cd,  and  J^g949  =  1.046  for  117Sn 

and  119Sn  0.  Based  on  this  information,  the  most  likely  nuclei  responsible  for  the 


strong  hyperfine  interactions  are  those  of  tin. 


With  the  magnetic  held  along  the  c  direction,  Figure  [29]  shows  that  the  central 
set  is  centered  at  352.9  mT,  and  composed  of  four  equally  intense  hyperfine  lines 
that  are  split  by  a  central  interval  of  0.63  mT  and  by  two  outer  intervals  of  0.74  mT. 
Upon  rotation  from  the  c  direction  towards  the  a  direction  in  the  b  plane,  the  cen¬ 
tral  set  moves  to  a  lower  held  reaching  a  minimum  of  340.8  mT.  Also,  the  central 
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Figure  29.  Hyperfine  lines  in  the  electron  paramagnetic  resonance  spectrum  of  the 
singly  ionized  sulfur  vacancy  in  SPS  with  the  magnetic  held  along  the  c  direction. 
The  data  were  recorded  at  a  microwave  frequency  of  9.3898  GHz  and  a  temperature 
of  45  K.  The  spectrum  is  photoinduced  with  442  nm  laser  light  from  a  cw  HeCcl  laser. 


interval  decreases  reaching  an  apparent  minimum  of  0.35  rnT,  and  the  outer  interval 
decreases  slightly  reaching  an  apparent  minimum  of  0.68  mT.  Upon  rotation  from  the 
a  direction  towards  the  b  direction  in  the  c  plane,  the  central  set  moves  to  a  higher 
held  reaching  a  maximum  of  354.2  mT.  Also,  the  central  interval  increases  reaching 
an  apparent  maximum  of  0.49  mT,  and  the  outer  interval  increases  reaching  an  ap¬ 
parent  maximum  of  0.77  mT.  For  rotations  in  the  a  plane,  the  spectrum  becomes 
asymmetric  indicating  splitting  due  to  inequivalent  orientations  of  the  defect. 

The  angular  dependence  of  the  four  lines  of  the  central  set  is  shown  in  Figure  [30] 
and  suggests  a  hyperfine  interaction  with  two  I  —  \  nuclei  with  a  100%  abundance. 
Nuclei  with  I  —  \  and  a  nearly  100%  natural  abundance  incude  1H,  19F,  31P,  89Y, 
103Rh,  169Tm,  107 Ag  and  109 Ag,  and  203T1  and  205T1  [3j  3].  Of  these,  3H  and  19F 
could  possibly  be  incorporated  into  undoped  crystals  through  crystal-growth  starting 


72 


materials,  and  31P  is  certainly  present.  Since  SPS  is  composed  of  P2S6  molecules 
where  the  two  phosphorus  atoms  form  bonds  with  the  six  sulfur  atoms  (see  Appendix 
0  ,  two  31 P  nuclei  are  likely  responsible  for  the  four  hyperhne  lines  in  the  central  set. 


Figure  30.  Angular  dependence  of  the  31P  hyperhne  lines.  Solid  lines  are  calcu¬ 
lated  at  9.39  GHz  using  the  spin-Hamiltonian  parameters  in  Table  |4}  Black  dots  are 
experimental  data  recorded  at  45  K  and  scaled  to  9.39  GHz. 


The  additional  outer  lines  in  the  central  sets  in  Figures  |27| and  28  are  attributed 
to  hyperhne  interactions  with  the  remaining  four  tin  nuclei  that  surround  the  P2Sg 
molecule  (see  Appendix  [D]) .  This  attribution  is  supported  by  a  simulation  of  the  cen¬ 
tral  set  in  Figure  [28}  The  simulation  was  carried  out  by  using  the  best-ht  g  parameter 
and  31P  hyperhne  interaction  parameters.  The  width  of  the  electron  paramagnetic 
resonance  lines  and  the  value  of  a  single  hyperhne  interaction  parameter  for  four 
equally  interacting  tin  nuclei  were  adjusted  until  the  simulated  spectrum  closely  ht 
the  appearance  of  the  recorded  spectrum.  Figure  [31] shows  the  recorded  spectrum  and 
two  simulated  spectra.  Figure  [3l](a)  shows  a  simulation  that  uses  only  the  best-ht  g 
parameter  and  31P  hyperhne  interaction  parameters.  Figure  [3l|(c)  shows  a  simulation 
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that  uses  the  best-fit  g  parameter  and  31P  hyperfine  interaction  parameters  with  the 
additional  hyperfine  interaction  parameter  of  37  MHz  for  the  four  equally  interacting 
tin  nuclei.  Figure |3l|b)  shows  the  recorded  spectrum.  As  can  be  seen  in  Figure  [3T|t  he 
spectrum  in  Figure  [28]  is  described  well  by  the  best-fit  parameters  with  the  additional 
four  equally  interacting  tin  nuclei. 


348  351  354  357  360 

Magnetic  Field  (mT) 

Figure  31.  Simulations  of  the  central  set  with  the  magnetic  field  along  the  b  direc¬ 
tion  compared  to  recorded  data.  The  simulations  include  the  calculated  g  parameter 
and  31P  hyperfine  interaction  parameters  for  the  magnetic  held  along  the  b  direction, 
(a)  Simulation  that  includes  the  calculated  g  parameter  and  31P  hyperfine  interac¬ 
tion  parameters  only,  (b)  Recorded  data,  (c)  Simulation  that  includes  the  calculated 
g  parameter  and  31P  hyperfine  interaction  parameters,  and  a  hyperfine  interaction 
parameter  of  37  MHz  for  four  tin  nuclei. 


Only  the  average  angular-dependence  data  of  the  high-held  tin  hyperfine  sets 
and  central  set  were  recorded  in  the  b  and  c  planes.  This  was  because  identihcation 
of  individual  lines  in  each  hyperfine  set  was  difficult  due  to  the  width  of  the  hyperfine 
lines  and  the  slight  splitting  due  to  the  isotopes  of  tin,  and  overlap  with  other  lines 
in  the  low- held  region  made  observation  of  the  low-held  tin  hyperfine  sets  difficult. 
Angular-dependence  data  in  the  a  plane  was  not  recorded  because  splitting  of  the 
spectrum  in  this  plane  made  determination  of  the  centers  difficult.  These  angular- 
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dependence  data  are  shown  as  black  dots  in  Figure  [32]  The  tin  hyperfine  interactions 
appear  to  be  nearly  isotropic  with  an  average  splitting  from  the  center  set  of  15.4  mT 
for  the  inner  hyperfine  set  and  22.7  mT  for  the  outer  hyperfine  set. 


Figure  32.  Angular  dependence  of  the  central  set  and  tin  hyperfine  sets.  Solid 
lines  are  calculated  at  9.39  GHz  using  the  spin-Hamiltonian  parameters  in  Table  [5] 
The  red  line  is  the  calculated  angular  dependence  of  the  g  matrix,  and  the  green  and 
blue  lines  are  each  the  calculated  angular  dependence  of  a  hyperfine  interation  with 
a  single  tin  nucleus.  Black  dots  are  experimental  data  recorded  at  45  K  and  scaled 
to  9.39  GHz. 


6.2  Spin-Hamiltonian  Analysis 

The  angular-dependence  data  presented  in  Figure  [30]  and  included  in  Appendix 


G  are  sufficient  to  estimate  the  principal  values  and  principal-axis  directions  of  the 
g  matrix  and  31P  hyperfine  interaction  matrices.  The  angular-dependence  data  pre¬ 
sented  in  Figure[32]and  included  in  Appendix [G] are  sufficient  to  estimate  the  principal 
values  and  principal-axis  directions  of  the  g  matrix  and  117Sn  and  119Sn  hyperfine  in¬ 
teraction  matrices.  The  width  of  each  electron  paramagnetic  resonance  line  is  about 
0.4  mT.  This  gives  a  mimimun  estimated  error  for  the  g  parameters  of  ±0.002,  and 
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for  the  hyperfine  interaction  parameters  of  ±2  MHz.  Through  observations  of  the 
angular-dependence  data  in  Figure  [32j  one  principal-axis  direction  is  determined  for 
the  g  matrix.  Using  this  observation  and  making  some  reasonable  assumptions,  the 
principal  values  and  principal-axis  directions  of  the  g  matrix  and  31P  hyperfine  inter¬ 


action  matrices  are  estimated  by  fitting  the  angular-dependence  data  in  Figure  [30]  to 
a  spin  Hamiltonian.  Using  some  reasonable  assumptions  and  the  measured  ratio  of 
the  nuclear  magnetic  moments  of  117Sn  and  119Sn,  the  principal  values  and  principal- 
axis  directions  of  the  g  matrix,  and  11  'Sn  and  119Sn  hyperfine  interaction  matrices 


are  estimated  by  fitting  the  angular-dependence  data  in  Figure  32  to  another  spin 
Hamiltonian. 


Angular-dependence  data  in  Figure  [32]  show  that  the  spectrum  has  a  low-held 
turning  point  along  the  a  direction  in  two  planes  of  rotation.  This  establishes  one 
principal- axis  direction  for  the  g  matrix  as  the  a  direction.  Since  the  crystal  has 
Ci h  symmetry  with  a  single  mirror  plane  normal  to  the  b  direction,  the  asymmetry 
observed  in  the  spectrum  for  rotation  of  the  magnetic  held  in  the  a  plane  is  likely 
due  to  inequivalent  orientations  of  the  defect  caused  by  the  crystal  symmetry.  These 
inequivalent  orientations  are  created  by  having  the  principal  axes  rotated  by  an  angle 
about  the  a  direction  away  from  the  b  and  c  directions. 

The  g  and  31 P  hyperfine  interaction  parameters  are  estimated  by  observing 
that  the  angular-dependence  of  the  31P  hyperfine  splitting  appears  to  behave  in  the 
same  manner  as  the  center  position  of  the  spectrum.  Therefore,  it  is  assumed  that 
the  g  matrix  and  the  31P  hyperhne  interaction  matrices  share  the  same  principal 
axes.  Based  on  this  assumption,  when  rotating  the  magnetic  held  away  from  the  c 
direction  towards  the  b  direction  in  the  a  plane,  one  orientation  of  the  defect  will 
have  a  spectrum  that  expands  and  moves  to  lower  held,  while  the  other  orientation 
of  the  defect  will  have  a  spectrum  that  contracts  and  moves  to  higher  held.  Due  to 
this  behavior,  a  single  31P  hyperhne  line  should  emerge  from  the  central  set  at  the 
low-held  side.  After  rotating  the  magnetic  held  30°  away  from  the  c  direction  towards 
the  b  direction  in  the  a  plane,  a  hyperhne  line  appears  on  the  low-held  side  of  the 
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central  set.  This  hyperfine  line  remains  observable  until  the  magnetic  field  is  rotated 
by  60°  where  the  line  begins  to  merge  with  the  central  set.  This  line  has  about  half 
the  intensity  of  a  31 P  hyperfine  line  when  the  magnetic  field  is  along  the  c  direction 
supporting  the  previously  asserted  assumption.  A  spin  Hamiltonian 

2 

Hs  =  fiB  S  ■  g  ■  B  +  ]T  (S  •  A<">  •  i(n)  -  sPpvi(n)  •  b)  ,  (62) 

where  gp  is  the  g  factor  for  31P  and  is  known  to  be  2.2634  [3],  is  fitted  to  the  angular- 
dependence  data  in  Figure  [30}  The  results  of  a  least-squares  fitting  routine  where  the 
g  and  31P  hyperfine  interaction  matrices  are  assumed  to  have  the  same  principal-axis 
directions  are  g  parameters  of  1.9700,  1.8945  and  1.9011,  31P  hyperfine  interaction 
parameters  of  18.6  MHz,  20.5  MHz  and  19.3  MHz  for  one  hyperfine  interaction,  and 
31P  hyperfine  interaction  parameters  of  28.1  MHz,  33.4  MHz  and  36.3  MHz  for  the 
other  hyperfine  interaction.  The  principal-axis  directions  are  along  the  a,  b  ±  12° 
towards  ±c,  and  c  ±  12°  towards  =F b  directions,  respectively.  The  standard  deviation 
of  the  frequency  data  from  the  fit  is  329  kHz.  This  standard  deviation  translates  into 
an  estimated  error  for  the  g  parameters  of  0.00007,  and  for  the  hyperfine  interaction 
parameters  of  0.7  MHz.  These  estimates  are  well  within  the  minimum  estimates. 
These  parameters  are  listed  in  Table  [4]  for  the  g  matrix  and  31P  hyperfine  interaction 
matrices. 

The  g  and  117Sn  and  117Sn  hyperfine  interaction  parameters  are  estimated  by 
fitting  the  angular-dependence  data  in  Figure  [32]  to  a  set  of  spin  Hamiltonians.  The 
angular-dependence  data  in  Figure  [32]  are  measurements  of  the  center  of  the  cen¬ 
tral  and  tin  hyperfine  sets.  Since  the  high-field  tin  hyperfine  sets  were  the  only  tin 
hyperfine  sets  measured,  the  spin  Hamiltonian 

Hs  =  gBS  ■  g  B  (63) 
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Table  4.  31P  hyperfine  interaction  and  g  parameters  for  the  singly  ionized  sulfur 

vacancy  in  SPS.  The  hyperfine  interaction  parameters  are  in  terms  of  h  and  are  in 
units  of  MHz.  The  estimated  error  in  the  g  parameters  is  ±0.002  and  in  the  hyperfine 
interaction  parameters  is  ±2  MHz. 


Parameter 

Principal  Value 

Principal- Axis  Direction 

g  parameters 

,9i 

1.8945 

b  ±  12°  towards  ±  c 

92 

1.9011 

c  ±  12°  towards  =F  b 

93 

1.9700 

a 

Hyperfine  interaction  parameters  for  31 P 

i4‘i 

18.6 

a 

i4‘i 

19.3 

c  ±  12°  towards  =|=  b 

i4‘i 

20.5 

— b  ±  12°  towards  =F  c 

i42,i 

28.1 

a 

i42i 

33.4 

b  ±  12°  towards  ±  c 

l42)l 

36.3 

c  ±  12°  towards  =F  b 

is  first  fitted  to  the  angular-dependence  data  of  the  central  set  in  Figure  [32]  to  calcu¬ 
late  the  principal  values  of  the  g  matrix.  The  results  of  a  least-squares  fitting  routine, 
where  an  initial  guess  for  the  principal-axis  directions  is  along  the  a,  b ,  and  c  direc¬ 
tions,  are  g  parameters  of  1.9685,  1.8942,  and  1.8991,  and  principal-axis  directions 
within  2°  of  the  a,  b ,  and  c  directions,  respectively.  The  standard  deviation  of  the 
frequency  data  from  the  fit  is  4.07  MHz.  This  standard  deviation  translates  into 
an  estimated  error  of  ±0.0008.  This  estimate  is  well  within  the  minimum  estimate. 
These  parameters  are  used  in  the  spin  Hamiltonian 


Hs  =  g,BS  g  B  ±  S  A  i  -  gSn/j,Ni  •  B, 


(64) 


where  gsn  is  the  average  g  factor  of  117Sn  and  119Sn  and  is  calculated  from  known  data 
as  —2.046  [1],  to  estimate  the  principal  values  of  each  117Sn  and  117Sn  hyperfine  inter¬ 
action  matrix.  The  results  of  a  least-squares  fitting  routine  to  the  angular-dependence 
data  of  the  inner  tin  hyperfine  set  in  Figure  [32]  are  hyperfine  interaction  parameters 
of  877.3  MHz,  910.5  MHz,  and  822.2  MHz  for  principal-axis  directions  along  the  a, 
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b,  and  c  directions,  respectively.  The  standard  deviation  of  the  frequency  data  from 
the  fit  is  5.46  MHz.  This  standard  deviation  translates  into  an  estimated  error  of 
±10  MHz.  This  estimate  is  greater  than  the  minimum  estimate.  The  results  of  a 
least-squares  fitting  routine  to  the  angular  dependence  data  of  the  outer  tin  hyperhne 
set  in  Figure [32] are  hyperhne  interaction  parameters  of  1334.0  MHz,  1283.6  MHz,  and 
1291.1  MHz  for  principal- axis  directions  along  the  a,  b ,  and  c  directions,  respectively. 
The  standard  deviation  of  the  frequency  data  from  the  fit  is  7.08  MHz.  This  standard 
deviation  translates  into  an  estimated  error  of  ±10  MHz.  This  estimate  is  greater 
than  the  minimum  estimate.  With  these  average  parameters,  the  117Sn  and  117Sn 
hyperhne  interaction  parameters  for  each  isotope  are  calculated  using  the  measured 
ratio  of  the  nuclear  magnetic  moments.  These  parameters  are  listed  in  Table  [5]  for 
the  117Sn  and  117Sn  hyperhne  interaction  matrices. 


Table 
sulfur 
are  in 


5.  11  'Sn  and  119Sn  hyperhne  interaction  parameters  for  the  singly  ionized 

vacancy  in  SPS.  The  hyperhne  interaction  parameters  are  in  terms  of  h  and 


units  o 


MHz.  The  estimated  error  is  ±10  MHz. 


Parameter 

Principal  Value  Principal- Axis  Direction 

Hyperhne  interaction  parameters  for 

117Sn 

l41}l 

804.4 

c 

40 

858.3 

a 

40 

890.8 

b 

40 

1255.8 

b 

40 

1263.1 

c 

l42)l 

1305.1 

a 

Hyperhne  interaction  parameters  for 

119Sn 

m 

840.0 

c 

140 

896.3 

a 

140 

930.2 

b 

40 

1311.4 

b 

40 

1319.1 

c 

40 

1362.9 

a 
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6.3  Discussion 


The  electron  paramagnetic  resonance  data  identify  a  defect  with  a  single  un¬ 
paired  electron  that  interacts  with  at  least  two  tin  and  two  phosphorus  nuclei.  The 
31P,  and  119Sn  hyperhne  parameters  give  information  about  the  distribution  of  the 
electron  wave  function.  This  information  along  with  the  crystal  structure  of  SPS  sup¬ 
port  the  assignment  of  this  electron  paramagnetic  resonance  spectrum  to  the  singly 
ionized  sulfur  vacancy. 

Using  the  hyperhne  interaction  parameters  in  Tables  [4]  and  [5]  along  with  the 
atomic  predictions  of  Morton  and  Preston  [55],  the  distribution  of  the  electron  wave 
function  among  the  phosphorus  and  tin  cations  is  estimated.  Since  both  of  these 
elements  are  post  transition  elements,  the  anisotropic  part  of  the  hyperhne  interaction 
matrix  can  be  separated  into  a  term  that  accounts  for  p  orbital  contributions,  and  a 
term  that  accounts  for  additional  contributions.  This  results  in 
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where  Morton  and  Preston  [55J  calculate  A  =  13306  MHz  for  an  electron  in  a  31P 
3s  orbital,  P  =  917  MHz  for  an  electron  in  a  31P  3 p  orbital,  A  =  —43920  MHz  for 
an  electron  in  a  119Sn  5s  orbital,  and  P  =  —1831  MHz  for  an  electron  in  a  119Sn  5 p 
orbital. 
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The  31P  hyperfine  interaction  matrix  for  the  weaker  hyperhne  interaction  in  its 
principal-axes  basis  is 


18.6  0  0 
0  19.3  0 

0  0  20.5 


J 

1 

0 

0 

2 

5 

0 

0 

1 

0 

0 

19.46 

0 

1 

0 

+  1.2916 

0 

2 

5 

0 

-0.35 

0 

-1 

0 

0 

0 

1 

0 

0 

4 

5 

0 

0 

0 

(66) 


where  the  signs  of  the  principal  values  are  chosen  to  give  positive  signs  for  A  and 
P.  The  31P  hyperhne  interaction  matrix  for  the  stronger  hyperhne  interaction  in  its 
principal-axes  basis  is 


A(2^  — >■ 


28.1  0  0 
0  33.4  0 

0  0  36.3 


1  0  0 

-loo 

5 

1 

o 

o 

0  1  0 

+  4.625 

0  -jj  0 

5 

-2.65 

0-10 

0  0  1 

i 

TtflLO 

o 

O 

0  0  0 

(67) 


where  the  signs  of  the  principal  values  are  chosen  to  give  positive  signs  for  A  and  P.  In 
total,  the  31P  hyperhne  interaction  matrices  can  account  for  6  +  L291g]tL625  « 

1%  of  the  electron  wave  function  in  phosphorus  3s  and  3 p  orbitals. 
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The  119Sn  hyperfine  interaction  matrix  for  the  weaker  hyperfine  interaction  in 


its  principal-axes  basis  is 


A(1^ 


-840.0  0  0 

0  -896.3  0 

0  0  -930.2 


-  888.83 


1 

o 

o 

-|  0  0 

5 

1 

o 

o 

0  1  0 

-  51.7083 

0  -jj  0 

5 

+  28.15 

0-10 

0  0  1 

i 

^liO 

O 

O 

0  0  0 

(68) 


where  the  signs  of  the  principal  values  are  chosen  to  give  negative  signs  for  A  and 
P.  This  hyperfine  interaction  matrix  accounts  for  Z%S99q  ~  2%  of  the  electron  wave 
function  in  a  silicon  5s  orbital  and  "^gj333  ~  3%  of  the  electron  wave  function  in  a 
silicon  5 p  orbital.  The  119Sn  hyperfine  interaction  matrix  for  the  stronger  hyperfine 
interaction  in  its  principal-axes  basis  is 


A^2)  — y 


-1311.4 

0 

0 


0 

-1319.1 

0 


0 

0 

-1362.9 


-  1331.13 


1  0  0 

-loo 

5 

1  0  0 

0  1  0 

-  39.7083 

o  -|  0 

5 

+  3.85 

0-10 

1 - 

o 

o 

i 

■^|LO 

O 

o 

0  0  0 

(69) 


where  the  signs  of  the  principal  values  are  chosen  to  give  negative  signs  for  A  and 
P.  This  hyperfine  interaction  matrix  accounts  for  «  3%  of  the  electron  wave 

function  in  a  silicon  5s  orbital  and  ~  2%  of  the  electron  wave  function  in 

a  silicon  5 p  orbital.  In  total,  the  119Sn  hyperfine  interaction  matrices  account  for 
approximately  10%  of  the  electron  wave  function  in  silicon  5s  and  5 p  orbitals. 

The  previous  analysis  accounts  for  only  11%  of  the  electron  wave  function  on 
neighboring  phosphorus  and  tin  cations.  Unfortunately,  hyperfine  lines  from  the  sulfur 
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anions  are  not  observed  because  33S  has  a  low  natural  abundance  of  0.75%  [3j,  making 
any  33S  hyperfine  interactions  difficult  to  observe.  Setting  this  aside,  it  is  assumed 
that  the  electron  wave  function  is  sparingly  dispersed  among  the  neighboring  sulfur 
anions  as  well.  Therefore,  the  electron  wave  function  appears  to  be  highly  diffuse  like 
the  electron  wave  functions  of  many  shallow-donor  defects. 

Since  the  sample  is  undoped  and  apparently  only  interacts  with  native  atoms, 
the  defect  is  likely  due  to  an  intrinsic  defect.  An  intrinsic  donor  defect  that  can 
possibly  have  the  observed  hyperfine  interactions  is  the  sulfur  vacancy.  Since  the 
defect  has  only  one  unpaired  electron,  the  sulfur  vacancy  is  likely  in  the  singly  ionized 
state.  As  described  in  Appendix  [DJ  the  two  phosphorus  cations  and  six  sulfur  anions 
form  a  hypodiphosphate  anion.  Therefore,  any  sulfur  vacancy  is  in  close  proximity 
(within  3.4  A)  to  two  phosphorus  nuclei  as  seen  in  Table  [6}  Each  sulfur  anion  is  in 
close  proximity  (within  3.3  A)  to  at  least  two  tin  cations.  Two  sulfur  anions  have  an 
additional  tin  cation  at  about  3.3  A,  two  other  sulfur  anions  have  an  additional  tin 
cation  at  about  3.5  A,  and  the  remaining  sulfur  anions  have  an  additional  tin  cation 
at  about  3.8  A.  It  is  the  locations  of  the  latter  sulfur  anions  (labeled  1  and  2  in  Table 
[6|  that  are  the  most  likely  candidates  for  the  location  of  the  sulfur  vacancy.  In  Figure 
|33|  one  of  these  sulfur  anions  is  encircled  with  a  circle,  and  the  other  is  encircled  with 
a  square.  They  are  each  connected  with  their  nearest  neighbor  tin  cations  by  a  solid 
black  line.  Since  the  a  direction  is  a  principal  axis  direction  of  the  g  matrix  with  the 
least  amont  of  shift  from  the  free  electon  g  factor,  it  is  interesting  to  note  that  the 
axis  formed  by  these  sulfur  anions  is  nearly  along  the  a  direction. 
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Table  6.  Distances  to  three  nearest  tin  and  two  nearest  phosphorus  neighbors  of 
each  sulfur  anion  in  the  hypodiphosphate  anion  of  SPS.  The  distances  were  calculated 
from  the  crystal  structure  parameters  of  Dittmar  and  Schafer  ra.  and  are  in  units 
of  A.  _ 


Sulfur 

Neighbor 

(Distance) 

1 

P 

Sn 

Sn 

P 

Sn 

(2.015) 

(2.882) 

(3.015) 

(3.356) 

(3.801) 

2 

P 

Sn 

Sn 

P 

Sn 

(2.018) 

(2.951) 

(2.990) 

(3.414) 

(3.853) 

3 

P 

Sn 

Sn 

P 

Sn 

(2.026) 

(2.802) 

(3.307) 

(3.395) 

(3.525) 

4 

P 

Sn 

Sn 

P 

Sn 

(2.022) 

(2.828) 

(3.080) 

(3.318) 

(3.397) 

5 

P 

Sn 

Sn 

Sn 

P 

(2.036) 

(2.943) 

(3.053) 

(3.269) 

(3.357) 

6 

P 

Sn 

Sn 

Sn 

P 

(2.023) 

(2.876) 

(3.171) 

(3.296) 

(3.344) 

C 


(c) 


Figure  33.  Possible  locations  for  the  sulfur  vacancy  in  SPS.  One  location  is  outlined 
with  a  square  and  the  other  is  outlined  with  a  circle.  Solid  black  lines  are  drawn  from 
each  location  to  their  nearest-neighbor  tin  cations.  A  hypodiphosphate  unit  as  viewed 
along  the  (a)  a,  (b)  b,  and  (c)  c  directions  is  shown  with  its  six  nearest  neighbor  tin 
cations.  Phosphorus  cations  are  blue,  tin  cations  are  green,  and  sulfur  anions  are  red. 
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VII.  Summary  and  Conclusion 

Observations  of  hyperfine  interactions  were  used  to  further  characterize  or  identify 
point  defects  in  three  wide-band-gap  semiconductors  of  particular  commercial  and/or 
military  importance.  These  wide-band-gap  semiconductors  are  zinc  oxide,  rutile  ti¬ 
tanium  dioxide,  and  stannous  hexathiohypodiphosphate  (SPS).  The  neutral  nitrogen 
acceptor  in  zinc  oxide  was  further  characterized  while  the  silicon  interstitial  impurity 
in  rutile  titanium  dioxide  and  the  singly  ionized  sulfur  vacancy  in  SPS  were  identified. 

The  neutral  nitrogen  acceptor  in  zinc  oxide  was  further  characterized  through 
the  observation  of  67Zn  hyperfine  interactions  in  the  electron  paramagnetic  resonance 
spectrum.  A  study  of  the  angular  dependence  resulted  in  67Zn  hyperfine  interaction 
parameters  for  axial  and  basal  zinc  neighbors  of  the  nitrogen  substitutional  impurity. 
These  6'Zn  hyperfine  interaction  parameters  along  with  the  14N  hyperfine  interaction 
parameters  provide  more  information  about  the  physical  and  electronic  structure  of 
the  defect.  The  principal-axis  directions  of  the  basal  67Zn  hyperfine  interaction  matix 
suggest  that  the  nitrogen  relaxes  away  from  the  axial  zinc  neighbor  and  into  the  plane 
of  the  basal  zinc  neighbors.  A  better  understanding  of  the  electronic  wave  function 
provided  by  the  data  in  this  dissertation  is  expected  to  aid  in  the  further  development 
or  validation  of  density  functional  theory  models  for  this  important  acceptor  in  zinc 
oxide. 

The  silicon  interstitial  impurity  in  titanium  dioxide  was  identified  through  the 
observation  of  4'Ti,  49Ti,  and  29Si  hyperfine  interactions  in  the  electron  paramagnetic 
resonance  spectrum.  A  study  of  the  angular  dependence  resulted  in  g  parameters 
and  two  sets  of  29Si  hyperfine  interaction  parameters.  The  hyperfine  interaction 
parameters  suggest  that  the  single  unpaired  electron  is  tightly  bound  to  two  silicon 
ions.  A  pulse  anneal  experiment  was  performed  to  determine  the  thermal  stability 
of  paramagnetic  defects  formed  after  exposure  to  near-band-edge  light.  A  sample 
was  reduced  by  heating  in  flowing  nitrogen  to  determine  the  change  in  charge  state 
between  the  non-paramagnetic  and  paramagnetic  states  of  this  silicon- related  defect. 
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These  experiments  suggest  that  the  silicon-related  defect  is  likely  a  deep  donor.  The 
47Ti  and  49Ti  hyperhne  interactions  along  with  the  known  crystal  structure  suggest 
that  one  silicon  is  located  at  a  titanium  site  and  the  other  is  located  at  a  neighboring 
tetrahedral  interstitial  site.  The  distance  between  the  two  silicons  is  very  close  to 
the  bond  length  between  two  silicon  atoms  in  acyclic  disilene  molecules,  suggesting 
that  the  silicons  may  share  a  double  bond.  This  will  have  to  be  confirmed  by  other 
experimental  methods  or  through  advanced  density  functional  theory  modeling. 

The  singly  ionized  sulfur  vacancy  in  SPS  was  identified  through  the  observation 
of  117Sn,  119Sn,  and  31P  hyperhne  interactions  in  the  electron  paramagnetic  resonance 
spectrum.  A  study  of  the  angular  dependence  resulted  in  two  sets  of  11 7Sn  and  119Sn 
hyperhne  interaction  parameters  and  two  sets  of  31P  hyperhne  interaction  parameters. 
The  hyperhne  interaction  parameters  suggest  that  the  single  unpaired  electron  wave 
function  is  highly  diffuse.  These  data  along  with  the  known  crystal  structure  suggest 
that  the  defect  in  its  paramagnetic  state  is  a  singly  ionized  sulfur  vacancy. 

Suggestions  for  future  work  follow  from  the  research  described  in  this  disserta¬ 
tion.  Hyperhne  interactions  associated  with  the  neighboring  67Zn  nuclei  should  be 
investigated  for  the  nitrogen  molecule  in  zinc  oxide  [54].  This  effort  will  complement 
the  present  study  of  67Zn  hyperhne  iterations  for  the  isolated  nitrogen  acceptor  in 
zinc  oxide.  The  silicon  interstitial  impurity  in  rutile  titanium  dioxide  needs  to  be 
further  studied  using  infrared  absorption  techniques.  The  silicon  interstitial  impurity 
is  expected  to  give  rise  to  a  local  vibrational  mode  in  the  midinfrared  region.  This 
will  provide  further  evidence  that  the  silicon  atoms  form  a  bond  as  described  in  the 
proposed  model.  From  here,  annealing  additional  samples  at  high  temperature  will 
determine  dissociation  energies.  Future  work  in  SPS  will  involve  obtaining  additional 
samples  from  our  crystal  growth  partner  at  the  Uzhgorod  State  University  in  Ukraine. 
A  portion  of  these  samples  will  be  annealed  in  a  vacuum  to  enhance  the  concentration 
of  sulfur  vacancies  while  others  will  be  annealed  in  a  sulfur  atmosphere  to  eliminate 
sulfur  vacancies.  Verifying  this  ability  to  modify  the  concentration  of  sulfur  vacancies 


in  a  given  sample  will  allow  device  engineers  to  optimize  the  photorefractive  behavior 
of  the  crystals. 

Electron  paramagnetic  resonance  spectroscopy  is  an  excellent  technique  for  de¬ 
termining  the  physical  and  electronic  structure  of  point  defects  in  semiconductors  and 
insulators  as  long  as  the  defects  are  paramagnetic  and  have  an  appropriate  spin-lattice 
relaxation  time.  This  technique  becomes  even  more  powerful  if  hyperhne  interactions 
are  present  or  when  coupled  with  other  spectroscopic  techniques.  It  is  a  necessary 
component  in  any  robust  materials-characterization  capability  that  guides  materials 
scientists,  crystal  growers,  or  systems  designers  and  engineers  towards  the  develop¬ 
ment  of  mission  critical  components  and  systems. 
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Appendix  A.  Angular  Momentum  and  Angular  Momentum  Operators 

For  classical  mechanical  systems  [73],  angular  momentum  of  a  particle  is  defined  as 

L  =  r  x  p,  (70) 


where  r  is  the  position  vector  from  a  specified  origin,  and  p  is  the  momentum  vector. 
For  a  system  of  particles,  the  total  angular  momentum  is 

J  =  ^r«xpi,  (71) 

i 

and  the  center  of  mass  is  defined  as 

(72) 

i 

where  M  is  the  total  mass  of  the  particles,  and  m  is  the  mass  of  a  particle.  The 
position  vector  of  a  particle  can  be  defined  as  the  sum  of  a  position  vector  in  the 
center  of  mass  frame  and  the  position  vector  of  the  center  of  mass, 


r  =  r7  +  R 


(73) 


where  r7  is  the  position  vector  in  the  center  of  mass  frame.  Equation  (73 )  is  substituted 
into  Equation  @  to  get 

J  =  J^(r7i  +  R)  x  Pi.  (74) 


Expanding  this  equation  and  taking 


dr  d  . 
p  =  m—  =  m—  (r  +  R) 
v  dt  dtK  J 


(75) 


results  in 


x — /  dr'i  \ .  c/R  d  \ — .  c/R 

J  =  2^  r  i  X  mi~JT  +  m*r  7  x  -IT  +  R  x  -J,  2_s  m*r  i  +  R  X  M 


dt 


(76) 


The  first  term  is  the  sum  of  the  angular  momenta  of  each  particle  in  the  center  of 
mass  frame.  The  second  and  third  terms  are  null  vectors,  since 

<77) 

i 

is  the  position  of  the  center  of  mass  in  the  center  of  mass  frame.  The  fourth  term 
is  the  angular  momentum  of  the  center  of  mass.  This  equation  is  generalized  by 
grouping  particles  (i)  into  multiple  centers  of  mass  (j)  such  that 

J  =  £  x  +  £  (Ri x  M^)  •  <78> 

With  this  equation,  the  total  angular  momentum  is  viewed  as  the  sum  of  internal 
angular  momenta  from  the  first  term  and  the  sum  of  external  angular  momenta 
around  a  specified  origin  from  the  second  term. 

For  quantum  mechanical  systems  m,  the  classical  definition  defines  relations 
among  a  set  of  position  and  momentum  operators, 

J  =  R  x  P.  (79) 


For  this  expression  to  be  of  use,  it  is  broken  up  into  its  component  parts, 


Jx 

Jy 


R  yPZ  ~  RZPy, 

RZPX  —  RXPZ:  and 


Jz  Rx  Py 


RyPX. 


(80) 


These  operators  represent  measurements  of  Cartesian  components  of  an  angular  mo¬ 
mentum.  Another  useful  operator  is  one  that  represents  a  measurement  of  the  mag¬ 
nitude  of  the  angular  momentum, 


p  =  j.j  =  3i  +  %  +  pt. 


(81) 


With  the  canonical  commutation  relations, 


\Rii  Pj  / ,  [ Pi ,  Pj  0,  and  \Pii  Pj 


0  where  i,j&{x,y,z},  (82) 


the  angular  momentum  operators  have  the  following  commutation  relations. 


\Jxi  Jy\ 

\Jy>  Jz] 

\Jzi  J x\ 


=  ihjz, 

zHJxi 

=  ihJy ,  and 
=  0  where  i  G  {x,y,z}. 


(83) 


These  relations  indicate  that  the  total  angular  momentum  operator  shares  eigenstates 
with  each  of  its  component  operators  and  that  only  one  component  of  the  angular 
momentum  can  be  simultaneously  measured  with  the  total  angular  momentum.  Tra¬ 
ditionally,  the  2;  component  is  chosen  as  the  component  to  measure  simultaneously 
with  the  total  angular  momentum,  and  its  operator  is  called  an  angular  momentum 
projection  operator.  The  focus  now  is  to  solve  the  following  eigenvalue  equations. 


J 2  |  nj:nz)  =  A  j  |  nj,nz) ,  and  Jz  \nj}  nz)  =  \z  \  rij,  nz) .  (84) 

This  is  accomplished  by  introducing  the  raising  and  lowering  operators, 

J+  =  Jx  +  iJy,  and  J_  =  Jx  —  iJy,  (85) 

and  analyzing  the  following  commutation  relationships. 


\JZ1  </+]  foJ+1 

[Jz,  J_]  =  —  HJ-,  and  (86) 

[J2, ,/,]  =  0  where  i  E  {+,  — }. 
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The  analysis  m  gives  the  following  results. 


J 2  \J, mj)  =  h2J(J  +  1)  |  J, mj) , 
•4  rrij)  =  hmj  \  J,mj) , 


(87) 


J+  |  J,  mj)  =  hsJj(J  +  1)  -  +  1)  |  J,  mj  +  1)  ,  and 


J-  |  J,mj)  =  h\J J(  J  +  1)  -  mj{mj  -  1)  |J,mj  -  1) 


where  J  €  is  the  angular  momentum  quantum  number  and  mj  G 

{— J,  —  J  +  1, . . . ,  J  —  1,  J}  is  the  angular  momentum  projection  quantum  number. 

There  are  two  flavors  of  angular  momentum  in  quantum  mechanical  systems. 
The  first  is  orbital  angular  momentum,  and  is  analogous  to  classical  angular  mo¬ 
mentum.  This  angular  momentum  is  typically  annotated  in  quantum  mechanics  by 
L.  The  second  is  spin  angular  momentum,  and  has  no  classical  analogue.  It  is  an 
intrinsic  property  of  a  particle.  This  angular  momentum  is  typically  annotated  in 
quantum  mechanics  by  S  for  an  elementary  particle  and  by  I  for  a  nucleus.  Since 
angular  momentum  can  be  seen  as  a  sum  of  internal  and  external  angular  momenta, 
the  total  angular  momentum  is  annotated  in  quantum  mechanics  by  J  such  that 


J  =  5>+5>- 


(88) 


i  j 


Since  all  of  the  angular  momentum  operators  have  a  common  factor  of  h,  it  is  cus¬ 
tomary  to  place  angular  momentum  operators  in  terms  of  h  such  that 


J  HJ. 


(89) 
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Appendix  B.  Crystal  Structure  of  Zinc  Oxide 

Zinc  oxide  is  a  binary  metal  oxide  composed  of  equal  parts  zinc  and  oxygen  by  atomic 
composition,  and  has  a  hexagonal  structure  with  four  basis  sites  as  shown  in  Figure 


of  length  a  that  are  separated  from  each  other  by  120°,  and  a  vector  (c)  of  length  c 
that  is  normal  to  these  coplanar  vectors.  The  plane  formed  by  the  coplanar  vectors  is 
called  the  basal  plane,  and  the  axis  formed  by  the  normal  vector  is  called  the  c-axis. 
For  zinc  oxide,  the  value  of  a  is  3.250  A,  and  the  value  of  c  is  5.207  A  [75] .  Two  of 
the  basis  sites  are  composed  of  zinc  cations,  and  form  a  hexagonal  quasi-close-pack 
structure  with  one  basis  site  at  the  origin  and  the  other  at 


34  Hexagonal  structures  are  commonly  defined  by  three  coplanar  vectors  (ai,  02,  03) 


-ax  +  -ay  +  -c,  where  x,  y  G  {1,  2, 3},  and  xj L-y. 


(90) 


A  hexagonal  close-pack  structure  has  c/a  =  2  a/6/3  ~  1.633,  but  for  zinc  oxide  c/a  ~ 
1.602.  The  other  two  basis  sites  form  the  same  hexagonal  structure  except  they  are 
composed  of  oxygen  anions  that  are  displaced  along  the  c-axis  by  0.3825c  =  1.992  A 
m-  In  this  structure,  the  zinc  cations  are  coordinated  with  four  oxygen  anions  and 
vice  versa  forming  slightly  deformed  tetrahedra.  A  hexagonal  close-pack  structure 
forms  regular  tetrahedra.  Since  each  atom  is  assigned  to  two  basis  sites,  there  are 
two  sets  of  tetrahedra  that  differ  by  a  rotation  of  180°  in  the  basal  plane.  The 
central  atom  is  coordinated  with  an  atom  along  the  c-axis  and  three  atoms  that 
form  a  basal  plane.  The  bond  angle  between  the  c-axis  atom  and  each  of  the  basal 
plane  atoms  is  approximately  108°  while  the  bond  angle  between  the  basal  plane 
atoms  is  approximately  110°.  A  hexagonal  close-pack  structure  has  bond  angles  of 
approximately  109°.  The  bond  length  between  the  central  atom  and  each  basal  plane 
atom  is  1.974  A  which  is  slightly  shorter  than  the  bond  length  along  the  c-axis. 


Hexagonal  quasi-close-pack  structures  are  composed  of  tetrahedral  and  octahe¬ 


dral  interstitial  sites  as  shown  in  Figure  35  Tetrahedral  sites  are  located  a  distance 
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a. 


[1  120] 


Figure  34.  A  unit  cell  of  zinc  oxide  with  the  zinc  nuclei  colored  gray  and  the  oxygen 
nuclei  colored  red. 

of 


above  and  below  (along  the  c-axis)  an  atom.  The  four  surrounding  atoms  form  the 
vertices  of  a  deformed  regular  tetrahedron  and  are  equidistant  from  the  interstitial 
site.  This  distance  is  1.978  A  in  zinc  oxide.  Three  of  the  surrounding  atoms  form  an 
equilateral  triangle  of  edge  a  that  is  coplanar  with  the  basal  plane,  and  the  other  site 
forms  a  line  with  the  center  of  this  triangle  that  is  collinear  with  the  c-axis.  Depending 
on  the  c/a  ratio,  the  tetrahedron  is  deformed  along  the  c-axis  by  varying  the  spacing 
between  the  basal  triangle  and  the  other  vertex.  Since  the  c/a  ratio  for  zinc  oxide  is 
close  to  2\/6/3,  the  tetrahedra  are  nearly  regular  with  edge  a.  With  this  information, 
zinc  oxide  is  thought  of  as  a  hexagonal  close-pack  lattice  of  zinc  atoms  with  half  of 
the  tetrahedral  sites  occupied  by  oxygen  atoms  and  vice  versa.  Octahedral  sites  are 
located  a  distance  of  c/2  above  and  below  (along  the  c-axis)  close-pack  voids  in  the 
unit  cell.  The  six  surrounding  lattice  sites  form  the  vertices  of  a  deformed  regular 
octahedron  and  are  equidistant  from  the  interstitial  site  by  a  distance  of 


(92) 
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This  distance  is  2.284  A  in  zinc  oxide.  Three  of  the  surrounding  lattice  sites  form 
an  equilateral  triangle  of  edge  a  that  is  coplanar  with  the  basal  plane,  and  the  other 
three  form  another  equilateral  triangle  of  edge  a  that  is  coplanar  with  the  basal  plane 
and  is  rotated  60°  from  the  other  triangle.  The  line  connecting  the  centers  of  these 
equilateral  triangles  is  collinear  with  the  c-axis.  Depending  on  the  c/a  ratio,  the 
octahedron  is  deformed  along  the  c-axis  by  varying  the  spacing  between  these  two 
equilateral  triangles.  Since  the  c/a  ratio  for  zinc  oxide  is  close  to  2\/6/3,  the  octahedra 
formed  by  this  lattice  are  nearly  regular  with  edge  a. 


Figure  35.  Unit  cells  of  zinc  oxide  showing  (a)  tetrahedral  and  (b)  octahedral 
interstitial  sites.  The  zinc  nuclei  are  colored  gray  and  the  oxygen  nuclei  are  colored 
red.  Interstitial  sites  surrounded  by  zinc  nuclei  are  colored  in  dark  red,  and  sites 
surrounded  by  oxygen  nuclei  are  colored  in  black. 
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Appendix  C.  Crystal  Structure  of  Titanium  Dioxide 


Titanium  dioxide  is  a  binary  metal  oxide  composed  of  one  part  titanium  and  two  parts 
oxygen  by  atomic  composition.  Anatase,  bookite,  and  rutile  are  the  three  polymorphs 
of  titanium  dioxide  that  are  stable  at  standard  temperature  and  pressure.  Rutile  is 
the  most  common  naturally  occurring  and  synthetically  derived  polymorph  and  is  the 
form  of  titanium  dioxide  studied  in  this  dissertation.  It  has  a  tetragonal  structure 
with  six  basis  sites  as  shown  in  Figure [36}  Tetragonal  structures  are  commonly  defined 
by  three  orthogonal  vectors.  Two  of  the  vectors  (01,02)  have  length  a  and  one  vector 
(c)  has  length  c.  The  plane  formed  by  the  vectors  of  length  a  is  called  the  basal 
plane,  and  the  axis  formed  by  the  vector  of  length  c  is  called  the  c-axis.  For  titanium 
dioxide,  the  value  of  a  is  4.594  A,  and  the  value  of  c  is  2.959  A  with  c/a  —  0.6441 
m  Two  of  the  basis  sites  are  composed  of  titanium  cations,  and  are  located  at  the 
origin  and 


2“I  + 


—do 

2 
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Four  of  the  basis  sites  are  composed  of  oxygen  anions,  and  are  located  at 


ud\  +  ua2, 


(1  -  u)ai  +  (1  -  u)a2, 

1  f1  \  ^  U 

-  +  u  J  ai  +  I  -  —  u  J  a2  +  -c,  and 

i  — MU  +  Q  +  uU  +  ic. 


(94) 


The  value  of  u  is  0.306  m-  The  titanium  cations  are  coordinated  with  six  oxygen 
anions  forming  a  slightly  deformed  regular  octahedron,  and  the  oxygen  anions  are 
coordinated  with  three  titanium  cations  with  all  fonr  atoms  in  a  {110}  plane.  The 
octahedron  is  regular  when  u  =  1  —  y/2/2  «  0.2929  and  c/a  =  2u  ~  0.5858.  The 
octahedron  is  defined  by  fonr  oxygen  anions  forming  a  rectangle  in  a  {110}  plane, 
and  two  oxygen  anions  forming  an  axis  normal  to  the  rectangle  and  intersecting  the 
center.  Since  each  titanium  cation  is  assigned  to  two  basis  sites,  there  are  two  sets 
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of  octahedra  that  differ  by  a  rotation  of  90°  in  the  basal  plane.  The  bond  length 
between  the  titanium  cation  and  the  planar  oxygen  anions  is  1.944  A,  and  between 
the  titanium  cation  and  the  axial  oxygen  anions  is  1.988  A. 


a2 

[01C 


Figure  36.  A  unit  cell  of  titanium  dioxide  with  the  titanium  nuclei  colored  turquoise 
and  the  oxygen  nuclei  colored  red. 

Titanium  dioxide  is  composed  of  tetrahedral  and  octahedral  interstitial  sites  as 
shown  in  Figure  [37}  Tetrahedral  sites  are  located  at 

1_ 

vai  +  va2  +  -c, 

(1  -  v)ai  +  (1  -  v)a2  +  ~c, 


(95) 


and 


where 


(96) 


The  value  of  v  is  0.176.  This  site  is  surrounded  by  four  equidistant  oxygen  anions 
which  form  the  vertices  of  a  slightly  deformed  regular  tetrahedron.  Two  of  the  oxygen 
anions  form  an  axis  parallel  to  the  c-axis  and  the  other  two  form  an  axis  parallel  to 


a  (110)  direction.  The  distance  to  the  oxygen  anions  is  1.703  A  in  rutile.  This  site  is 


also  surrounded  by  three  titanium  cations  which  are  contained  within  a  {110}  plane 
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along  with  the  interstitial  site.  Two  of  the  titanium  cations  form  an  axis  parallel  to 
the  c  axis  and  the  other  is  displaced  from  the  interstitial  site  along  a  (110)  direction. 
The  distance  to  the  c  axis  titanium  cations  is  1.871  A,  and  to  the  other  is  2.104  A. 
Octahedral  sites  are  located  at 

-(ai  +  c),  and 

1  (97) 

-(a2  +  c). 

This  site  is  surrounded  by  six  oxygen  anions  and  six  titanium  cations  with  each  set  of 
atoms  forming  the  vertices  of  a  slightly  deformed  regular  octahedron.  These  octahedra 
are  similar  to  the  octahedra  formed  by  the  nearest-neighbor  oxygen  anions  around  the 
titanium  cations,  in  that  they  are  defined  by  four  atoms  forming  a  rectangle  which  is 
perpendicular  to  an  axis  formed  by  the  other  two  atoms  and  intersecting  the  center. 
The  oxygen  rectangle  forms  a  {2,  1—^,0}  plane,  and  the  titanium  rectangle  forms 
a  {100}  plane.  In  rutile,  the  distance  to  the  axial  oxygen  anions  is  1.665  A,  to  the 
planar  oxygen  anions  is  2.227  A,  to  the  axial  titanium  cations  is  2.297  A,  and  to  the 
planar  titanium  cations  is  2.732  A. 


(a)  (b) 


Figure  37.  Unit  cells  of  titanium  dioxide  showing  (a)  tetrahedral  and  (b)  octahedral 
interstitial  sites.  The  titanium  nuclei  are  colored  turquoise,  the  oxygen  nuclei  are 
colored  red,  and  the  interstitial  sites  are  colored  black. 
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Appendix  D.  Crystal  Structure  of  Stannous  Hexathiohypodiphosphate 

Stannous  hexathiohypodiphosphate  (SPS)  is  a  metal  thiophosphate  composed  of  two 
parts  tin,  two  parts  phosphorus,  and  six  parts  sulfur  by  atomic  composition.  The 
tin  cations  have  an  effective  charge  state  of  2+  while  the  thiophosphate  anions  have 
an  effective  charge  state  of  4—.  The  prototype  for  the  thiophosphate  anion  is  hy- 
podiphosphate.  This  anion  is  composed  of  two  phosphorus  cations  that  share  a  single 
bond.  Each  of  the  phophorus  cations  have  an  effective  charge  state  of  4+.  Each 
phosphorus  cation  also  shares  single  bonds  with  three  chalcogen  anions  that  have 
an  effective  charge  state  of  2—.  In  the  case  of  hexathiohypodiphosphate,  all  of  the 
chalcogen  anions  are  sulfur  anions.  Each  set  of  three  sulfur  anions  roughly  forms  an 
equilateral  triangle.  The  plane  formed  by  each  triangle  is  roughly  perpendicular  to 
the  axis  formed  by  the  two  phosphorus  cations.  This  axis  roughly  intersects  the  center 
of  each  triangle.  The  triangles  are  rotated  from  one  another  by  a  180°  rotation  along 
this  axis.  Each  set  of  three  sulfur  anions  is  at  the  periphery  of  the  hypodiphosphate 
anion  at  roughly  the  same  distance  from  their  nearest  neighbor  phosphorus  cations. 
Each  hexathiohypodiphosphate  anion  is  surrounded  by  six  tin  cations  that  form  a 
highly  distorted  regular  octahedron  as  seen  in  Figure  |38j 
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Figure  38.  The  crystal  structure  of  stannous  hexathiohypodiphosphate,  showing 
two  hexathiohypodiphosphate  anions  and  their  six  nearest-neighbor  tin  cations  as 
viewed  along  the  (a)  a,  (b)  b ,  and  (c)  c  directions.  Phosphorus  cations  are  blue,  tin 
cations  are  green,  and  sulfur  anions  are  red.  The  box  outlined  in  black  is  a  unit  cell. 


The  crystal  structure  is  monoclinic  with  twenty  basis  sites  as  listed  in  Table  [7] 
[72].  Monoclinic  structures  are  commonly  defined  by  three  vectors.  The  vector  b  is 
orthogonal  to  vectors  a  and  c.  Vectors  a  and  c  are  not  orthogonal  with  an  angle  of 
separation  f3.  For  SPS  the  magnitude  of  vector  a  is  9.378  A,  vector  b  is  7.488  A,  and 
vector  c  is  6.513  A.  The  angle  f3  is  91.15°  [72] . 


Table  7.  Basis  for  the  crystal  structure  of  SPS  in  terms  of  the  lattice  vectors  [72] . 
The  additional  10  basis  sites  are  found  by  adding  \  to  a  and  c,  and  multiplying  b  by 
-1. 


Element 

a 

b 

c 

Sn 

0.5270 

0.3856 

0.7224 

Sn 

0.0279 

0.1245 

0.7870 

P 

0.6836 

0.8608 

0.8108 

P 

0.8174 

0.6447 

0.6908 

S 

0.4898 

0.7511 

0.8491 

S 

0.7788 

0.9462 

0.0750 

s 

0.6942 

0.0517 

0.5921 

s 

0.0149 

0.7437 

0.6497 

s 

0.7147 

0.5633 

0.4288 

s 

0.8050 

0.4480 

0.9025 
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Appendix  E.  EPR  Spectroscopy  Data  of  the  Nitrogen  Acceptor  in  ZnO 

The  angle  theta  is  from  the  c-axis  to  the  the  [1010]  direction.  The  angle  phi  is  from 
the  [1010]  direction  to  the  [21 10]  direction.  Calculations  use  a  Plank’s  constant  value 
of  6.62606957  x  10-34  Js  [32],  a  Bohr  magneton  value  of  9.27400968  x  10-24  J/T  [32], 
and  a  g  value  for  the  MgO:Cr  signal  of  1.9800  [38].  The  corrected  data  in  bold  in  the 
basal  zinc  table  were  corrected  using  the  MgO:Cr  calibration  data  in  the  axial  zinc 
table,  since  these  data  were  taken  with  the  same  experimental  setup. 
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9.402031E+009HZ 

Field(MgO:Cr) 

0.3392931T 

Correction  Factor 

0.999930711 

Defect 
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Frequency 

(GHz) 

Field  (G) 

Corrected 
Field  (mT) 

Theta 

(Degrees) 

Phi  (Degrees) 
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90 

0 
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340.060 

90 

0 
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3404.281 

340.405 

90 

0 

9.394845 

3401.426 

340.119 

90 

240,120 

9.394845 

3404.701 

340.447 

90 

240,120 

9.394976 

3397.609 

339.737 

90 

5 

9.394976 

3400.854 

340.062 

90 

5,125 

9.394976 

3404.188 

340.395 

90 

5,125 

9.394976 

3401.933 

340.170 

90 

245 

9.394976 

3404.981 

340.475 

90 

245 

9.394853 

3397.617 

339.738 

90 

10 

9.394853 

3400.918 

340.068 

90 

10 

9.394853 

3404.330 

340.409 

90 

10 

9.394853 

3402.200 

340.196 

90 

250 

9.394853 

3405.535 

340.530 

90 

250 

9.394853 

3400.265 

340.003 

90 

130 

9.394853 

3403.512 

340.328 

90 

130 

9.394948 

3397.867 

339.763 

90 

15 

9.394948 

3401.140 

340.090 

90 

15 

9.394948 

3404.542 

340.431 

90 

15 

9.394948 

3402.576 

340.234 

90 

255 

9.394948 

3405.755 

340.552 

90 

255 

9.394948 

3399.860 

339.962 

90 

135 

9.394948 

3403.135 

340.290 

90 

135 

9.394895 

3398.073 

339.784 

90 

20 

9.394895 

3401.349 

340.111 

90 

20 

9.394895 

3402.685 

340.245 

90 

260,140 

9.394895 

3406.029 

340.579 

90 

260 

9.394895 

3399.413 

339.918 

90 

140 

9.394876 

3398.419 

339.818 

90 

25 

9.394876 

3401.655 

340.142 

90 

25 

9.394876 

3402.962 

340.273 

90 

265 

9.394876 

3398.975 

339.874 

90 

145 

9.394876 

3402.277 

340.204 

90 

145 

9.394868 

3398.641 

339.841 

90 

30,150 

9.394868 

3401.885 

340.165 

90 

30,150 

9.394868 

3402.866 

340.263 

90 

270 

9.397276 

3320.726 

331.990 

0 

0 

9.397276 

3327.254 

332.643 

0 

0 

9.397276 

3333.549 

333.272 

0 

0 
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Frequency(MgO:Cr)  9.397601E+009  Hz 


Field(MgO:Cr)  0.3391937T 

Correction  Factor  0.999752459 

Defect  Axial  Zinc 


Frequency 

(GHz) 

Field  (G) 

Corrected 
Field  (mT) 

Theta 

(Degrees) 

Phi  (Degrees) 
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3304.348 

330.353 

0 

0 

9.397276 

3317.319 

331.650 

0 

0 

9.397276 

3330.418 

332.959 

0 

0 

9.397301 

3306.702 

330.588 

10 

0 

9.397301 

3319.553 

331.873 

10 

0 

9.397301 

3332.499 

333.167 

10 

0 

9.397409 

3313.662 

331.284 

20 

0 

9.397409 

3326.092 

332.527 

20 

0 

9.397409 

3338.522 

333.770 

20 

0 

9.397475 

3324.812 

332.399 

30 

0 

9.397475 

3336.970 

30 

0 

9.397475 

3347.962 

30 

0 

9.397504 

3339.088 

333.826 

40 

0 

9.397504 

3349.683 

334.885 

40 

0 
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Frequency(MgO:Cr) 

9.402031E+009  Hz 

Field(MgO:Cr) 

0.3392931T 

Correction  Factor 

0.999930711 

Defect 

Nitrogen 

Frequency 

(GHz) 

Field  (G) 

Corrected 
Field  (mT) 

Theta 

(Degrees) 

Phi  (Degrees) 

9.394845 

3415.959 

341.572 

90 

0 

9.394845 

3419.242 

341.901 

90 

0 

9.394845 

3422.672 

342.243 

90 

0 

9.394976 

3416.012 

341.578 

90 

5 

9.394976 

3419.286 

341.905 

90 

5 

9.394976 

3422.729 

342.249 

90 

5 

9.394853 

3415.983 

341.575 

90 

10 

9.394853 

3419.221 

341.898 

90 

10 

9.394853 

3422.660 

342.242 

90 

10 

9.394948 

3416.002 

341.577 

90 

15 

9.394948 

3419.280 

341.904 

90 

15 

9.394948 

3422.678 

342.244 

90 

15 

9.394895 

3415.980 

341.574 

90 

20 

9.394895 

3419.259 

341.902 

90 

20 

9.394895 

3422.661 

342.242 

90 

20 

9.394876 

3415.965 

341.573 

90 

25 

9.394876 

3419.240 

341.900 

90 

25 

9.394876 

3422.674 

342.244 

90 

25 

9.394868 

3415.937 

341.570 

90 

30 

9.394868 

3419.199 

341.896 

90 

30 

9.394868 

3422.670 

342.243 

90 

30 
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9.397601E+009  Hz 
0.3391937T 


Frequency(MgO:Cr) 

Field(MgO:Cr) 

Correction  Factor  0.999752459 

Defect _ Nitrogen 


Frequency 

(GHz) 

Field  (G) 

Corrected 
Field  (mT) 

Theta 

(Degrees) 

Phi  (Degrees) 

9.397276 

3336.950 

333.612 

0 

0 

9.397276 

3366.036 

336.520 

0 

0 

9.397276 

3395.089 

339.425 

0 

0 

9.397301 

3338.913 

333.809 

10 

0 

9.397301 

3367.606 

336.677 

10 

0 

9.397301 

3396.271 

339.543 

10 

0 

9.397409 

3344.702 

334.387 

20 

0 

9.397409 

3372.185 

337.135 

20 

0 

9.397409 

3399.672 

339.883 

20 

0 

9.397475 

3353.746 

335.292 

30 

0 

9.397475 

3379.272 

337.844 

30 

0 

9.397475 

3404.691 

340.385 

30 

0 

9.397504 

3365.129 

336.430 

40 

0 

9.397504 

3387.927 

338.709 

40 

0 

9.397504 

3410.726 

340.988 

40 

0 
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Appendix  F.  EPR  Spectroscopy  Data  of  the  Silicon  Interstitial  in 


TiOo 

Frequency  data  is  in  GHz  unless  otherwise  stated.  The  data  headings  label  each  line 
with  their  assigned  position  in  the  spectra.  “C”  stands  for  center,  “H”  stands  for 
high  held,  “L”  stands  for  low  held,  “I”  stands  for  inner  hyperhne,  and  “O”  stands 
for  outer  hyperhne.  The  data  in  bold  in  the  (001)  data  were  measured  with  the  Hall 
probe.  The  magnetic  held  measured  for  the  MgO:Cr  signal  with  the  Hall  probe  is 
3385.69  G.  Calculations  use  a  Plank’s  constant  value  of  6.62606957  x  10-34  Js  [32] ,  a 
Bohr  magneton  value  of  9.27400968  x  10~24  J/T  [32],  and  a  g  value  for  the  MgO:Cr 
signal  of  1.9800  [38]. 
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30 
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50 
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60 
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3478.785 
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3519.829 

3532.646 

3560.227 

80 
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3505.710 
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3532.899 

3561.630 
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3476.059 

3505.621 

3519.541 

3533.363 
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LO 

LI 

C 

HI 

HO 

[001] 
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9.562257 

9.562257 

9.562257 
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9.562017 

9.562017 

9.562017 

20 
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9.561770 
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9.561770 

30 
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9.561273 

9.561273 

9.561273 

9.561273 

40 

9.561737 

9.561737 

9.561737 

9.561737 

9.561737 

50 

9.561376 

9.561376 

9.561376 

9.561376 

9.561376 

60 

9.561210 

9.561210 

9.561210 

9.561210 

9.561210 

70 

9.562017 

9.562017 

9.562017 

9.562017 

9.562017 

80 

9.562554 

9.562554 

9.562554 

9.562554 

9.562554 

[100] 

9.563654 

9.563654 

9.563654 

9.563654 

9.563654 

Corrected  (mT) 

LO 

LI 

C 

HI 

HO 

[001] 

350.804 

352.556 

354.253 

10 

350.590 

352.537 

354.409 

20 

350.168 

352.475 

354.698 

30 

349.653 

352.360 

354.988 

40 

349.150 

352.268 

355.307 

50 

348.628 

352.136 

355.560 

60 

348.166 

350.833 

352.010 

353.171 

355.784 

70 

347.840 

350.652 

351.944 

353.225 

355.983 

80 

347.621 

350.532 

351.899 

353.251 

356.124 

[100] 

347.567 

350.523 

351.915 

353.297 

356.210 
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9.402333E+009  Hz 
3.393154E-001T 
0.999897112 


Frequency(MgO:Cr) 

Field(MgO:Cr) 


Correction  Factor 
Plane _ (110) 


NMR  Data  (G) 

LLO 

LU 

L 

LHI 

LHO 

HLO 

H 

HHO 

[001] 

3514.267 

3531.744 

3548.616 

3514.267 

3531.744 

3548.616 

10 

3508.647 
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3515.604 
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20 
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3550.926 
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30 
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40 
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3528.433 
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50 
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3495.640 
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60 
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70 
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80 
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3494.039 

3531.936 
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3454.603 
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50 
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60 
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70 
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LU 
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H 
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355.053 

351.524 

353.266 
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20 
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352.393 

355.056 

353.499 

355.181 

30 
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354.009 

355.685 

40 
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349.433 
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351.940 

352.807 

354.529 

356.190 

50 
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348.024 
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351.017 
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354.963 
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60 
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347.014 
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353.751 

353.884 
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357.218 

70 
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346.140 
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353.409 

354.409 

356.095 

357.712 

80 

341.805 

345.583 

347.480 

349.368 

353.157 

354.636 

356.412 

358.019 

[110] 

341.591 

345.425 

347.337 

349.249 

353.089 

354.874 

356.546 

358.153 
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9.575892 
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Appendix  G.  EPR  Spectroscopy  Data  of  the  Sulfur  Vacancy  in 


SrioPoSc 

The  angle  theta  is  from  the  c  direction  towards  the  a  direction.  The  angle  phi  is  from 
the  a  direction  towards  the  b  direction.  The  columns  labeled  with  a  g  are  the  center  of 
the  spectrum.  The  columns  labeled  with  a  P  and  a  number  are  phosphorus  hyperfine 
lines  where  a  lower  number  indicates  a  lower  field.  The  columns  labeled  with  Sn  and 
a  number  are  the  center  of  a  set  of  lines  associated  with  a  tin  hyperfine  interaction 
where  a  lower  number  indicates  a  lower  field.  The  columns  labeled  with  Sn,  a  number, 
P,  and  another  number  are  phosphorus  hyperfine  lines  within  a  set  of  hyperfine  lines 
associated  with  a  tin  hyperfine  interaction  where  lower  numbers  indicate  a  lower 
field.  Calculations  use  a  Plank’s  constant  value  of  6.62606957  x  10~34  Js  [22],  a  Bohr 
magneton  value  of  9.27400968  x  10-24  J/T  [32],  and  a  g  value  for  the  MgO:Cr  signal 
of  1.9800  (38]. 


109 


MgO:Cr  Frequency  9.39691500E+09  Hz 

MgO:Cr  Field  3.39152000E-01  T 

Correction  Factor  0.999802 

Lowest  Field  Phosphorus  Hyperfine  Line 


Phi 

Theta 

Frequency 

Field  (G) 

Field  (mT) 

90 

30 

9.389342 

3520.415 
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90 

40 
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50 
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90 

60 
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3525.815 
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MgO:Cr  Frequency  9.39553800E+09  Hz 

MgO:Cr  Field  3.39106400E-01  T 

Correction  Factor  0.999790 
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